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Die Grundgleichungen fiir den Verbund bei Stahlbeton-Rechteckplatten 
Von R. Trostel 


1. Einleitung. Die gréf8eren Schwierigkeiten, die bei der Untersuchung des Plattenverbundes 
gegentiber dem Problem des Tragerverbundes auftreten, sind im wesentlichen dieselben, die die 
Theorie der elastischen Platten von der gewéhnlichen Balkentheorie unterscheiden. Die grund- 
satzliche statische Unbestimmtheit des Plattenproblems zwingt auch im vorliegenden Falle, den 
‘Lésungsweg tiber Verschiebungsgleichungen zu suchen!. Dabei gehen wir ebenso wie in der Theorie 
elastischer Platten vor, wo wir vermége geeigneter Verformungsannahmen letztlich auf Verschie- 
bungsgleichungen fiir die Mittelflachenverformungen gefiihrt werden. Erst nach Kenntnis des Ver- 
schiebungszustandes kénnen wir die am Verbundelement wirkenden Gesamtschnittlasten berech- 
nen, die wir nun noch im Sinne von Satiler? fiir eine Spannungsberechnung auf die jeweils auf 
den Beton- bzw. Stahlquerschnitt entfallenden Anteile (VerteilungsgréBen, UmlagerungsgréBen) 
aufzuschliisseln haben. 

In den nachfolgenden Untersuchungen, die sich in Abschnitt 3 auf den elastischen Anfangszu- 
stand und in Abschnitt 4 auf zeitlich beliebig verinderliche Vorginge mit Bericksichtigung des 
Maxwellschen Materialverhaltens des Betons beziehen, werden zundchst die am Verbundelement 
wirkenden Gesamtschnittlasten bzw. die Teilschnittlasten in Abhangigkeit der Plattenver- 
zerrungen (Dehnungen u. Kriimmungen der Mittelflache bzw. geeigneter Schwerebenen) aus- 
gedriickt. Nach Elimination der Verzerrungen folgen damit Beziehungen zwischen den Teil- 
schnittlasten und den Gesamtschnittlasten. Fir die Bestimmung der Gesamtschnittlasten stehen 
die Gleichgewichtsbedingungen am Plattenelement zur Verfiigung, aus denen wir dann letztlich 
die Verschiebungsgleichungen des Problems folgern kénnen. Um die Rechnung nicht unnétig zu 
komplizieren und die Anzahl der ideellen Querschnittskonstanten méglichst klein zu halten, wird 


die fiir Stahlbeton ohnehin kleine Querkontraktionszahl vernachlassigt. Als Beispiel schlieBt in 


Abschnitt 5 u. a. eine Untersuchung des Dischinger-Effektes bei Stahlbetonplatten an. 

Die grundsatzlich mégliche Reduktion des Problems auf drei Verschiebungsgleichungen weist 
den Weg zur Untersuchung auch anderer Flachentragwerke (Schalen). Auch das Problem vor- 
gespannter Flachentragwerke laBt sich in dieser Weise behandeln. 


a. 


Abb. 1. 


s _ 2. Erklarung der verwendeten Bezeichnungen. Es bedeuten (Abb. 1) 
baw. F,, die in einer Schnittflache « = konst. bzw. y = konst. (Schnittlange 4x = Ay = 1) 
freigelegten Bewehrungsquerschnittflachen, 


we Die im Spezialfall der rotationssymmetrisch belasteten Kreisplatte vom Verfasser, Bautechn. 36/7 (1959) 
Ss 263 R diragebetie Reduktion des Problems auf eine unbekannte Kraftgré8e benutzt zuvor ebenfalls Verfor- 
al 


hmen. ; : _ 
Fe Sauter, Theotio der Verbundkonstruktionen (Spannbeton, Stahltrager in Verband mit Beton), Berlin 
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F,~ baw. F,y die je Schnittlangeneinheit freigelegten Betonquerschnittflachen, wobei Fi + Fox 
=F, = F,, + hg = Fy 1 = d sind (d = Plattendicke), 

F;,, baw. F;y die je Schnittlangeneinheit anfallenden ideellen Gesamtquerschnittflachen, bestehend 
aus den Betonquerschnitt- und den n-fachen Stahlquerschnittflachen (n = E,/E,), 
also F;, = Fy, +n Fy =d+(n—1) Fx, Fiy = Foy +n Fy =d+(n—]1) Fry, 

©,. baw. G,y die Schwerpunkte der Betonflachen F,, baw. Fy y, 

6.x baw. G.y die Schwerpunkte der Stahlquerschnittflaichen F.,, bzw. F.y, 

Gix baw. Gy die Schwerpunkte der ideellen Gesamtquerschnittflachen F;, bzw. F;,, 

Jy baw. Joy die (Eigen-)Tragheitsmonente der Betonflachen F,,, bzw. F,, hinsichtlich ihrer eige- 


nen Schwerachsen ©) bzw. Gy, 2. B. Jie a es a! 
Fh x 


Jox baw. te y die Eigentragheitsmomente der Stahlquerschnittflachen F,, bzw. F., hinsichtlich 
ihrer eigenen Schwerachsen ©,x bzw. Gey, zB. Jog = f 2x Mes 
Fox 
Jix baw. Jiy die Tragheitsmomente der ideellen Gesamtquerschnittflachen F;, baw. F;5 hin- 
sichtlich ihrer ideellen Schwerachsen ©;, bzw. ©;,. Es sind 
Ties, 9) = Tope,y) + 0 Tete, y) + Sbee,9) Fo(esy) +m Sete, 9) Fete, y) « 
Dabei bedeuten 
Spx baw. spy die Abstande der Betonquerschnitts-Schwerpunkte ©,, bzw. Gp), von den ideellen 


Gesamtquerschnittsschwerpunkten ©;, bzw. ©;,, 
Sex baw. s,.y die entsprechenden Abstainde der Stahlquerschnittsschwerpunkte. 


Weiterhin sind 


s, bzw. s, die Abstande der ideellen Gesamtquerschnittschwerpunkte ©;, bzw. ©;, von der 
Plattenmittelflache z = 0, 
z, bzw.z, plattennormale Ordinaten, gemessen von den jeweiligen Gesamtquerschnittschwer- 


punkten ©;, bzw. Giy aus, also z, = z— s,, Zy = %— Sy, 
yn, baw.7,  plattennormale Ordinaten, gemessen von den jeweiligen Beton- bzw. Stahlquerschnitt- 
flachenschwerpunkten. 


Es gelten noch folgende Identitaten: 


In, = fm, 4% = Sn, t= Jn,, Ff, =, 
(Fx) bg rs 5 eee) wy 
[a4,= fudiin f sdf,<0, (a) 
(Fx) Fh x ex 
jy ef, = —-f zdf, +n f-2 df, =0, 
Cee ri pes hate 
gn etre el ee Biss Ree re 
also 
Sheng = Soy ens Sb F, =NS, sa 
und fee ae (b) 
s, d = — l 2, df, =— f 4 Uf, — f %_ Ue = s,, Fug —Sex Fox =(n—1)5,, F, 
(x) (Fix) Fiz) 6 6 ( ) x 
sowie _ 8, d= (n—1) 5, Fy. ¥ 


3. Lastverteilung zur Zeit £ = 0; Ermittlung der Verteilungsgréfen. Wir betrachten hierfiir das _ 
in Abb. 2 skizzierte Plattenelement mit den Kantenlangen 1, das die positiven Koordinaten- und 
Spannungsrichtungen angibt. Die Plattenschnittlasten, worunter wir die > Spannungsresultierenden* 
je Schnittlangeneinheit verstehen, ergeben sich dann zu 


Ue Ox d, x? Ss pea x d Pa) = 5 = 
are Ore a Re yr aay. fe 
Or. [tes dfs G2" fede. m| 
(Fx) ¥ (Fy) s y (2) 


Os fo, 2 df, , Moxy = las z df, , ™M yx = Tyee Oo > Moy = df. 33 
# ey y eye if, Oy eee if, (3a) 


SCANATIEN nm 
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Die Plattenmomente beziehen sich hierbei jeweils auf eine in der Plattenmittelflache liegende Achse 
was der Index 0 zum Ausdruck bringen soll. Als rechentechnisch einfacher erweist es sich aber auch 
hier, die Momente auf die jeweilige, durch den zugehorigen ideellen Gesamtschwerpunkt ©,, bzw. 
G;, verlaufende Horizontalachse zu beziehen, also die Plattenmomente in der Form a 


BAe sing: —— f0, 2d], . M; xy = [%< Baap. « Tg rea SOL sg) <I ae z.. 3b 

,) ee oo eee 
zu definieren. Die nach (1), (2) und (3) definierten Schnittlasten sind in ihrer jeweiligen positiven 
Wirkungsweise in Abb. 3 eingetragen. 


ie 

1 Om, } 
Peet ee 
Gy ardy 


Abb. 3. 


Grundlage fiir die Ermittlung des Zusammenhanges zwischen den Schnittlasten und den elasti- 
schen Verformungen der Platte sind die Hookeschen Gesetze sowie die verallgemeinerte Bernoullische 
Hypothese, die den Verschiebungszustand der Plattenpunkte letztlich 
allein auf die gemeinsamen lotrechten Verschiebungen (Biegeflache) Bee OS: 
und die waagerechten Verschiebungen der Mittelflache reduziert. Wir N ip on 
setzen fur die Komponenten des Verschiebungszustandes mit den lot- | id 
rechten Verschiebungen w(x, y) 


: : ow 
in x Richtung: U = Uy — 257» 
5 : Ow 
in y Richtung: pf a Bie ae re 


Hierbei stellt der zweite Anteil (Abb. 4) jeweils den Anteil der Biegedeformationen dar, wahrend 
U(x, vy) baw. v(x.) die Verschiebungen der Plattenmittelflache in x- bzw. y-Richtung bedeuten. 
Dann sind ~~ 


__ ou wo) ty Boe er: AS 
&s 9x Ox Tigi gti, mh Se 
dv OU o?w 4 
= = == == —ZHx 
Ey Oy dy z By? E0y y? ( ) 


du dv dug , O% o?w ae 
= — ——_ => —— Un 
Le By -- a By +. ae Ee aeey Voxy xy? 


und aus den Hookeschen Gesetzen erhalten wir bei Vernachlassigung der Querkontraktionen die 
Normal- und Schubspannungen 


E E 
6, = Ee, =E (&021—2X_) » 0, =Ee, =E(eoy —2%,) » Try = Vey = gy (Youy — 22% ay)» (5) 


wobei fiir die Betonspannungen E = E, und im Stahlbereich E = E, zu setzen ist. Diese Ansatze 
laufen bei den Normalspannungen im Betonbereich auf die Vernachlassigung der Beton-Quer- 
kontraktion hinaus; im Stahlbereich sind sie, sofern wir die Bewehrungslagen F,,, und F,y als nicht 
miteinander fest verbunden ansehen, sogar richtig. Bei den Schubspannungen enthalt der gewahlte 
Ansatz im Stahlbereich auch die Vernachlassigung der Stahlquerkontraktionszahl, was wir aber im 
Hinblick auf die damit verbundene geringere Anzahl ideeller Gesamtquerschnittswerte in Kauf 


nehmen wollen. 
1 * 
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Nach Einsetzen von (4) in (5) erhalten wir nach Integration gema8 (1) bis (3) unter Beachtung 
der in Abschnitt 2 angefiihrten Beziehungen die Schnittlast-Verzerrungsrelationen 


n,= fo,4f,= f 0,4f.+ J 9, af. =F, | } eos Se eh eles — 22) Oe 
(Fr) ) (Fe x) ) 


(Fox ("ox ex 
ad E0x (Fy, + n F, x) 4 jzdfs-n J =dfa)) = Bs F;, (€0.x — Sq %x) (6a) | 
(Fox (Fe x) 
und entsprechend? 
Ty Lege Bie [32 — 5s, te) ; 3 (6b) 
Ey 128 LA (6c) 
BR NTS y “xy}? 

ny = E, Fiy (&0y — Sy %y) » (6d) | 


sowie 


Mo x = fo,2 Uf, = J o,% df, + ado S (co. —2%,) 2 df, +n J (con — 2 %4) # Ue 
(Fx) (Fox) F (ox) F 


ex ex 


bx ex (Fox) 


= feos ( f zdf,-+n I = dhe) — He fo df Ea) aa zara ~Je0= 7) (7a) : 


und entsprechend 


™0 x y = By [Fixes 822 — Soon ry); (7b) | 
V0 « ] 

mM yx = By (Foy 9,2” — Seay tes) (7c) | 
mo y = b (Fiy Sy E07 — Jidy Xy) - (7d) i 


Hierbei sind Jip x bzw. Jioy die Tragheitsmomente der ideellen Gesamtquerschnitte F;, baw. F; 

-hinsichtlich der durch die Plattenmittelebene verlaufenden Bezugsachsen. Besonders einfach |} 
werden die Zusammenhange zwischen den Momenten und den Biegeverzerrungen, wenn man die 
Momente auf die durch die ideellen Gesamtquerschnitts-Schwerpunkte ©;, bzw. iy verlaufenden 
Achsen bezieht. Man erhalt 


Ms ~¢ = fo,%, &, = J 6, (@— s,) df, = fo,zdf, —-s, fo, @, = Mox, —— S$, N, 
(Fx) (Fx) (Fx) 


= — E, (Jioz — 52 Fiz) 4 = — B, Sint : (8a) |} 
und entsprechend 


Ms xy = Ey, Jiztey > Ms yx = — BE, Jiy zy» Msy = — EB, Jiy x, » (8b, c, d) 


wobei die J;, bzw. Jiy die Tragheitsmomente der ideellen Gesamtquerschnittsflichen hinsichtlich 
ihrer ideellen Schwerachsen ©;, bzw. Giy bedeuten. Auch die Membranschnittlasten kénnen wir in 
einfacher Weise durch die Membranverzerrungen in den dieideellen Gesamtquerschnittschwerpunkte 
enthaltenden horizontalen Plattenebenen ausdriicken. Aus (4) erhalten wir 


ee” = bxl = oy ay ete Pare ayo = Eylens, = €0y — Sy %y » ( ) 
“ 9 
VED = Yeyle = oy =Yoxy Sy Uxy » yey = Yayle = sy Slay ae Sy Yxy 
und damit aus (6) 
(x cs S E s s 
ng = Ey Fix ee”, he ai in Vo), Bey = = a iy VED » n, = E, Fy, (10 a—d) 


Fur die Ermittlung der VerteilungsgréBen driicken wir nun noch die auf die reinen Beton- bzw. 
Stahl-Teilquerschnittflachen wirkenden Schnittlasten durch die Membran- und -Biegeverzerrungen 
der ideellen Gesamtquerschnitte aus. Es ist 


z=, — (Sox — $,) = re (Soy — Sy) =). ol (Sex 1 Sz) = Ney oe (Sey a8 Sy) ? (11) 


_, | Die zugeordneten Membranscherkrafte sind also strenggenommen nicht gleichgro8. Siehe hierzu die Aus- 
fiihrungen am Ende dieses Abschnittes. 
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und damit nach (4) und (9) 
Es — 605 —— [n, = (Sb oe Sz) | TE CO om nn Hg Vip. ——- Sb x) ty = glx) 5 Che ——= Sb x) Hy 


und entsprechend 
(sy) 


by = by eit Si 9) ey ee \—_(n, + Sey) ty (12b) 
we = yes) —2 Wi. — Sy PS) Pie ye —_ ; (ee + a Uxy > (12 c) 
= yep —2 (n, , — Spy) Xxy = yy —2 (Ney + Sey) Hay - (12d) 


Mit den Hookeschen Gesetzen (5) erhalten wir damit fiir die auf die Teilquerschnitte wir- 
kenden Membranschnittlasten und die Momente, jeweils bezogen auf die Teilquerschnittschwer- 
punkte unter Beachtung der Relationen (a) und (6) schlieBlich nach Integration, 


J df, = By Fou ee? + sox)» Teg Af GO, Of, = BPog (ee — $e %z) sen) 
bx 


Bie oy Uf, = Ey Fry (26? + s5y xy) » Ney = fo, fp E Foy (es — Sey % y) 
(Foy) (Fey) 


(8.2) 
tay =f toy fe By Foa( 
(6 x) 


(Sx) 
Y, 
+ S1a%). Nexy = Jf tarde = E, Fa ee su) ? 
Fex 


hats ves my 
Ny y « = Ty x y= Ey BE + siy ter) Ne yx = freeads. = E. F, fe 
(Foy) Fe ) g 


~ saya] 
(13) 
Sas 04%, , Ux = — Ey Jon re» sen SS JE ORUIES df= — E, Teghes : 
bx) ex, 

aa fo Bo Nyy df, = — E, Joy xy» Mey = f265 1, df, = — E, Soy tte 

( ¥by) (Fey) 
be yi is Txy Ny x df, = — Ee Tox Axy » Me ey f Txy Ne x df, =—E e Seeteae 

(Fox) (Fe 
Mp yx = J ty Nyy df, = — Ey Joy Xxy » Meyx = J tyeNey df, = — Bi; Jace 

(Foy) ; Fey) 


wobei die durch einen Querstrich gekennzeichneten Tragheitsmomente die Eigentragheitsmomente 
des Beton baw. Stahlteiles, jeweils bezogen auf die zugehérige Teilschwerachse bedeuten. Driickt 
man hierin die Verzerrungsgréfen ¢ und x nach (8) bzw. (10) durch die resultierenden Schnittlasten 
der Gesamtquerschnittflachen aus, so hat man die VerteilungsgréBen zur Zeit t = 0 in Ab- 
hangigkeit von den Gesamtschnittlasten ermittelt: 


= Foe Foz Sta . = wie ae Bae Ne 
bx =F a $x 9 CF paar % $x 9 
Lee Tix FF, Tix 
see Foe Soe t Fa Faas AR 
Nb xy =F May Zz Ms x y 2 Nexy—=— a Nxy Zit — Msx y 9 
é Ce 7 Six Bs if Bis Tix 
F es Sb y ek, nF oy Sey 
Nb yx = yea Ms y x 9 Le aaa My x i syn? 
F J 
4 i : (14) 
a ee Foy Soy ‘ nF, A mB ey Sey 
(OG ae ae wee He Ly esd rusts 
KF, Jiy a J; 
dpe Ee me J, me e Joy’ 
ii —— mm => b a ? by — ? 
bx De sx? bxy Tie sxy 2 yx J. syx y¥ Jiy sy 
Ws Os ae nJ, nde sa 
Mex = Ts MN, x 9 Mex y = Bie ee? Me yx StI cesta Mey = Jiy sy* 


Die Gesamtschnittlasten der Platte werden dabei nach der Theorie der orthotropen Platten | 
bzw. Scheiben berechnet. Die maBgebenden Gleichungen gewinnen wir, indem wir in die Gleich- 
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gewichtsbedingungen fir das Plattenelement* (Abb. 3, hier wurden der Ubersichtlichkeit 
wegen nur die Zuwachse eingetragen) 


on on 
SU ee ge afee td =0, 15 
SK, =0: = oT hee (15a) 
on. on, 
Laas xy af vey b 
2K, = 0: aging ay 1 Py Y Cee) 
oq oq, 
ey Come . Rec ae, =— 15 
» K, =0: eee (5c) 
SN p. == Oe omg xy | ony —q,=—0 (15d) 
= Ax . ax i] dy dy ? 
om (a) 
SMG gs eh Oe cong eee (15e) 


Ox oy 


die Schnittlast-Verzerrungsrelationen (6) bzw. (7) einsetzen, womit wir letztlich drei gekoppelte | 
Differentialgleichungen fiir die Verschiebungen der Mittelflache erhalten. Zunachst hat man aus 
(15d, e) 


am ». omy 7% j omy zy dmg y 
de ~~ Gx eect dy ~ bx ay” (ita 
so daf nach EHinsetzen in (15c) 
am , 2 em, 
y 
pat Tb te py Monet Moye) tet Pa 8 (16¢) 


entsteht. Mit (7) wird daraus 


e2 e2 o2 
gaz Wio x %x] Te by Se 2 hoy) Heyl aye io y Hy] 


ac Pehue 1 Fog) 22 |G, Fie ee) (16d) 
oe oye Ux Vix Ox dx Oy x fix y * iy 2 ay? 'y iy €0y] » 


und insbesondere bei konstanten Steifigkeiten 


o*w o*w o*w 
Jive aa (Seon  Je09) bx By? Sees ay 
No Mee Ee 6 @ (Au, a 8, Fy a @ [Ou av 
Sige F | ding + Bs (Coan Fe by (A alee > yf ee 
Hinsetzen der Gleichungen (6) in die beiden ersten Gleichgewichtsbedingungen liefert 
Ox (Fix £0 x) aa ay (a te) =z Dy (s,, F;., Hx) > by. (sy ee Y%zy) + E, => 0 Py (17a) 
0 Yoxy é a) (7) 
Fe (Fe 82) + 2 iy 005) — Bee, Fetes) — 2 (6, Fay) + B= 0 (17b) 
und insbesondere bei konstanten Steifigkeiten 
_ Bug Fy @ Ou Ovo Pp é Ow ow 
Fis Ox? tt 2 oe 1 ; ’ ess Is Fis ag + Sy ee =20,, (17e) 
_ 2 vy Fi, a OUg Bvy Py 7) o*w ow 
ot alg tS) +R— ema tym <0. ara) 


Wir kommen nun auf die Momentengleichgewichtsbedingung hinsichtlich der Normalachse z 
und damit auf die Frage der Gleichheit der zugeordneten Membranscherkrafte Nyy baw. ny, zuriick, 
Zur Klarung dieses Sachverhaltes ist es notwendig, die Schubverformungen eines Vetbundelementes 
der Dicke dz (bestehend aus einem Beton- und einem Rundstahlanteil, der angenahert als Recht- 
eckquerschnitt angesehen werden kann) einer genaueren Betrachtung zu unterziehen 2. In Abb. 5 


' Die Momentengleichgewichtsbedingung hinsichtli ; 2 ; 
Ae fag zleicl gung hinsichtlich der z-Achse lassen wir zunachst unberiicksichtigt. 
Mit dieser Schnittfiihrung folgen wir der wblichen Betrachtungsweise, die in der Theorie der diinnen Biateon 


angewendet wird, indem man zur Herleitune der Schni I i 
( g der Schnittlast-Verzerrungsrelationen die Spannungen o. 
Txy als einen ebenen Spannungszustand bildend ansicht. s re ee see 
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ist ein in «-Richtung bewehrtes Verbundelement und dessen Schubverzerrung y = yy + y, skizziert. 
-Wahrend die im Betonteil entstehenden Schubspannungen auf Grund der konstanten Schubver- 
zerrung y ebenfalls konstant sind und 


E 
%=Gyr~ Sy (18a) 


betragen, miissen wir im Stahlbereich eine veranderliche Schubverzerrung und damit auch einen 
veranderlichen Schubspannungsverlauf voraussetzen. Dies zeigt folgende Uberlegung. 


‘s 
| dbgrhe 
as 

a 
Bb rey 
| 
| | 
L 

lz 
mes 


Da die Betonspannungen 1, langs der Grenzflichen AA’ und BB’ stetig in den Schubspannungs- 
verlauf t.y im Eisen tibergehen miissen (t, 4 = T,p = 1%), andererseits jedoch der Schubmodul G, 
des Stahles wesentlich gréBer als derjenige des Betons ist, sind wegen 


TeA TeB Tv ob = 
Vet eB a G, rhe G, Se G, =< Ce SY. 


die Schubverzerrungen an den Randern des Eisens kleiner als die benachbarten (konstanten) Be- 
tonverzerrungen y. Damit jedoch der Verschiebungszustand beim Ubergang vom Beton zum Stahl- 
teil langs der Grenzflachen AA’ und BB’ stetig ist, miissen die geringeren Randschubverzerrungen 
des Eisens durch groBere Schubverzerrungen im Innern des Stahlbereiches ausgeglichen werden. 
Die gestrichelt in Abb. 5 eingetragene konstante Schubverzerrung ist also im Bereich des Eisens 
nicht vertraglich. Die tatsaichlichen Verzerrungen und damit auch die Stahlschubspannungen 
werden eher den in Abb. 5 durch die ausgezogenen Linien markierten Verlauf haben. Ein mit y ver- 
anderlicher Schubspannungsverlauf tT, , , fordert jedoch aus Gleichgewichtsgriinden am Element be- 
kanntlich die gleichzeitige Wirkung von Normalspannungen ©, und 0,, wobei letztere an den 
Randern AA’ bzw. BB’ verschwinden miissen, da der angrenzende Betonkérper zufolge seiner 
vorausgesetzten Deformationsfigur allein unter der Schubbelastung 1, stehen kann. Hier liegt also 
ein zweiachsiger Spannungszustand mit den Randbedingungen 


1) Rand 4a'ly = 2 


Texy =%> Gy=0, U=Uy—WZ> also Beg t= Ose buwi nO One te US Oy 1) XS 

; eds): 

2) Rand BB (y= 4): 
a h = 

ayy. Oy = 9,4 u=Uuy t+N>? also ¢, =0. bzw. 6,=0,. v= vy +yo% 


vor, und die Aufgabe besteht darin, die an den Randern x = 0 bzw. x = Ax aufzubringenden 
Normal- und Schubspannungen G,, und Tex zu bestimmen, die die genannten Randbedingungen 
befriedigen. Die strenge Lésung dieser biharmonischen Randwertaufgabe ist jedoch fur die Klarung 
der Frage nach der Gleichheit der zugeordneten Membrankrafte nicht im Einzelnen notwendig, da 
wir hierfir lediglich Aussagen fiir die Spannungsintegrale tiber die Eisenquerschnittflachen be- 
_nétigen. Dabei benutzen wir die aus den Randbedingungen hinsichtlich u langs der Rander AA’ 
und BB’ folgende Beziehung: 


h ow = : 
u(x 2)—4(=—9} = | Bah (18b) 
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die wir mit der Naherungsannahme, ne 


PY Vy + Vo % (18c) 


innerhalb des gesamten Stahlbereiches sei (was sicherlich eme sehr gute Naherung darstellt, 
da langs der nahe beieinanderliegenden Rander AA’ und BB’ diese Bedingung exakt gilt), in 


+ he/2 + he/2 + h,/2 : + he/2 
ou = ou Cie = eS 
Yih, + yohe= yh, = | . +1) 5 | (He ats Bhi | Yexy dy = E, | Tony AY » 
hele hel = hy? —he/2 


also mit (18a) schlieBlich in 
+ he/2 

— ef? 
iiberfiihren kénnen. Damit gewinnen wir einerseits fiir die resultierende erforderliche Schubkraft 
im Eisenquerschnitt df,, = h, dz 

+ he/2 
dT,,, =z J Toxy dy = dent h, =n% dfex» (19) 
—h,/2 

womit sich zeigt, daB auch hinsichtlich der Schubspannungsverteilung ein Stahlquerschnitt nahe- 
rungsweise wie ein n-facher Betonquerschnitt wirkt, was die Richtigkeit der Formeln (6b, c) und 
(7b, c) auch unter dem Gesichtspunkt dieser eingehenderen, die Bernoullische Hypothese im 
Stahlbereich verlassenden Deformationsbetrachtungen unterstreicht, und andererseits unter Be- 
achtung der Gleichgewichtsbedingung 


a6, x or, xy 
eek i : 

os . . oo, Xs . . 

fiir die infolge der Spannungszuwachse es Eisenquerschnitt entstehenden Schnittlastzuwachse 
+ he/2 + he/2 
Cy ie COR ae CED aie + he/2 
Fy (dite x) = dz | = dy = — dz | 5 1 dy =—dz [tory] yo : (20a) 
—he/2 —he/2 
+ he/2 + he/2 + h/2 

é an a a, x — 7 ace es = + h/2 = 

Zldiie,,) =—de | Fay =a FOF = de VF tony] ee — | tony dp. (20b) 
—h,/2 —he/2 —he/2 


Unter Beachtung der Randbedingungen fiir die Eisenschubspannungen sowie mit (18d) folgt daraus 


O nes 
aq (dMex) = — de (%{—%) =0, (20c) 


lope E, 
dq (din, 2) = dz {t,h,—n1,h,} =—dz(n—1) th, =—(n—1) 1, df,. =—(n—1) 5¢y df. (204) 


Der die Eisenverformungen erzwingende Spannungszustand langs der Rander x = 0 bzw. 
x = Ax hat dementsprechend die Eigenschaft, daB er neben der nach (19) berechneten Schubkraft 
in den Eisenquerschnitten resultierend keine Langskraftzuwachse, sondern allein Momentenzu- 
wachse um die z-Achse nach (20 d) hervorruft. Mit 
Voy = Veo Ae, teas 


erhalten wir also fiir den Momentenzuwachs um die z-Achse langs einer Schnittlangeneinheit 


Ayes 
a saa ager 1) = | oe — 24, Xz) df, = —(n—1) Sb (ylea = 2 $5 %xy) Fox, (22a) 
und entsprechend, wenn me Bere in der y-Richtung liegen, 
“BF = (9-1) Bt) 2 aay ta) Fay in) 


Die vorausgesetzte Verzerrungsfigur nach Abb. 5 ist also allein durch Schubspannungen bzw. 
Scherkrafte nicht realisierbar, sondern nur bei zusitzlicher Wirkung von Momenten um die 
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_ 2-Achse, die auf die einzelnen Bewehrungsstahle wirken. Unter diesem Gesichtspunkt lautet die 
Momentengleichgewichtshedingung 


DM gg =O 


fiir das Plattenelement von Abb. 3 unter Hinzunahme der gestrichelt eingetragenen Momenten- 
zuwachse um die z-Achse 


om.,, om, . 
(n,y dy) dx —(n,, dx) dy + ( ac as) dy el a ay) dx == 05 


om, om 
Ney + = (ny. : a ==) he 


also 


Ox 
Mit (10b, c) und (21a, b) erhalt man daraus 


= E 
ves) [Fi:.—(n—1) F, «| ae YD [Fiy—(n—1) F.,] + (n—lE, tay (Soa: lax — Sey I, 4) —=10 


also mit 


Oe F.. = hy, +n F,,—(n—1)F.,=F,, + Fiz =d=F,,—(n—1) F,, 
un 


(n—1) 5,4 Fox =s,d bzw. (u— I)ishy Fey = s,d 
nach (b) schlieBlich 


1 : : 
x3 (o> a?) + Hxy Sa Sy) =0. 


Beachtet man hierin noch (9), so erkennt man, da diese Gleichgewichtsbedingung identisch 
erfillt wird, ohne die Gleichheit der zugeordneten Membranscherkrafte bzw. Torsionsmomente 
fordern zu miissen. Im allgemeinen wird man die zugeordneten Groen jedoch gleich setzen, da die 
Unterschiede insbesondere bei schwacherer Bewehrung vernachlassigbar klein sind. 


4, Die zeitabhangigen Vorginge (Ermittlung der Umlagerungsgréfen). a) Allgemeines. Die 
- Lésung des zeitabhangigen Problems wird grundsatzlich analog den Uberlegungen des vorange- 
- gangenen Abschnittes 3 aufgebaut. Nach Formulierung des maBgebenden Spannungs-Verzerrungs- 
Gesetzes (analog den Hookeschen Gesetzen (5)) kommen wir, wiederum unter Heranziehung der ver- 
_allgemeinerten Bernoullischen Hypothese (4) nach Integration iiber die Gesamtquerschnittflachen 
_ [analog Formel (6), (7), (8), (10)] zu den Zusammenhingen zwischen den (hier allgemein zeitlich 
langsam veranderlich vorausgesetzten) Gesamtschnittlasten und den (nunmehr ebenfalls zeitlich 
veranderlichen) Plattenverzerrungen ¢ und x. Weiterhin werden auch durch entsprechende Inte- 
grationen iber die Teilquerschnittflachen [analog (13)] die Zusammenhange zwischen den Teil- 
schnittlasten und den Plattenverzerrungen ermittelt, so da man dann [analog (14)] durch Kombi- 
nation der bei diesen Schritten gefundenen Beziehungen endgiiltig einen Zusammenhang zwischen 
den Teilschnittlasten und den (zeitlich veranderlichen) Gesamtschnittlasten erhalt. Nach Ab- 
-spaltung der VerteilungsgréBen verbleiben dann allgemeine Ausdriicke fiir die UmlagerungsgréBen. 

Fiir die zeitlich veranderlichen Plattenverzerrungen ¢ bzw. x, auf die sich das Problem letztlich 
zurickfiihren laBt, werden dann wiederum aus den Gleichgewichtsbedingungen am Plattenelement 
[(15), (16), (17)] entsprechende Beziehungen gewonnen, die sich jedoch nun als ein System dreier 
gekoppelter partieller Integro-Differentialgleichungen reprasentieren. Sie lassen sich auf ein System 
dreier simultaner partieller Differentialgleichungen reduzieren. 

Da der Beton neben den elastischen Verformungen zu Belastungsbeginn durch das Kriechen 
und Schwinden auch zeitabhangige Verformungen erleidet, miissen wir anstelle der Spannungs- 
Verzerrungsgesetze (5) im Betonbereich das durch die Kriech- und Schwindverformungen verall- 
gemeinerte Gesetz 


Oy (9) pa 
a(—) =e +0+E | = 


? 
beck eyiee a eet dal FD heel? 22 
= E, Pe) uy Sov) in Poo y Le: E, dp ( Jos : Pro ? ( a 


1 Siehe Funote 2 von Seite 1. 
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heranziehen, wobei ¢,,, bzw. g,, die Endwerte der Schwindverkiirzung bzw. der Kriechzahl be- 
deuten. Zur Auflésung von (22a) nach op,(y) differentiieren wir zunachst nach gy, womit die Dif- 
ferentialgleichung 

ie 


d ; d ah 
+ onslg) =? Z (04x07) = BA leva) + 2g 


entsteht, deren Paras mit der pgaaeee Cc 


Oox(Q) = B, les.) + ae em ? f eyx(Q) e” dy Cee 


lautet. Einsetzen in die urspriingliche ieee (22a) liefert fiir die Integrationskonstante 


Ps 
C=—E,~ + E, [f eox(y) e? dp]p =o 


so daB wir letztlich 


Oop) = E, 


? 
cog 
&x(y) —e ” J Bi AX OL Fry ele =») (22b) 
n= 2 i 
erhalten. Anstelle der Hookeschen Gesetze (5) benutzen wir dementsprechend im Betonbereich die 
SpannungsVerzerrungsrelationen } 


Pp 
€ 
oO, = E, |e, —e” feds a—en], (23 a) 
~=0 za 
? 
é, 
o, =F, epee) A] i es ae) : (23b) 
~4=0 = 
- ? 
c= he Je Vey in. (23¢) 
= 
im Stahlbereich bleiben die Hookeschen Gesetze erhalten: 
eens > 5, hao ae Se (2344) 


Samtliche Verzerrungen lassen sich dabei wieder, da wir auch fiir die zeitabhangigen Vorgange 
die verallgemeinerte Bernoullische Hypothese voraussetzen, durch die Membranverzerrungen in 
den, die ideellen Gesamtquerschnitts-Schwerpunkte enthaltenden horizontalen Plattenebenen bzw. 
durch deren Kriimmungen z. B. nach (12) ausdriicken. 

b) Zusammenhang zwischen den Gesamtschnittlasten und den Plattenverzer- 
rungen. Indem wir nach (1) baw. (3b) durch Integration tiber die Gesamtquerschnittflachen aus. 
den Spannungen wieder die Plattenschnittlasten berechnen, gewinnen wir letztlich nach Einsetzen 
von (23 a—c) fiir den Betonbereich bzw. mit den Hookeschen Gesetzen (23 d—f) im Stahlbereich einen 
Zusammenhang zwischen den Gesamtschnittlasten und den Plattenverzerrungen x bzw. ¢. Manerhalt 


es ale dfg= fo. df, + fie, df= : J ee for dy+- Pepe ae 9|aiern fee ut. 


(Fox) (Fex) (Fox) (F ex) 
me a é, nach (12a) unter Beachtung von (a) und (b) 
¢ é 
n, = E, je%) Fin — Fox | f eX (ex) + sy, x,) dy — "ia (1 — e— ol 3 (24a) 
also entsprechend : 
(sx) iy (Sac) 
ee r Vx 
nee nf ce set [| a + 562% uf. (24b) 
(sy) g (sy) 
Ee Vx eS Vx 
Dye 1, ae x Fy, — Fyy € sl (= + Sy ey) . (24c) 
e é 
ny = E, ey Fiy— Fyy |e—® | e G” Peee Ney = 9 24d 
i ih y Hy) OX ’ (24d) 


1 In (23c) entfallt der Schwind-EinfluB, d. h 
Boe seme = ae ib ua a ein homogener Betonkorper unter alleinigem Einflu8 des Schwin- 
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sowie 


2 =) 04% df, , J Ont df. + J et df. 


9 
mip oe &,—e—? | eX @, ease: “seo oe *) %, Of, 1 | a 
ae) nae 


: (Fe x) 


und wiederum unter Beachtung von (12a), (a) und (b) 


WP 


Ms x =F oes fetta. ya 
0 


Q 
e-? ic (ef) Eros s v_) dy ——= (1 we - ; 


0 


(25a) 
also entsprechend 
& Q P y(sa) 
Many = — By Textey — Suu e—# J oF Yay ts Foee* fo 2 + 5443 dys, (25b) 
: ae) 
Ms yx = — By \Jiy Hey — Joye o0 f et dy — spy Fy e— le (PEE + 565 ey) dye » (25c) 
tr) 
= é, 
WE, == Fifa Tey =o { etx ty dy — Spy Fase o f exit? Yt sy, 2 Hy) dy — 2 (1—e— olf 
6 0 Poo d 
(25d) 


Die Auflésung dieser Gleichungen nach den Plattenverzerrungen erfordert die Lésung jeweils 
zweier gekoppelter Integralgleichungen, da in den jeweiligen Gleichungen [(24a), (25a)], [(24b), 
(25b)], [(24c), (25¢)] und [(24d), (25d) ] stets zwei unbekannte VerzerrungsgréBen ¢ und x enthalten 
sind. Wir betrachten zunachst allein (24a) und (25a), um hieraus die Verzerrungen ¢(*x) und x, 
zu berechnen. Der analoge Aufbau der iibrigen Gleichungen erlaubt dann eine Verallgemeinerung 
der fiir diesen Fall gefundenen Resultate. Wir multiplizieren zunachst (24a) mit s,, und addieren 
zu (25a). Es entsteht schlieBlich 


(Sx) __ Ts 2 Ny, dss x“ Je 


ef 2 ee OS 6 © [ex dy. 26 
- Es, IBS T BF, Spx Fix Six Yee 4=0 ze - 


Setzt man dies in (24a) ein, so folgt letztlich eine Integralgleichung zur Berechnung von x,. Man 
erhalt anit den Abkiirzungen 
Ue Fox WG Jos Fox Te 
a = ge Ee Ox = py 
fe Fix Six : Fix Six 
nach kurzer Zwischenrechnung 


(27a, b) 


0 LO 


P P x 
Hy — (1 — dex + Ox) ev | x, dy -- Ob oof i enn, | dy 
0 


- bh Fs 
lok ne dys (Le?) | ws 
e Je n,, dy = ( e— ?) y 
(27c) 


@ 
1 Fox aie x Sox Fox 
peed ee — e* m, x d ee 
Ey Jix (m. Pies é Bie Pec Jy. 


Ux 0 
2 
Fiihrt man in dieser Gleichung die durch e~ ? = [e?()] zu beschreibende Operation durch, so gewinnt 
man schlieBlich die Sp eee 
e+ (Ll + ee — dip) pt tex Me mR ay aC y) = ae +2 tea. (27d) 


a S 


Ihre allgemeine Lésung setzt sich zusammen aus der Lésung der homogenen Gleichung von 
(27d) 
HH = C, e lx? -+- C, e2x?, (28a) 


-wobei die charakteristischen Exponenten 7}, und 2, Lésungen der quadratischen Gleichung 


z+(1 + be — Xx) Vy + %ex = 9 
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sind, und 


LA ORE seme le + Oe x» —Ob x) 5 /( ae Bex — bx)? —4 dex] ‘ (28b) 


betragen, und dem partikularen Integral der vollstandigen Gleichung (27d), das man zweckmabig 
nach der Methode der Variation der Konstanten ermittelt. 
Insgesamt haben wir demnach als allgemeine Lésung 
ty = Cy @ ix? ae G @ 2x? —— eae for e 2%? dp —o'** [yt @. 1x? dy . (28c) 
24 guarie 
und wenn wir hierin noch y* durch die urspriingliche rechte Seite y der Integralgleichung (27¢) 
ausdriicken, 


Hy = G | Seed Cc, ee? | y eu + Vex)" 22 fg "8P y dg Ot Me) ete? { &—*? do. (28d) 


ges ee 1% Vo" Vix 


Zur Reduktion der bei der Integration entstandenen GréBen C, und C, wird die Lésung (28d) 
in die urspriingliche Integralgleichung (27d) eingesetzt. Man erhalt schlieBlich nach einiger Rech- 
nung mit Beachtung von 

Vix V2x = Nex » Vix + Vee = — (1 + H%ex— Mb x) » 1 sien + Ob, = 2 = Vix + Vox 9 


: Oe = 1 + Vix Vex + Vix + Vex = (1 | Vix) (i + v2) 
die Beziehung 


lc, ead ae (fe oe dg), 2 i “bx Xx (1 — (1 =e Vix) 0) Cin 


TE viz eerie 
ce! ae Vox) EL 1— a, ae he Xb x — Pr 
+|G, + Be Carel (fe ydp), 6} | Taye (eae ¥ ( (1 +n9¢)/¢ ?—0, 
woraus 
sect + Vix oe ee ae dl + Vex)? eae. 
o ae Vo" —Vix [ye dp), Ce "erase Vox Vix (f« v dp), 
und damit nach (28d) fir die Lésung der Integralgleichung (27c) 
p ? 
Me = = (1 rar ee? fe —Vex% dy — ce! ae ee? [ye "dy (28e) 
6 0 


gefunden wird. Driicken wir hierin yp letztlich durch die urspriinglichen Gréfen m,,, n; und das 
Schwindglied aus, so folgt endgiiltig 


4 
1 ee eee Fy, frosts 
agi met het : ‘fa = leh Big — F.. (ms + 520) 6 1st dy 

ple x Be & F * Ses 

= Fra eer z *fia tT ¥2x) Mm Le (ms x vm) ee i 
Fy, Shu F 1+y, ie 
S00 aes ee YIxP Vou vax? | 
Tix Poo On)" ; EE Vo 1x —Vax) a : ) z (29) 


so dafs wir hiermit aus (26) schlieBlich nach langeren Umformungen 


(Sx) 1 diya 
Ex me lie: ; pt B + ae 


a ? 
ont Ve — Spal Six) %2 xP — Von X 
yl eee e Jn. e dy 


= QP 
1 eve Vo xl Ji x) @ ix? i! n ey er dy 
Vyx Vex é vs 


— 


ena Fis She I Sh Die ova f — I ah i oe e xP — 2x 

E, Fin Jix Vix Vox ae ates x Vox fmee es 

ye “sss ee 1 Jon 1x@ 1 Jy 22D) |” 
PES. Gon ve ie we (er) (1 toe Fs) ‘) (39) 
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erhalten. In gleicher Weise ergeben sich die ubrigen Zusammenhan 


und den Schnittlasten. Die Formeln fiir x. und ey) z. B. gehen aus den Formeln (29) und (30) 
hervor, sofern man dort iiberall an Stelle der Indices x den Index y schreibt. Auch die Formeln fiir 
die Schubverzerrungen sind analog aufgebaut. Beispielsweise gehen die Formeln fiir Yr 


aus (29) und (30) hervor, wenn man dort %, durch x,y, ef*) durch iS y@x), n, durchn,, sowie 
: : 2 oxy v 
m,, durch m,xy ersetzt und den Schwindanteil Null setzt. 


ge zwischen den Verzerrungen 


(sx) 
y und yee 


c) Zusammenhang zwischen den Teilschnittlasten und den Plattenverzerrungen. 


Durch Integration der Spannungen bzw. der Spannungsmomente iiber die Teilquerschnitte er- 
halten wir die Teilschnittlasten. Es ergibt sich 


me = { on df, = E, [ fee feed + Sea—eol at 


(Fb x) (Fo x) ~4=0 
2 & 
=E, Fis | + sxex,—ee { (ell? + sary) ot dy Hegiey) Nei 
iO es 
und enstprechend 
(sx) 2 (sx) 
Nay = Eb, La = + Shx%xy — & ? I a + 54%.) e*% dy ’ (31b) 
7—0 : 
(sy) P inplSy) 
ze Vx 
Nyx = E, Foy soe + Shy %xy —e_” fl aT + Sy Hes] oa,| ; (31c) 
%=9 
n é 
my = EF, Foy ev) Shy Hy — @—* f (e°9) + Spy %y) e% dy aera, —e—?)|, (31d) 
%Z=0 
9 
ae [o.m:df.=—E, ive net | omar) (32a) 
(Fb x) KO 
und entsprechend 
2 
y=0 
P 
Mp y x SS faye [rey e* in) (32c¢) 
y=0 
? 
my = — FE, Joy (,— ie e” in) (32d) 
x¥=0 
sowie 
Nex = foo, df, Ss E, fe — *, tt) df, ae E Le (el? aay Se 2) (33 a) 
(Fe x) (Fe x) 
und entsprechend be 
Vee 
roay = Bi Fea a — sates) (33b) 
' (sy) 5 
mega =B, Fay ("aE —siyHey) (30) 
Tig ae tee ee (ef?) — Sey Hy) » (33d) 
sowie : 
Mex = foes df, =—E, Jex Xx» (34a) 
(Fe x) 


und entsprechend . 


Mex y ae og eys Mey x =— E, Jey May Mey =—E, Jey xy- (34 b—d) 


4 3 
ey 
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Mit (24), (25), (31), (32), (33) und (34) erhalten wir noch folgende Indentitaten: 


_ Foz Mz — Moe — Mex ce Eb) a Mie ete ee 
Web 3 aig nn ee Se Nex = F n, Ee ’ 
ix Sha T Sex ix Sox T Sex 
Ny n et Are ECE, n n Fox Ney 4 Wied Soe eee a eey, F 
ease oa Specs i iad ASE a 0 eel She + S, 
x bx ex ix bx ex 
in ees : Meyvn yn 1 Mey x oa CEG os Moye 7 Mb yx Meyx 
= pe J es = 5 
by x Fi, "y % Sby + Sey ? Jere yi& F,y coma Sby + Sey 
rs _ ey Mm, |, — Miya Ty * ts te Mi, ig —— WY, 
oar 9 ey let | 
y Fiy y Sby T Sey F;,y y Sby 1 Sey 


d) Zusammenhang zwischen den Teilschnittlasten und den Gesamtschnittlasten. 
Mit den durch die Gesamtschnittlasten nach (29) und (30) (sowie den hier aus Platzgriinden nicht 
niedergeschriebenen analogen Gleichungen) ausgedriickten Membran- und Biegeverzerrungen € 
und x der Platte, kénnen wir nun, indem wir die Plattenverzerrungen ia (31) bis (34) einsetzen, den 
Zusammenhang zwischen den Teilschnittlasten und den Gesamtschnittlasten angeben. Wir be- 
trachten auch hier aus Platzgriinden allein die Normalkrafte und Biegemomente in der x-Richtung, 
da fiir die itbrigen Gro8en wieder véllig analoge Formeln gelten. Einsetzen von (29) und (30) in 
(31a), (32a), (33a) und (34a) liefert schlieBlich nach einiger Rechnung 


a a P 
hie Joa Vix VIX F, —11 x7 
= 4 if ——— x 1xX% 
Mp x Tie t aC ae e ( Se Vix) m Fr (ms x + Six na) e dy 
6 

¢ 


eS a ate ee, | la + V2") Ms ~— = (ms + Spx no) a dad | 


Vix — Vox 
0 
é, She % y 
es F, oo bx *bx e axP ix? 35 
Po Vix — Vex ( )s ( ) 
= =. P 
Rides Decl eae rae F, Ever 
eka x = 1 = * 1%% 
aoe Sex a elie Oss ( + M1) mee F;, (rms a a) : u% 
0 
P 
L+ Vex vex” [ F, —» 
ee ee ee 1 ogy Ed as Qu 
re 7, e | (1 + vax) Ms x F;. (m,, +- Spx n.) e dy 
0 
i Soe Wes 1) ery ; 
ea F,; $00 ex “ex 1x J] — 1%” ou C8 (|e Y2%P 
See CES ae igor itd (36) 


F Sha Ef 
bx __ “be ~ bx 


Sage vi Six ie 
= a Q 
(L+ 1x) Sex + Vix Sox soe Fee —1%% 
SiG eee ayers (1 + 12) Oe, ee =} Sy¢ 2,))e= *** dy 
ey S e 
htt + Vex) Sex + Vax Jox = Le ? —VYoxX 
Tix Sex + Se) Cin —Pox) 6 : ee Be nen (mss ae oe Me) oe dy 
ye FoalP (9,2) wd, Ig ue (1-92) 2 Jee + Vax Ji 
E,F yx f We x. = 1x. ex lx JY¥bx ee 2x ex Vox bx 72x? 
+ eke Pp a Vin“ — Vox Jak M a GS. e 2 (37) 
nF gy, Tt Fox Sex 
Ne x ee F;, n, + Ji Ms x 


P 


prow (1 +91.) 0 Jon + Vix Jee Vix P Fy» —V1Ix 
Jix (Sex Fie Spx) (Vie —Vex) : . ES 1) ee UEP (ms « xs oe ng) ‘ rig dy 


0 


Q 
' (l +22) Jon + Vox Fe Ura Z Fy see Ok 
+ ee eae | (CL ted mee — BE (me + soa me) 7 dy 


rs Fog! Bie { (L294) Beg Puede ee (L+ 2.) Jee +Vox Jp 
Sy A Ty | i 1x ex 1x Ybx« 10 ae & e Vox % 72% 
B, be Pao f s Poel Fie ( JixVix a Tre ; e ‘| , 38) 
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Dabei stellen die unterstrichenen jeweils ersten Ausdriicke (vgl. Abschn. 3, Formel (14)) die Ver- 


teilungsgréBen dar, wahrend der Rest dementsprechend der allgemeine Ausdruck fiir die jeweilige 
UmlagerungsgréBe ist. / 


e) Die Verschiebungsgleichungen des Problems. Mit (24), (25) baw. mit (31) bis (34) 
sind sdmtliche Plattenschnittlasten durch die Plattenverzerrungen ausgedriickt worden. Daher 
bendtigen wir zur endgiiltigen Lésung der Aufgabe Beziehungen fiir die Plattenverzerrungen bzw. 
fiir die Verschiebungen der Mittelflache. Analog Abschn. 3 gewinnen wir sie aus den Gleichgewichts- 
bedingungen (15a—e) am Plattenelement}. Beachten wir, daB 


M0x = Mz x = Sells Moxy = Msxy + Sy Mey » 
Moy = My + Sy Ny » Myx = Msyx + Sy Nyx 
ist, so erhalten wir zunachst aus (16c) 
om @2 02m. 2 
Sx sy ra) e? 
Ber a Sa By By Mey + m, x2) aL ay? oe axe (s, n,) a Ox ay (s, Nx y = Sy Nyx) 
2 
+ yz (Sy my) + Pp, = 0. (39a) 


Hierin werden die Beziehungen (24) und (25) eingesetzt, wobei wir die Membranverzerrungen der 
ideellen Schwerebenen nach MafSgabe von (9) in die Mittelflachen-Membranverzerrungen umrech- 
nen. Man erhalt dann nach einiger Rechnung die Integro-Differentialgleichung 

oe a? ? 
Px (Jiox Hy) Gsoy [(Jiox =e Jioy) | shyt (Jioy %y) 
ete 2 { & @ @ 
=e” | e & (Jt0x %z) + Ox dy [(Jtox + Jeoy) %ey] + dy? (Ivo )| dy 
: @2 Ge ROSE Se ee @2 
— Ps — F, . S, €0 x) 2 x 4.) + aya (Fiy Sy €0y) 


E, 2 
oe Dy x Pe ey so by a 
| 5 eee + ya (avy Boy eoy) dy 


Ox Oy 


= oe? 
1 6 | e% = (aox Fux 0x) + 
} 
oon 2 62 
a (es) Ze (tox Fx) + pye (My F,,)| ° (40) 
Dabei bedeuten J,ox bzw. Jyoy die Tragheitsmomente der Betonquerschnittflachen F,,, baw. F, E 
hinsichtlich der, durch die Plattenmittelflache verlaufenden Achsen und a,, bzw. ap y die Schwer- 
punktsabstande der Betonflachen von der Mittelflache. Fiir die Membranverzerrungen der Mittel- 
flache erhalten wir aus den beiden ersten Kraftgleichgewichtsbedingungen (15a, b) die maBgeben- 
den Beziehungen. Einsetzen von (24) in (15a, b) liefert, nachdem man wieder die Membranver- 
zerrungen der ideellen Schwerebenen nach MaBgabe von (9) in die Mittelflachenverzerrungen um- 
gerechnet hat, nach einiger Rechnung 
@ 


: Dif ee cg wes: a CS eae 
x=0 


Ox Oy 


Pp 


(7) 0 (7) (7) 
Ee, Fig) — Ee Piy 85 tes) ——? | ot 1E (ase Four) + & (ony Foy res)} dy 


4=0 
Ce oF, . ap 
7) aoe aid Ox ieee (41) 
: Q 
x é Y0x (2) 
2(r.2822) + 2 (riyeoy) —eme { o* [2 (ma) + 2 eyes} dy 
Ly ae : 
Enea ee (on i, e-? fee Pere a Ps x} dy 
Bee ae xe fix xy ay ap Uy aap . ax x y dy y xy. 
3 ue 
ay Se ny AP a 42 
prguate a6 ?) By 10037 are (42) 


1 Die zeitliche Veranderlichkeit der Lasten sei hinreichend klein, um das Problem als ein statisches ansehen 
zu k6nnen. 


4 
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Zur Lésung dieser Gleichungen wird man stets auf die zugehérigen Differentialgleichungen 
zuriickgehen, die man dadurch gewinnt, daB man mit e? multipliziert und hernach nach g differen- \ 
ane den Spezialfall einer zur Plattenmittelflache symmetrischen Bewehrung (a x = yy = S, | 
= s, = 0) zerfallt das Gleichungssystem, wie auch in der gewohnlichen Plattentheorie, in eine |j 
Gleichung fiir die Mittelflachendurchsenkungen w und in zwei gekoppelte Verschiebungsglei- | 
chungen fiir die Scheibenverschiebungen uy und v9 

e2 

Pa 


2 e2 
(Jon %e) + seg Wow + Joy) tay] + pe Tioy Hy) 


dy= (43 a) |, 


p 
A 2 e2 
AY rs | e% fe (Joo 2%). ++ iy [(Joox + Jooy) %xy] + gye Joo y %) 
6 
QP 


a ) F. ee la 0 F,, 
# (Fis000) +2 ($2 7045) 8 r |e (xen) +2 (42 yong) te 
0 


Belg gt) See OS. ah) 


a (Fix a a a (Fox a ; 
5 (Fer .05) + ay Fixe) |. tiglet roma iy Phy eas) 


oes ea oF, y Py 
ate is ee (43 c) 


und speziell bei konstanter Bewehrung in 


btw o4w o*w 


Jiox 5a + (Fiona lies) Ox? Oy? t Jioy aa 


P 
o* o4 e 
ema ie Ive, et ne Seay. ae dy =e (44a) 


é 


sowie in 


Q 
2uy Ey 8 [du , av ee 07Ug yy 7) e } + Pe 4 
Bis Ga 13 ae 9 ine ha ay \ay + @x/| ° ie aa Sa 

ft) 


Fa 


v 


2 0x 


2 Lox Ox 


? 
du dv Guy is Fi, 8 [dug ed 6v9| : Pipe: 
gto ia: ai (Bt a) + Fy Bi tact = 0. (te) 


oy 


5. Beispiele: Die konstant und zur Mittelflache 
symmetrisch bewehrte, allseitig gelenkig gelagerte 
Platte unter Quer- und Axialbelastung. a) Allge- 
meine Lésung beikonstanter Randbelastung. 


aK Von der, unter einer lotrechten Belastung 
PA*,¥) = > D4G,, sinja = sin ka (45) — 
j=1lk=1! se b 
{ und den konstanten Randdruckkraften D, und D, 


stehenden Rechteckplatte der Abb.6 wollen wir || 
voraussetzen, da sie nur im Koordinatenursprung || 
unverschieblich gelagert sei, so daB die durch das_ | 
Kriechen und Schwinden bedingten Verformungen || 
der Platte in ihrer Ebene lings der Schneiden- 
lagerungen x = 0, x = a, y = 0 und y = b nicht be- 
hindert werden und somit auch keine in der Plattenebene wirkenden zusitzlichen Auflagerlasten | 
langs der Schneidenlager mobilisiert werden. Mit | 


tos ¥) = ug) = Ag) +25, vex 9) = 09) = Bl) -y (46a) | 
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ist dann ei die geometrischen Randbedingungen sowie die Verschiebungsgleichungen! (44h, c) 
hbefriedigender Membranverschiebungszustand gefunden, so daB wir aus (24) mit? 


: (sx) aes ou, s 0 
fe" = Ce = GP = Alp), ey” = ry = FY = BY), 
(sx) IR sy) OUy ov ae 
Ysy Vee Yory dy a bx : 
schlieBlich 
pe F, 
fies, A(p) = Fi, 9 |e A(y) dy + ae T= -] 
0 co 
; Fe (46b) 
n, = E, | Fiy Bip) — Foy oo [fo B(x) dy + Te (1 — 1) ; 
0 co 
(OA See 
erhalten. ; : 
Die dynamischen Randbedingungen langs der Schneidenlager 
n,(0,y) =1,(a, y) = — D,, n(x, 0) = n,(x, b) =— D, 
kénnen mit (46b) nur erfillt werden, sofern 
n, =—D,, Ny aes (46c) 


ist, womit fiir die zundchst unbestimmt gebliebenen Funktionen A(y) und B(g) aus (46b) zwei 
Integralgleichungen hervorgehen. 


Wir betrachten allein die erste Gleichung weiter und erhalten aus n, = — D, die Integral- 
gleichung 
’ 
F, x &s D 
Fi, A(p) — Fox 9 [ et A(y) Wa (Ue es area e (47a) 
Po E, 


0 


und nach Anwendung der Operation e—” % [e?( )] unter Beachtung von F;, — Fy, =n F. x 


dA n Ee Fi, x é, D 
hid Alo WMS eral 47b 
Se aeR yo seme hia Yan = Fy Fig ee 
mit der Lésung 
nFex 
ee = F;x v Bye et D, 
ewe e. Saks te wk, ee) 


Einsetzen der Lésung in die urspriingliche Integralgleichung (47a) liefert schlieBlich nach 
Koeffizientenvergleich hinsichtlich e” 
y BY fae 
n Fy Poo BE, 0 Fie Fin” 


so daB endgiiltig 


nF. x nFex 
Eye os ( $5. Fig :) ‘De le ( oh Fyx 4 
MO) a F Pao l—e Sek bis are l—e 9 (48a) 
und entsprechend® 
Zs OT aaibooey ple Fi PLOT EL RS 2 Ay ess F; 4 
B(g) = Rune Poo 1 e€ ys; E, Fiy I + nF ey il e x ( 8b) 
hervorgeht. 


Um den Einflu8 der Membrankrifte auf die Biegung der Platte beriicksichtigen zu kénnen, 


: fiigen wir, ebenso wie bei der Untersuchung von Plattenbeulungsproblemen zu der gegebenen 


- Belastung p, noch die lotrechten Komponenten der (nunmehr am verformten System angreifenden) 


1 Hierin sind p, = p, = 0 zu setzen. aN 
2 Bei symmetrischer Bewehrung sind s,, = Spy = Sx = Sy = %x = [by = 0, und die ideellen Gesamt- 


: querschnitts-Schwerpunkte liegen in der Mittelflache. 


3 Dieses sind dieselben Funktionen, die auch beim Balken das durch Bewehrung behinderte Schwinden 


und die Verkiirzung durch eine konstante Druckkraft beschreiben, vel. z. B. FuBnote 2 von Seite 1. 


2 
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Membrankrifte hinzu. Dies fiihrt, wie bekannt, im allgemeinen Falle dazu, das man in der Platten- |) 
gleichung anstelle der Querbelastung p, mit der ideellen Querbelastung 
0?w 0?w 


= dw 
Ps = Ps + 2 a + (ngy + Nyx) ax dy + Ny ay? 


rechnen mu. In unserem Falle ist also 

Sater D 0*w o*w 

Pe Pee Ox2 DS oy? 
in (44a) einzufiihren, so daB wir endgiiltig fiir die Biegeflache die Integro-Differentialgleichung * 
aw dw Ce me 0-w Dy ew 


r ° f | | 
Jix axt t (Jix + Jiy) ax? Ay” t Jiy dy | BE, ox? ' E, ey? 
“2 4 
ee = 64a = = o4w S atw| De of och dj heme MeN Be Ye y 
ay {fe ree, + (Jou + Joy) id by? Spy yt { OX = 2, ON) ea (49a) | 


0 
zu lésen haben. 
Die Randbedingungen des Problems umfassen die geometrischen Aussagen 


w(0,y) =0, Ww Gs ¥) — OF, wx, 0) = 02 w(x, b) = 0 (49b) | 
sowie die dynamischen Bedingungen 
m0.) =O. m, (a, y)=0, m,,(x, 0) =0, Msy(x,b) = 0. (49¢) | 


Mit Riicksicht auf (25), woraus mit s,, = spy = 0 eine alleinige Abhangigkeit der Momente von den 
zugehérigen Kriimmungen folgt, lassen sich die Bedingungen (49c) auf 


6?w ?w 
*A0>¥) = 5a | =0, AO Y) = gaa| = 0% (49a) | 
? 2 d | 
Saks ay? |e 0 ° a dy? |.6 


zuruckfihren. 
Zur Lésung von (49a) benutzen wir den, die Randbedingungen (49b) und (49d) befriedigenden 
Ansatz 


w= > SOXy) snja— snka”. (50) |} 
j=Hlk=l a b 


Hinsetzen in (49a) liefert nach Koeffizientenvergleich hinsichtlich sin j a — sin kx ~ das System 
der Integralgleichungen 5 . 
() ? 
Qi #8 H*{ ly) ed tie 
0 


Gj k 
FA] 


(51a) 
wobei zur Abkiirzung 


;: peat [7222 Ey Jix E, 72 k2 24 k2 a2 [IP 2? EB, J; 
= fe = Del ee lates pee eae D,), (51) 


a a a? b? b2 b2 


(b) 54S 72 Jg2 — — k4* — 
Af) = Bot |B Ste ae ox by) + Edy “(Sey 


gesetat ist. Mit der weiteren Abkiirzung 
©) 4) __ j2a?/jeatn By Jew PP ke at 7 7 k2a® (ha? n Ey Je 
Aik =) —M ik = a — D, ape 0 Ey (Sex + Jey) + ie ( a ae o) 
(51d) 


wobei n Jen = Jix oe Jos baw. nJey = Jey Ji y die n-fachen Tragheitsmomente der (symmetri- 
schen) Bewehrung hinsichtlich der durch die Mittelflachen verlaufenden Achsen bedeuten, erhalten 


: eres : cd] 
wir aus (51a), nachdem wir die Operation aa [e?()] angewendet haben, die Differentialgleichung 


5 yD? 
Hadi ah 


1 Da wegen der Symmetrie der Bewehrung samtliche Schwerpunkte in der Mittelflache liegen, ist J; = J 
Jioy = Jas Sion = Jeu Jhoy = Soy: Nig 


ix? |i 
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19 
deren Lésung 
a 
Were Sine? 
D; ke C. ‘7k ate he 
Ae) 
jk 
lautet. Nach Einsetzen in die urspriingliche Gleichung (51a) wird nun die Konstante CG, zu 
() 
Ge eee Ki k 
jk 4) roy 
jk “jk 
bestimmt, so daB wir endgiiltig als Lésung 
(e) 
eee 
ey Gj k | Aik \ iy 
jk Ai) 10) ee (Sle) 
s J 
haben. 
Nach (50) ist somit die Lésung fiir die Biegeflache gefunden 
2 ( ah 
2h Sie aie a) a de y 
w= po) She ey eee ee sin je sin ke, 52 
1 1 ay i A) ety: ee ( ) 
und wir erhalten die Plattenschnittlasten im Sinne von (25) (mit.s,, = sy, = 0) 
ie ” 
m,, —=— E, { Ue — Jb x re? ee “| 
0 
(e 
(6) eee) Sil 
me E, ee se ie Gk Tir i k —= Jy x Mh at) : . fo es y 
= ee Ile. 5* sin j 7 -— sink 7 —, 
a a 2, A) Bae ave) J IRE Ts Aes: 
und entsprechend 
if ag 
- qe) —_ () (i) 
cE, sik Jig4 jk Jt | A eae y ) 
eee AE hy Proms Ee: J cos ja— coska—, 
be ab a oe a? be hd ay a b 
AO) 
Be en () ae (i) Ses 
7 Ey IR GRY Jit = Joy 4h ( 4 he 
Mm, 4 = — ; Jiy + l—e Jj cos ju—coska7, 
yx ab 1 2 a) oy, a) J b 
(¢) 
J 
2 (6) a) wen 
ae "dik Jae Jy 4 | ay Rape ee ¥ 
Ms y 5" = a z Jiy ae) ewes Sn he Ia 
J 


_Dabei stellt das unterstrichene Glied jeweils den Anteil dar, der unter Zugrundelegung rein elasti- 
scher Verformungen entstehen wiirde. Der jeweils zweite, von y abhangige Ausdruck gibt den Anteil, 


der durch die zeitabhangigen Verformungen des Betons bedingt ist. 


Wir betrachten noch zwei Spezialfalle, bei denen wir zwei grundsatzliche Erkenntnisse ge- 


winnen. 
b) Die Platte ohne Axialbelastung. In diesem Falle ergibt sich aus (51h, c, d) 
(i j2. 772 k?2 7? A ey sa I pe ke? 7? Ey Jiy 
AY =(%3 a9 b2 ( a2 ae b2 2 
=i): 2m?  heae\ (2, hay Jie 
ai a ( 7D on b2 a a b2 ’ 


70 = L ae eae hs s: be Pes. 
je ? 


a b? a b? 
yi * 


(54a) 
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also = ee see 
ut Ee dese eee (54b) | 
qa ORE RT + Ma J, 1 
yk Ux ly 
und insbesondere 
dns on UU psa Slee tees Om De | 
A —T1 —— 1g Vee Fey) Jy —J,, M12 Vey — Fea) | 
— = So ee ae : = »  (S4e) 
at RL eta a) ay On gen) 
tS 


woraus wir erkennen, daB® fiir Jer = Jey (Platte isotroper Biegesteifigkeit) sich die | 
Plattenschnittlasten zeitlich nicht andern, also die aus der titblichen Theorie ela- jj 
stischer Platten sich ergebenden Werte beibehalten werden. Sind jedoch die Hisen- | 
Traigheitsmomente in den beiden Richtungen verschieden, so tritt eine Umlagerung — 


der Schnittlasten im Laufe der Zeit ein, und zwar in dem Sinne, da die Richtung |f 


mit der gréBeren Eisensteifigkeit mit der Zeit groBere Momentenbeanspruchungen | 
erhalt, wahrend die Querschnitte mit der kleineren Hisensteifigkeit gleichzeitig eine | 


Entlastung erfahren. 
Um die GréBe dieser Lastumlagerungen, die aus der ungleichen Be- |} 


ae wehrung herrithren, abschatzen zu kénnen, betrachten wir eine qua- |} 

; | SNE , dratische Platte mit a =b =I, die (was den extremalen Fall darstellt) jj 
; : ee Pa in y-Richtung unbewehrt ist. Fiir eine Belastung I | 

2 7 | 

t , aS 


; Eade y 
1 Pz = sing — sina 7 


erhalten wir dann mit 

0 = 6e(n—1)p,, 
wobei u, = F,,/d der Bewehrungsprozentsatz und ed der Abstand der (symmetrischen) Hisen- |} 
einlagen von der Mittelflache sind (Abb. 7), schlieBlich mit J = d?/12 | 


qitsina — sina > nop 
w = : ey pee, ee 
t@E,Ja+e | * no 


P1+2 1 ——_—*? —|\| . : 
<a ae FTE Fa) cine sina, 


1+2¢@ 
ae a “ soe F Ea = 0 ( ee Ho=itra)| e082 e rend 7 
Ms yx = My yx = — ne i + -6 WFD Cos 7 — COS 7 = F 
my = my = 2 ae entre sina | sing 7, (55) | 
und weiterhin aus (32) und (34) | 
ae Get te on nay 
Nin S se wae = 1 4 ee ( a ~w=itra)| sin 7 > sin 7 + ; 
a ieayite’ oo mage 
Mexy = =F aoe 5 + a ( —oe=tra)} cos 7% = Cos 7% *. 


Im allgemeinen betragt die Plattenbewehrung maximal 1%. Setzen wir ferner Emax © 0,4 und |i} 


fe 10, so wird : 


Omax = 6-0,16-9- 0,01 & O71: 


Ninax 
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Mit diesen Zahlenwerten erhalten wir dann aus (55) die in Abb. 8 dargestellten Kurven, die eine 
thebliche Belastung der bewehrten Querschnitte zugunsten einer Entlastung der unbewehrten 
y-Richtung zeigen. In Abb. 8c, d sind die jeweils auf den Beton- bzw. Stahlteil entfallenden Mo- 
mente in ihrer zeitlichen Abhangigkeit dargestellt. Die auSerdem strichpunktiert eingetragenen 


mes(P49), Mes[th) rel ho) _reyalvbg) macltbq)  many(t0e) real PEg),_Mexylsig) 
Tsa(z} 0) Mszey(4l,0) My, LE) Meyzo(t,1,0) Moa, Lyb,0) Meezy|,1,0) Men|¥9) ™Mesy| 11,0) 


War i: Ler aaee: Cc 
Abb. 8, 


Kurven sind diejenigen, die sich ergeben wirden, wenn mit den anfanglich vorhandenen, aus der 
ublichen Theorie elastischer Platten berechenbaren Momenten 


l l Pe er Gee eo 
mse» Oy 0) it, Als l, 0) Ag? ieee 0 
die Teilschnittlasten berechnet worden waren ohne Riicksicht darauf, da die Gesamtschnitt- 
lasten m, , bzw. m, x y selbst zeitlich veranderlich sind. Es ergeben sich fiir diesen Fall aus (35) und 
(36) nach einiger Rechnung 
2nep 


me il SES psd! 68 l= 2 0. Gee ee) 
ine “Ge p) = —Fves(l l, ~) See (a1) (l 6)” 9 


2nep 
pls a sees a os aie ( ins  @=iittaa)| 
Thee(t > +39) =— Fire BS les aya ee 2no S : 
Dabeierkennen wir, da die tatsachlich verbleibenden Momente im Beton- bzw. Stahlteil erheblich 
erofer sein kénnen als die unter der Annahme konstanter Gesamtschnittlasten berechneten Werte. 


c) Die mit p, = q: sin ™* sin sowie durch konstante Randdruckkrafte D, und D, 


belastete Platte. Der Dischinger-Effekt. Aus der allgemeinen Lésung (52) und (53) verbleibt 
allein das erste Reihenglied 7 = k = 1 und wir erhalten mit 


Dee 
bat (Pe + 32) 
sowie 
eA (al (ef) 28h (Nhe 
poo = | + (5) | f+ 525) | = ee HCG) | ae) | | 
a Ey Jy x a diy a 3 ee mt Ey Jy y ey Jo x a 56 
B= Ff 3 (e) |i ie ie =e aa ee a i fe a hal (56) 
p25 [Gif eel] belli 
schlieBlich aus (52) und (53) 
ie: Ba) too ( — oy fea "| sin 7 * sina oe 
WTO l=7 Rays ee 
an ; = a __ 4e0—8 
, See oe oh 4 Choa Soaks) + Ka8(y_, =F sina Zainab, 
Ee ogi 8) [2 ma a 
— = Aco — 9 
me: : =F ? 
q =H | (Ji:4o0 Juehs) + Ie, 9 —6 ee cosa, (57) 
= 6ab4;,— 9) — 4,9 —9 
‘ = = 4e0— 9 
on aS qz LE, Te ee Mehl ths", AGO "heart coent. 
>. bh, —o (1 Lio a 
E = = Aca 8 
a qu Ey | 7 4 Uokom Soh + Ke8( ios "| sin ~ sin >. 
a) = 0, —9) [2 Ay —8 a J 


\\ 
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Wir erkennen hieran, da8, zundchst einmal den Ausdruck A, 9 —6 betrachtend, der zeitabhangige | | 
Teil bei wachsenden Axiallasten, also mit wachsenden Werten 6 gréBer wird. Wahrend fiir A, 9 —0 > 0} 
das zeitliche Anwachsen schwicher als linear mit y sich vollzieht, haben wir fiir 7.9 —d < 0 
(Aio —6 > 0 vorausgesetzt) ein mit gy starker als lineares Anwachsen. Der Grenzfall A, 9 — 6 = 0) 
fiihrt, wenn man beachtet, daf nach der L’Hospitalschen Regel i\ 


Ago -=2 P eee 
SS P 10 
lim Le: tae = lim ieee ‘ = td 
(Aeo — 6) > 0 Ag — 9 (eo —6) +0 1 Aip — 0 


ist, zu einer linearen Zunahme der Verformungen und der Gesamtschnittlasten mit g. Hier er- 

geben sich also dieselben Erscheinungen, die man auch bei der genauen Untersuchung des Langs- | 
kraft-Biegungs-Problems fiir Stabe! beobachten kann. Wahrend der Ausdruck /,9 —0 nur bei| 
den zeitabhingigen Verformungen von Einflu8 ist, wirkt sich der jeweils im Nenner stehende 
Ausdruck /;9 —6 auch auf die elastischen Verformungsanteile aus. Bei kleinerwerdendem A;9 — 0 | 
ist das Anwachsen der zeitabhingigen Verformungen und damit der zugehérigen Schnittlasten | 


weitaus starker als das der elastischen Anteile. | 
Der Grenzfall 
Ao —6 =0 (58) 


bedeutet ein Versagen der Platte, und dies entspricht gerade ihrer Beulbedingung. Auf diesen 
starken Einflu8 einer Axialbelastung, die in der GroéBenordnung der Eulerlast liegt, auf die zeit- |} 
abhangigen Verformungen und Schnittlastzuwachse hat erstmals Dischinger? hingewiesen, der |} 
diesen Effekt an einem Stabe ohne Bewehrung untersuchte. Die zweidimensionale Verallgemeine- 
rung dieser Untersuchungen, also die unbewehrte Betonplatte, ergibt sich aus (57), wenn wir 
ex = Jey =9, Siz = Jiy = Sez = Joy = J = #12, also 


Zio =Av0 ase, = i. = () | <4 Ae ==") (59a) i 


b 
setzen. Mit der Abkiirzung 


a\ il 
a n# By a+ (2) 
a 


mee ess 2 (59b) 
“oem 
folgt aus (57) schlieBlich 
qa* 1 v — : ; 
cat Hay | a tea eae sina — sina, 
1+(2 ? 
b 
2 ? 
Meg =e (=) er sina — sina +, 
1+ (— 
(60) 
qa 1 y == 
Ms %y = Ms y ¢ = nb | ieee "cos — cosa 2, 
1+ 3) 
b 
q b 1 y = 


Zum Vergleich werden die Formeln von Dischinger fiir einen beidseitig gelenkig gelagerten Stab iI 


angegeben. Fiir ei ee aoe ee 3 ; ve : : | 
oe Pa eine Querbelastung q(x) = q) sin z ; hat man bei gleichzeitiger Wirkung einer | 


1 Hieriiber erscheint in Kiirze eine Arbeit vom Verfasser im Osterr. Ing.-Arch. | 
® F, Dischinger, Bauingenieur, 18 (1937), S. 487, 538 u. 595. i} 


yey \! 
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die Durchbiegungen und das Biegemoment 


j4 oe stage ; P eae 
w(x) = aE aT f +(e =A sina =, M(x) =; — eat sin 


6. Zusammenfassung. Das Problem des Plattenverbundes la{t sich, nachdem man bei der Dar- 
stellung der maSgebenden Beziehungen weitgehend in Analogie zur Theorie elastischer orthotroper 
Platten vorgeht, letztlich auf ein System dreier simultaner Integro-Differentialgleichungen fiir die 
Mittelflachenverschiebungen zurickfihren, die man fiir konstante Steifigkeiten mit Hilfe von 
Produktansatzen lésen kann. Im allgemeinen Falle der verainderlichen Bewehrung wird eine exakte 
Lésung kaum noch zu erreichen sein. Man kann in solchen Fallen nach der Methode Sattler mit 
Hilfe angenommener Zeitabhangigkeiten die Gleichungen zunichst einmal auf reine, nur Ab- 
leitungen nach x und y enthaltende Differentialgleichungen zuriickfiihren, fiir die dann z. B. im 
Sinne des Ritzschen Verfahrens eine Naherungslésung aufzufinden ware. 

AbschlieBend sei noch bemerkt, da mit der gelungenen Darstellung des Zusammenhanges 
zwischen den (i. a. zeitabhangigen) Schnittlasten und den Plattenverzerrungen nunmehr auch im 
Rahmen der Loveschen Naherung der Schalenverbund allgemein untersucht werden kann. 
Durch geeignete Einfiihrung der Vorspannstahlquerschnitte in die ideellen Gesamtquerschnitts- 
werte kann man auch das Problem der vorgespannten Flachentragwerke erfassen. 


(Eingegangen am |. Februar 1960.) 


Anschrift des Verfassers: Privat-Dozent Dr.-Ing. R. Trostel, Berlin-Charlottenburg 2, HardenbergsstraBe 34, 
Lehrstuhl fiir Mechanik. 
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Uber gesteuerte Anheizvorginge bei Zylindern* 


Von R. Ansorge 


1. Einleitung. Eine Reihe von Arbeiten der letzten Jahre befaBt sich mit der theoretischen Be- 
handlung gesteuerter Anheiz- und Abkihlungsvorginge. Das sind Vorginge, bei denen die Grenz- 
schicht des einen Kérper K umgebenden Mediums nach einem gewissen Programm aufgeheizt oder 
abgekiihlt wird. Gefragt wird nach der dadurch bedingten, mit der Zeit veranderlichen Temperatur- 
verteilung in K. Die Integration der diese Vorginge beschreibenden Anfangs-Randwertprobleme 
darf vom Standpunkt der reinen Mathematik unter der Annahme temperaturunabhangiger physi- 
kalischer GréBen (Warmeleitfahigkeit, spezifische Warme usw.) als gelést gelten’. Die Lésungen 
erscheinen in Form unendlicher Reihen, die numerisch oft nur mit groBem Arbeitsaufwand aus- 
werthbar sind. Das Problem wird naturgemaf einfacher, wenn man die Zahl der auftretenden Orts- 
koordinaten herabdriicken kann. ; 

In den technischen Anwendungen tritt haufig dic Frage nach der Temperaturverteilung in 
Vollzylindern bei gesteuerten Anheiz- (oder Abkiithlungs-)-vorgangen auf. Diese Aufgaben wurden 
theoretisch bisher meist unter Voraussetzungen behandelt, die es nach Einfiihrung von Zylinder- 
koordinaten r, py, z erméglichten, das Problem unabhangig von y und insbesondere von 2 zu be- 
handeln. Solche Voraussetzungen sind: vélliger Temperaturausgleich innerhalb des Zylinders vor 
Beginn des Anheizvorganges, nur zeitabhangiges Aufheizprogramm des AuSenraumes, unendliche 
Lange des Zylinders. Zur Begriindung der letztgenannten Voraussetzung wird die Erfahrungstat- 
sache herangezogen, daf fiir den Fall groBer Zylinderlange gegen den Zylinderdurchmesser die L6- 
sung des ebenen Problems (Abhangigkeit der Lésung nur von r und der Zeit t) die gemessene 
Temperaturverteilung gut wiedergibt. Das kann natiirlich nur fiir Zylinderpunkte gelten, die nicht 
in unmittelbarer Nahe der Deckflachen liegen. 

Es erhebt sich die Frage nach der GréBe der EinfluBzonen der Deckflachen und bei gegebenem Zy- 
linderdurchmesser nach der minimalen Zylinderlange, fiir die der Zylinder bei gegebener zulassiger 
Abweichung noch als unendlich lang betrachtet werden kann. In der vorliegenden Arbeit wird 
versucht, darauf fiir einfache, in der Technik haufig vorkommende Spezialfalle eine befriedigende 
Antwort zu geben. 

Da iiberdies der Kinflu8 der meist nur schwer vorauszubestimmenden Warmeiibergangszahl a 
interessiert, werden die numerischen Untersuchungen fiir verschiedene GréBenordnungen des 
dimensionslosen Parameters A/a ro (siehe Ziff. 2) durchgefiihrt. Dabei mége a den gesamten Warme- 
iibergang durch Leitung, Konvektion und Strahlung beriicksichtigen. 


2. Verwendete Bezeichnungen. Es soll bedeuten 


] 
A eae cad| die Warmeleitfahigkeit, 


0 ca | die Dichte, 

1 
c ; - al die spezifische Warme, 

i 

O eemeGad die Warmeiibergangszahl, 
T [Grad] die Temperatur im Zylinder, 
T, [Grad] die Zylinderanfangstemperatur, 
Tz [Grad] die Temperatur im Zylinder endlicher Lange, 
T., [Grad] die Temperatur im unendlich langen Zylinder, 


* Auszug aus der gleichnamigen, von der Fakultat fiir Natur- und Geisteswissenschaften der Bergakademie 
Clausthal seneieee Dissertation des Verfassers; Referent: Prof. Dr. phil. H. Konig, Korreferent: Prof. Dr. 
rer. nat. M1. Menzel. 


1 H. Jackel, Ing.-Arch. 26 (1958), 146. 
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@® [Grad] die Temperatur der Grenzschicht des den Kérper umgebenden Mediums 
(AuBenraumtemperatur), 

A den Laplace-Operator, 

Jos Ji die Zylinderfunktion erster Art (Besselfunktionen), 

To [cm] den Zylinderhalbmesser, 

2% [cm] die Zylinderlange, 

= = | eine Abkiirzung und demgema 

LZ. einen Zylinder mit 2) =x ry, 

T.,, [Grad] die Temperatur im Zylinder Z,,. 


3. Die Lésung des Anfangs-Randwertproblems fiir den Vollzylinder. Ruft eine Warmezufuhr 
oder ein Warmeentzug in einem beziiglich der Warmeleitung homogenen und isotropen Kérper 
innerhalb des betrachteten Temperaturbereiches keine Anderung des Aggregatzustandes hervor, und 
entsteht oder vergeht im Kérper selbst keine Warme, so werden die gesteuerten Anheizvorgange 
durch 

oT A A (eT 7 

ee ee (Do Er, Parents fir t>0 (3.1) 
beschrieben. Dabeisoll F darauf hindeuten, daB die Werte langs der Kérperoberflache zu nehmen 
sind, wahrend n die nach aufen weisende Flachennormale darstellen soll. 

Bei der nun folgenden Betrachtung des Vollzylinders soll zwar auch von einigen der eingangs ge- 
nannten vereinfachenden Annahmen Gebrauch gemacht werden, namlich 


1. « hat fiir alle Punkte der Zylinderoberflache den gleichen Wert, 


2. T4 = konst, d. h. vélliger Temperaturausgleich innerhalb des Zylinders vor Beginn des Auf- 
heizens, 


3. Das Aufheizprogramm @ ist nur zeitabhangig, also 0 = P(t) 
(wir haben es dann mit einem zylindersymmetrischen Problem zu tun), jedoch behalten wir die Ab- 
hangigkeit von z bei. Dann lautet (3.1) 


2 2 
erate ae Bes Tg ast 
Ot oc 


ér2 | or @r Oz? 


yal 
(3.2) 
ae A. T (ror % t) = D(t) | fiir t>0. : 


Fiir den Fall T, -+ ©(0) kann von der Lésung des Problems (3.2) beit = 0 nur verlangt werden, 
daB sie im Innern des Zylinders mit T,, iibereinstimmt, wahrend sie fiir t > 0 auch den Temperatur- 
verlauf auf der Oberflache richtig beschreiben soll. 

Es bedeutet nur eine Verschiebung der Temperaturskala, wenn wir T, =O0setzen. Bestimmen 


_wir T aus (3.2) mit T4 = 0, so wird der Fall T, == 0 dadurch erledigt, daB man T durch T— Ty 
und @(t) durch O(t) — Ty ersetzt. 


Mit der von Kneschke1 eingefiihrten, von einem Parameter t abhangenden Einflu8funktion 


G(r, 2, T, t) versuchen wir den Ansatz 


T=t 


t 
Tro, = C(r-7:0, 0) | lee G(r. z,%)t—7)| dr —E(r, % #, 0). (3.3) 
0 


Dieser erfiillt bereits die Forderung T(r, z,0) = 0. Er erfiillt auch die weiteren Gleichungen der 


_ Aufgabe (3.3), wenn man 


G(r, 2, T, t) = — P(t) g(r, t) h(z, t) (3.4) 
A, Kneschke, Ing.-Arch. 24 (1956), 77. 
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setzt und g und h folgenden Bedingungen unterwirft: 


(3-3) 
und 


(3.6) | 


Das System (3.5) kann durch den Ansatz 


2 
ou 


ee) Oy ay 
g(r, t) = 2 Ay, u,(r) € reas 
gelést werden, wenn man die u,(r) aus der Eigenwertaufgabe 


” i 5 C) eaiee 
Un tun + (7 u,=0, | | 
(3.7) 


y ee 
a Un(To) + u,,(To) =0, 
u,(0) = 0 
Doreunett und dann die Koeffizienten A,(n =1,2,...) nach dem Entwicklungssatz! so wahlt, 


da Dy A,,u,(r) im Innern des Zylinders mit der Entwicklung der Funktion f(r) = 1 nach den 


is Gifanktiones der Aufgabe (3.7) iibereinstimmt. Es ist 
=e (0, = 
To 


Dabei stellen die «, die positiven Lésungen der Eigenwertgleichung 


A 
J(@) =—— wW J,(@) (3.8) 
aT 
dar. 
Fur A, erhalt man den Ausdruck 
A 

eae 

ar, 


A 


(3.9) 


n 


1+ 2 x2 zal Fane 


Da wir spater neben x auch den Paramter J/« ry variieren und iiberdies zu dimensionslosen 


Koordinaten iibergehen wollen, schreiben wir statt g(r, t) in Zukunft g ae f und er- 
halten ae age ra 


2 
“7p % To \.OCT? 


A r 
ene 2 a 
8) ar, r, SO Sel 
s(z, sik at \= Sy : lari ia 0883 (3.10) 


1 Courant-Hilbert, Methoden der mathematischen Physik Bd. I, Berlin 1924. 
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Auf gleiche Weise lésen wir (3.6). Es ergibt sich 


2 z 
h («5 3 Aes i} S ee cos ( A a! 
{GT o: lpeeOetes js Pa cos (Qn %) sin (Q,,2%) + Qn “ph (3-11) 


Die EKigenwerte Q,, sind dabei die positiven Lésungen der Gleichung 


A 
Mit (3.3), (3.4), (3.10) und (3.11) lautet somit die (eindeutige) Lésung von (3.2) (die wir sogleich 
fiir beliebiges T4 = konst formulieren) abkiirzend 


Mr, 2.2) =O) + [T.— FO I6 (3 2, nea) h(x fee ) 


2 Aaa? \ pare 2 
ao To QCTo a1 To QCcro 


t 


ED) fhe tA Tigh aaa 
| dt (ci, Acca (t ») (is ot, aoe (¢ )) ae. (3.13) 
0 


Dies gilt fiir T, +: O(0) beit = 0 nur fiir das Innere des Zylinders. 
Die Lésung des ebenen Problems (Unabhangigkeit von z beim beidseitig unendlich langen Zylin- 
der) erhalt man aus (3.13), wenn man dort h= 1 setzt. 


4. Die Lésung des ebenen Problems als Grenzfall der Lésung fiir den Zylinder endlicher Lange. 
Die in der Literatur angegebenen Beispiele zur Ermittlung der Temperaturverteilung in zylindri- 
schen Kérpern (Vollzylinder, Rohre usw.) bei gesteuerten Anheizvorgaingen behandeln fast stets 
das ebene Problem. Die Verfasser gehen dabei (z. T. unter Berufung auf die Erfahrung) von der 
Voraussetzung aus, daf fiir groBe x die Lésung des ebenen Problems die gesuchte Lésung fiir den 
endlich langen Zylinder hinreichend genau wiedergibt (wenigstens fiir die von den Zylinderenden 


- hinreichend weit entfernten Zylinderpunkte, z.B. den Zylindermittelpunkt). Genauer formuliert 


heiBt dies: Es wird unterstellt, daB die Lésung im Falle des endlich langen Zylinders fiir x —> oo 
gegen die Lésung des ebenen Problems konvergiert, und zwar gleichmaBig innerhalb des betrach- 
teten Temperaturbereichs, d. h. vor allem gleichmafig fiir alle t aus 0 < t < t) mit beliebig grobem, 
jedoch als gegeben zu betrachtendem ty. 

Wir wollen diese Annahme anhand der Lésung (3.13) wenigstens durch einen entsprechenden 
Beweis fiir die Punkte der Zylindermittelebene (z = 0) stiitzen. Wegen der Beschranktheit von 


e(z pee 4 ‘) und 


2 
CRity = teen (One 


t 
d®(t) Ane? A ex 
| dx ® ares ST: t ))ae 


0 


A 
fir 0 <t < ft geniigt es, zu zeigen, dab h G 7 3205 


cer ) mit x —> oo gleichmabig fiir alle nicht 
0 
negativen t < t, gegen 1 konvergiert. Das ist in der Tat der Fall. 


Zum Beweise vermerken wir zunachst, daB nach (3.12) die Werte Q,, x (m =1, 2,...) aufgefabt 


werden kénnen als die Schnittpunkte der Geraden y = : fore mit y = ctg x. Wegen ae >0 
folgt daraus 
Gn Onn e, (4.1) 

also 
; cos (Q,, x) sin (Q,,%) >0. (4.2) 
Weiterhin erkennt man sofort, dah 
| sin (2, %)| > [sin (Qm41%)]- (4.3) 
ist. Aus (4.1) und (4.2) folgt 
: cos (Q,, #) sin (Qn %) + Qm%< C08 (Qm +1 %) sin (Qm 41%) + (ORI (4.4) 


Ingenieur-Archiv 
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A 


re Fs : 
Beachtet man Q,, << Qm 41,80 resultiert aus (4.1) bis (4.4), daB h(x ae 0, ae t) fiirt = Oeine 


alternierende Reihe darstellt, bei der die Betrage der Glieder eine eigentlich monoton fallende Folge 


bilden. Daraus folgt 
sa el 


MstP . — Q2 ——_t! 
Ay 1 = asin (5) ee 
bs C08 (Qqm%) Sin (Qn 2%) + Qn% | 
m=M+1 
. — 22 SSE 
2 |sin (2), . 1 %)} : M+ perk 2 a 
= Ta 4 SE, S = at Le ae. Ge a > 
= cos (Qqyy , 1%) sin (Qy 4 1%) + Qy 41% Qy 4.1% “Ma 


somit sicher A}; .1< ¢ fir M = a = M,(e), und M, hangt nicht ab von x und t. 


a ) konvergiert folglich gleichmaBig fiir alle x und t > 0. Wir 
LG 

haben daher die Méglichkeit, den Grenziibergang x —> 0o in dieser Entwicklung im wesentlichen 
gliedweise durchfiihren zu kénnen. Hierzu werden noch die Grenzwerte 


lim 2, Frome Oe a lim 2% 


Die Entwicklung h(x; ae Sil) 
Rire 


u%—> © x—> CO 4—> CO 
fiir beschrankte m < M und 0 <t < t bendtigt. Aus (4.1) folgt 
lin: 2 = lim —2;, = 0 (4.5) 
x—> 0 x—> CO 


fir beschrankte m und ¢. Gleichung (3.19) ergibt 
ctg (Q,, x) =4.9,= te (- — Lf On) 
0 


also 
2m—1 A 
as a — Q,,% = arc te(s,, 2) ‘ 
wobei wegen 
2m—l1 


RK 3 


Tt — Dy <a 


rechter Hand der Hauptwert zu nehmen ist. Mit (4.5) folgt daher 
2m—1 
2 


(4.6) 


lm 2% = 
x—> 


fiir alle beschrankten m. 
Wegen (4.5) und (4.6) existiert somit zu gegebenem positiven ¢ bei gegebenem M und f fiir alle 
m< Mundt < | ein 
% = #,(M, €, ty) 
derart, daB fiir alle x > x, die Beziehung 
M 2 sin (Q,,% Gtx 4G gaa se 
hes Onan hom Gee oe to ats ee crest 


{yeaa 


besteht. Auf Grund der Leibnizschen Reihe ist 


Poet pote yn |Silxe 


4 
cs 


firm > M,(e), also fiir alle positiven t < t, 


co 


~ 2 sin (Q,,%) —22 Ree 1 
cos (2,2) sin (Qn) + Qn is nee 


m=1 


<aeorers 


wenn Man nur 


x > my (M(e), E, ty) = x(€, ty) 
wahlt mit 


M(e) = Max {M,(e), M,(e)}. 
Damit ist unsere Behauptung bewiesen. 
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Unklar bleibt natiirlich weiterhin, was bei der Approximation der gesuchten Temperaturvertei- 
lung in einem, hinreichend langen“ Zylinder durch die in einem beidseitig unendlich langen Zylinder 
unter ,,hinreichend lang“ zu verstehen ist. 

Diese Frage wollen wir in den Fallen konstanter AuBenraumtemperatur* und zeitlich linearer 
AuSenraumaufheizung numerisch untersuchen. Wir stellen deshalb zunichst die zugehérigen Lé- 
sungen als Spezialfalle von (3.13) zusammen. 


5. Die Lésungen in den Fallen konstanter Aufenraumtemperatur und zeitlich linearer Aufen- 
raumaufheizung. Bei konstanter AuBenraumtemperatur 


P(t) = 0 
liefert (3.13) fiir den Zylinder endlicher Lange 
A r A (ers A 
Te = 9 + (Ta— 9) 8(55 ee ii (us = gen’): (5.1) 


Entsprechend erhalten wir fiir den beidseitig unendlich langen Zylinder 


To =O+4(T, na(Z jas i). (5.2) 


? 2 
a1 To Q@CcTo 


Bei zeitlich linearer AuBenraumaufheizung 


P(t) =Ot+@, 
ergibt (3.13) 


’ 2 
aM %) ?OcCr 


t 
A r A A Zz A 
eee BORO ALS. a ates Arora ‘ Dro 
@ | 8(en ih gea")*(Sinin goa") * (5.3) 
0 


Te = Ot +O) + (Ta—)a( Fs : alist owe ) 


bzw. 


t 
A r A Ayregst A 
Tap = Ot + Oy + (Ta— 8 (5 Es ceat\—9 | (seo gente GM 
t 
Der Ausdruck (5.4) kann fir grofe t geschlossen dargestellt werden; denn fihrt man die Integra- 
SS 
tion aus und vernachlassigt die Glieder mit e *eerg , so erhalt man 


2 “+ 
T,.(r, t) = Ot + BD, — % = zo DS “i : o,\2 fir groBe t. (5.4) 
Di (rs on 2) 1000 (2) 


Mit (3.9) steht rechts die Reihe 
r 
CO A, Jo (1 a 
Se a: 
"2 
%(r) geniigt im Innern des Zylinders wegen der dort bestehenden Beziehung 


= r 
DA h(on7)=1 


n=1 


der gewéhnlichen Differentialgleichung 


FB") +H) +r =0 ; 
mit den Randbedingungen 


AF (04) + Blro) = 0, 0) =O. 


et, sondern als ,,natiirlicher Aufheiz- 


“2 ee nn) . Sofi nt heizvorgang gewert 
* Dieser Fall wird haufig nicht als gesteuerter An gang g sing Steasruns, wealialb 


vorgang“ bezeichnet. Doch liegt auch im Konstanthalten der Umgebungstemperatur 


wir diesen Fall terminologisch nicht ausklammern. 
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Die Lésung ist 


; __ [To A i 
(7) = ae ic 7) 4 ? 
so da wir wegen der Eindeutigkeit der Lésung bei zeitlich linearer AuBenraumaufheizung fiir groBet ||P 
den geschlossenen (und von T unabhangigen) Ausdruck 


Ooc|[ro [To A Te Ab) 
je roe ee ae 


T..(7,t) =Ot+G, 


erhalten. 

Nach einiger Zeit ist also die Differenz zwischen AuSenraumtemperatur und Zylinder- 
temperatur nur noch von r abhingig, fiir ein bestimmtes r somit konstant. Diese Erscheinung tritt | 
natiirlich auch beim Zylinder endlicher Lange auf, doch ist hier die nur von r und z abhangige Tem- 
peraturdifferenz nicht geschlossen darstellbar. 

Die Lésungen (5.2), (5.4) und (5.4b) fiir den beidseitig unendlich langen Zylinder finden sich im 
wesentlichen schon bei anderen Autoren: (5.2) steht fiir } = 0 z. B. bei Gréber-Erk1, (5.4) und || 
(5.4b) wurden fiir T., = ®, von Kneschke angegeben. 


6. Der zeitliche Verlauf der Temperatur im Zylindermittelpunkt. bzw auf der Zylinderachse. 
a) Die dennumerischen Untersuchungen zugrunde gelegten Werte der Parameter | 


Alarg und x. Um nun den Einflu8 von x auf den zeitlichen Temperaturverlauf im Innern des | 


Zylinders studieren zu kénnen, wollen wir die Mittelpunktstemperatur (r = z = 0) fiir verschiedene 


x beikonstanter AuBenraumtemperatur und zeitlich linearer AuBenraumaufheizung auftragen. Wir | 


wahlen x = 1, 2,4, 8,00. Von der Veranderlichen t gehen wir zu der dimensionslosen Verander- 
lichen 
A 


2 
(Oi r 


— 


tiber. Eine Anderung von A, @, ¢, fy bedeutet dann in t-Richtung lediglich eine Anderung des MaB- | 
stabes. 

Von wesentlichem Einfluf8 auf die Temperaturverteilung ist weiterhin die meist nur schwer 
vorauszubestimmende Warmeiibergangszahl a. Neben x variieren wir deshalb den (dimensionslosen) 
Parameter A/x ry), dem wir die Werte 0,1; 1; 10 beilegen. In diesen Bereich diirften neben vielen 
anderen Fallen z.B. bei etwa 0,5 em < rpg < 4.cm und den (bei Luft als Umgebungsmedium) vor- 
kommenden Warmeiibergangszahlen der GréBenordnung nach Erden, Baustoffe (ohne Hélzer) und 
dgl. fallen 2. 

Da vornehmlich die Temperaturerhéhungen T,; — T 4, bzw. T,, — T 4, interessieren, werden 
wir statt der Mittelpunktstemperaturen die ebenfalls dimensionslosen GréBen 


Tz (0, 0,t)— T4 Ay ¥ _ AL = 
B—T, = 1—8\ar, 3+ h oy Oey 
Mpipta—e(2ia) | 

0 


o—T, 


bei konstanter AuBenraumtemperatur und 


a b 
T,(0, 0,1) — T = Tt T 
£(, 0, t) A =| [ts (gn307)h (54,507) er, baw. 
0 0 
0 


bzw. 


(6.1) | 


@ ets 
_ (6.2) 


“< t 
To(0,t)— Ta _ Byes 
pe et 


Cae F 
bei zeitlich linearer AuBenraumaufheizung auftragen. 
Dabei haben wir im zweiten Fall wie Kneschke ®, = T'4 gesetzt. Wir gehen also von der haufig 


praktisch erfiillten Voraussetzung aus, da8 vor Beginn der AuBSenraumaufheizung schon volliger 


-Temperaturausgleich des Zylinders mit seiner Umgebung stattgefunden hat. (6.2) erhalt man dann 
durch Integration von (6.1). 


1 Gréber-Erk-Grigull, Grundgesetze der Warmeiibertracun > Berlin/Gotti Heidelb 
* A. Schack, Der industrielle Warmeiibergang, Diisseldorf 1329. a ee 
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Unter Tz, bzw. T,,, (die Argumente lassen wir im folgenden fort) verstehen wir in Zukunft 
den Temperaturverlauf im Zylindermittelpunkt baw. auf der Zylinderachse. 


b) Die Eigenwerte fiir verschiedene Kombinationen der Parameter Alary und x. 
Um die Ausdriicke (6.1) und (6.2) numerisch auswerten zu kénnen, miissen zundchst die in den 
Funktionen g und h auftretenden Eigenwerte w, und Q,, fiir die verschiedenen Parameterkombina- 


3 A 3 
tionen (« = ermittelt werden. 
0 


Da die Entwicklung von h bei groBem % sehr langsam konvergiert (insbesondere bei kleinem t) — 
fiir z = 0 im Grenzfall x = 00 wie die Leibnizsche Reihe —, steht zu erwarten, da zur Erzielung 
einer gewissen Genauigkeit mit wachsendem x immer mehr Reihenglieder und damit Eigenwerte 


. ‘ enue A lee: 2 
heranzuziehen sind. Andererseits ist h («3 poe 0,¢) eine monoton fallende Funktion von i, die sich 
0 


bei kleinem ¢ und groBem x nur wenig von 1 unterscheidet, so daB z.B. bei x — 8 im Rahmen un- 
serer Rechengenauigkeit fiir t < 1,8 keine Funktionswerte h berechnet zu werden brauchen. Nun 
nimmt mit wachsendem t die Zahl der benétigten Reihenglieder schnell ab. Folglich wachst die 
Anzahl der zu berechnenden Eigenwerte 2, mit wachsendem x nicht so stark, wie zunachst erwartet 
werden mu8. Da die Kigenwerte w, und Q,, auch bei anderen Aufheizprogrammen auftreten und 
damit von allgemeinerem Interesse sind, stellen wir die in dieser Arbeit benétigten Werte in Ta- 
belle 1 und Tabelle 2 zusammen 


Tabelle 1. Eigenwerte w,, 


n | 
4 | 1 | 2 | 3 4 
= | | | ma 
0,1 / 2,17950 5,03324 1,95688 10,93633 
1 | 1,25578 4,07948 7,15580 10,27100 
10 | 0.44168 3,85771 7,02983 10,18329 


Tabelle 2. Eigenwerte Q,, 


1 1,42887 | 4,30580 | 7,22810 | 10,20025 

2 0,74806 | 2,24574 | 3,74770 | 5,25583 | 6,77099 | 8,29320 0,1 
4 0,38374 | 1,14948 | 1,91616 | 2,68335 | 3,45121 | 4,21986 | 4,98939 | 5,75984 

8 0,19392 | 0,58178 | 0,96967 | 1,35758 | 1,74554 | 2,13357 | 2,52167 | 2,90985 

1 0,86033 | 3,42562 | 6,43729 | 9.52932 | — 

2 0.53844 | 1,82180 | 3,28916 | 4,81478 | 636115 | 7,91680 1 
4 0,31615 | 0,98379 | 1,70350 | 2.45297 | 3,21694 | 3,98841 | 4,76411 | 5,54241 

8 0,17473 | 0,52829 | 0,89079 | 1,26186 | 1,63927 | 2,02093 | 2,40544 | 2,79189 

1 0.31105 | 3,17309 | 6.29906 | 9,43527 

2 0,21642 | 1,60197 | 3,15742 | 4,72297 | 6,29113 | 7,86034 10 
4A 0,14825 | 0,81589 | 1,58653 | 2,36675 | 3,14953 | 3,93334 | 4,71769 | 5,50233 

8 0,09888 | 0,42180 | 0,80093 | 1,18859 | 1,57870 | 1,96984 | 2,36148 | 2,75343 


ce) Der Verlauf der Temperaturzunahme im Zylindermittelpunkt bei konstanter 
AuBenraumtemperatur. Aus (6.1) folgen fiir die Parameterpaare (%, A/a 19) verschiedene 
Funktionen von i, die fiir kleinere # zunichst stark ansteigen und sich dann asymptotisch dem Wert 1 

‘nahern. Aus diesem Grunde berechnen wir die nun zu ermittelnden Funktionswerte anfangs iiber 
einem engeren, ab ¢ = 2,0 iiber einem weiteren Raster. Wir wahlen folgende Abszissen t: 


Ge Onde 040.6 -.5, 1,8062,0- (3,0. 4,02 13,0: 


Diesem Bereich 0 < t < 13 entspricht z. B. beieiner Temperaturleitfahigkeit //o ¢ = 0,003 cm? sec}, 


wie sie etwa der GroBenordnung nach bei Erden (z.B. gepreBtem Al,0,-Pulver) vorliegt, und 
ry = 0,5 em der Bereich 0 min St < 18 min. 


Ingenieur-Archiv 
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Die mit Hilfe der in Abschnitt b) angegebenen Eigenwerte an den genannten Stellen berechneten 
Es A Nes) = : 3 33 Fi 5 
Werte der Funktionen g ie 3 0, | und h Ho 0,z) finden sich, soweit sie fiir die weiteren 
xT T 


Untersuchungen benétigt wurden, im Anhang. Ausihnen ergeben sich nun nach (6.1) die Tabellen 3, 
4.und 5. Graphische Darstellungen der hier tabellierten Funktionen finden sich im Anhang. 


Pi 7y Dt ae ee ee 
Tabelle oe es Py bzw. Fy as Fs ° fiir To 0,1 
Ph { 
aa 1 2 A 8 oo 
t 
0,0 0 0 0 0 0 
0,1 0,1284 0,1000 0,1000 0,1000 0,1000 
0,2 5022 4010 3998 3998 3998 
0,4 8693 TI4T 7656 76056 7656 
0,6 9664 9195 9094 9093 9093 
0,8 9914 9718 9650 96049 9649 
1,0 9978 9902 9866 9864 9864 
12 9994 9966 9948 9948 9948 
1,4 9999 9988 9980 9980 9980 
1,6 1,0000 9996 9993 9992 9932 
1,8 1,0000 9999 9997 9997 9997 
2,0 1,0000 9999 9999 9999 
3,0 1,0000 1,0000 1,0000 1,0000 
4,0 1,0000 1,0000 1,0000 
Teen: 1 Pavey ge ae? ae 
Tabelle 4. o—T,” bzw a fiir at 1 
0,0 0 0 | 0 : 0 0 
0,1 0,0299  —-0,0232 0,0232 0,0232 — 0,0232 
0,2 1728 1301 1298 1298 1298 
0,4 4666 3658 3580 3580 | 3580 
0,6 6637 5510 5314 5314 | 5314 
0,8 71884 6862 6584 6581 | 6581 
1,0 8668 | 7822 AD ES 71506 7506 
2 9163 8494 8193 8181 8181 
1,4 94.73 8961 8689 8673 86073 
1,6 9669 9284 9051 9032 9032 
1,8 9791 9507 9314 9294 9294 
2,0 9870 9660 9505 9485 9485 
E 3,0 9987 | 9947 9905 9894 9894 
4,0 9999 9992 9982 9978 9978 
5,0 1,0000 9999 9997 9996 | 9995 
6,0 1,0000 1,0000 9999 9999 = 9999 
7,0 5 1,0000 1,0000 1,0000 1,0000 
8.0 ; 1,0000 1,0000 1,0000 


d) Die Abweichung der Temperaturerhéhung im Mittelpunkt des endlichen 


Zylinders von der des unendlichlangen Zylinders bei konstanter AuSenraumtempe- 
ratur. Wir sind nun in der Lage, fiir den Zylindermittelpunkt die relative Abweichung der Uber- 
temperatur T; — T', des endlich langen Zylinders Z, von der des beidseitig unendlich langen Zylin- 


ders Zo, also von T,, — Ty, in Prozenten von T,, — T, fiir jedes der benutzten x anzugeben. 
Wir bilden demnach fiir ¢ + 0 die Ausdriicke 


(Ui 1 ad Be) 
aoe 100 


und erhalten und Benutzung von (6.1): 


A = A = 

Ne a= 3 h ze BE 

(Tp —T4)— (Ts — Ta) 99 — elena] 
Leo deg a 


—1]100. (6.3) ff 


bs 
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oy aie o— 


Tabelle 5. oes bzw. Poo 14 = ithe Jae 10 
Dy a 


To 


0.0 0 0 0 0 0 
0,1 0,0034 0,0026 0,0026 0,0026 0,0026 
0,2 0222 0164 0163 0163 0163 
0,4 0740 0539 0524 0524 0524 
0,6 1262 0938 0887 0886 0886 
0,8 1757 1338 1237 1235 1235 
1,0 2224 | 1730 1573 1570 1570 
1,2 2665 2111 1899 1892 1892 
1,4 3081 2479 2217 2203 2203 
1,6 3473 2831 2523 2501 2501 
1,8 3843 3168 2820 2788 2788 
2,0 4192 3489 3108 3064 3064 
3,0 5662 4887 4403 4294 4293 
4.0 6760 5986 D474 5308 5305 
5,0 7580 6849 6348 6145 6137 
6,0 8192 7525 7056 6835 6821 
7,0 8650 8057 7629 7404 7384 
8,0 8992 8474 8091 1874 7848 
9,0 9247 8803 | 8464 8260 8230 
10,0 9437 9060 8763 8576 8543 
11,0 9580 9262 9004 8836 8801 
12,0 9686 9420 9199 9049 9014 
13,0 9766 9545 9355 9223 9189 


Da die Werte von g und h mit wachsendem ¢ kleiner werden, steigt die Anzahl der angebbaren 


Stellen der Quotienten ae mit wachsendem ft. Wir erhalten die Abweichungen (iiber 0,1%) wie 


sie in Tabelle 6 dargestellt sind. 


Tabelle 6. (Pp Tale Ga FA) 100 [%] 


Pagel’, 


0,1 28,4 — teas Fisk Geek eT slice an eee bes | oe 
0,2 25,6 0,30 meme Poe ae coef | 0,2 es hh oer eee cGy eh 
0,4 13,5 1,19 PEGG SO 23 I ERS CT a a ee eee 
0,6 6,27 sid EG OS OO) S369 ch noi — AO Ae Ora 0,1 = 
0,8 2,74 Cte ste 198 £OG NOES a ve ay Deeg as 0.16" | 
1,0 ibs 08th en aS 4,21 | 0,1 SE ATT ale 10229170, 19. hice 
1 0,470 U,186 bis Shs or 89,0 B03 F045 se) =e 40 Oude sO -0,37 =| — 
1,4 0,189 = ee 9,22 | 3,32 | 0.18 |- —_4 39,9 |-12,5 | 0,64 | — 
1,6 = 2s eS ee 7,05. |-2,79 .|.0,213)) .— 1936,9. | 13,2 | 0,83 | — 
1,8 = = hea a 536-270" | 0999 | a7 oe iso LAS | 
2,0 ae = see Nee 4,06 | 1,85 | 0,217; — | 36,8 | 13,9 | 1,44 | — 
3,0 o Se ee 0,946 | 0,545/ 0,115} — | 31,9 | 13,8 | 2,56 | — 
4,0 es ae == ies 0,208] 0,139; _— | — 1 27,4 | 12,8.) 3,19 | — 
5.0 = ae eee is Ss eee ong il.6+ | 3,44.) 0,13 
6.0 = Fe iC at Sioa) ame VO 10;8 99): 3:45°"|. 0,21 
70 a re Sie Ae ae Siete imal Teo aie 9,1 B21) 3532-10527 
8.0 av ces oS = ae ei a6. 17,98) 310-1, 0,33 
9.0 ia aS Sri ee ia Sake eee 2.41 6,96) 2,84 *| 0,36 
10.0 ihe = = ieee ee ee eel 10,5-7| 6,05) 235881 0,39 
11.0 a oa = ee = ee eae eet ore.85 |. 5.23.1 2.31 | 0,40 
12.0 cae = ae ane = ee 7-46 | 4,51:1"2,05 |} 0,39 
13.0 aid yee Poa set a ey ee 0:28. 3,88 | 1,81: | 0,38 


eo 
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Die Werte zeigen auBer dem zu erwartenden Ergebnis, daB bei wachsendem x und festem a || 
(A und r, denken wir uns gegeben) die Abweichungsmaxima fallen, folgende Erscheinungen: . 


1. Die Betrage der Abweichungsmaxima steigen mit fallender Warmetbergangszahl beifestem x. |} 


2. Die Stelle des Abweichungsmaximums verschiebt sich mit wachsendem x bei fester Warme- } 
iibergangszahl in Richtung wachsender #. 


3. Die Stelle des Abweichungsmaximums verschiebt sich mit fallender Warmeiibergangszahl bei 
festem x in Richtung wachsender t. i! 


Je gréBer also die Warmeiibergangszahl ist, desto kleiner darf bei gegebener zulassiger Abwei- | 
chung das kleinste x sein, bei dem der hier behandelte Anwarm- (oder Abkiihl-)-Vorgang noch als || 
ebenes Problem durchgerechnet werden kann. Fiir den Zylindermittelpunkt und den von uns be- || 
trachteten Bereich des Parameters A/a ry kénnen diese kleinsten Werte x = xj, aus Tabelle 6 | 
abgeschatzt werden. La®t man z.B. eine Abweichung von héchstens 1% zu, so diirfen bei kleinem | 
« etwa die Falle x < 6 nicht mehr mit dem Fall x =o identifiziert werden. Doch darf man wohl }| 
sagen, daB die unteren Grenzen fiir x tieferliegen als man gemeinhin annimmt. 

Nun interessiert natiirlich bei vielen Aufgaben nicht nur die Mittelpunktstemperatur, sondern || 
die Temperaturverteilung im ganzen Zylinder. Hier ist gewiB die Annahme berechtigt, daB fiir |i} 
Zylinder mit einem gréReren x als dem hinsichtlich der Mittelpunktstemperatur bestimmten min | 
nennenswerte Abweichungen von der Lésung des ebenen Problems in solchen und nur in solchen 
Zylinderpunkten zu befiirchten sind, deren Abstand von einer der Deckflachen kleiner ist als rp xin. 


e) Der Verlauf der Temperaturzunahme im Zylindermittelpunkt und ihre Ab- 
weichung von der des unendlich langen Zylinders bei zeitlich linearer Aufen- 
raumaufheizung. Die Werte von (6.2) fiir verschiedene t bei verschiedenen Kombinationen 
(x, A/a ro) kénnen nun aus den Werten der Tabelle 3 bis 5 leicht durch numerische Integration be- | 
rechnet werden. Mit diesen Werten kénnen anschlieBend fiir t > 0 die uns wesentlich interessieren- || 
den relativen Abweichungen 


Sa ee f (L—g-h) de 
CS Ramee) Saal Ue £06 == che z : 14100 [%] (6.4) 


Co ad t 
fd-sd 


ermittelt werden. Wendet man jedoch auf (6.4) den erweiterten Mittelwertsatz der Differential- ||} 
rechnung an, so erhalt man 


1— g/—— ;0, € )-h (x; —3;0,€ | 
ae er Gap a mers, 
ro — Te) 199 = ae —1]100[%]  (6.4a) j 
(imi 2 


mit 0< = &(t) <i und erkennt daraus sofort, das die Maximalwerte von (6.4) nicht gréfer sind |} 
als die entsprechenden maximalen Abweichungen in dem zuvor behandelten Fall konstanter AuBen- |} 
raumtemperatur (diese Feststellung kann natiirlich ebenso fiir jeden anderen Punkt des Zylinders |} 
bewiesen} werden). Wir wollen deshalb darauf verzichten, fiir alle bisher benutzten Werte von |} 
A/x rq die Temperaturkurven hier wiederzugeben. Wir rechnen lediglich noch den Fall A/a ry = 1 
fiir r = z = 0 mittels bekannter Quadraturformeln durch. An diesem Beispiel kénnen wir den Ein- 
flu8 von x untersuchen und werden sehen, daf bei gegebenem x der Maximalwert von (6.4) die |} 
GréBenordnung des Maximalwertes von (6.3) auch wirklich erreichen kann. 


Bei diesen Betrachtungen bestimmen wir 


0,2 
f (l—gh)d 
unter Benutzung der frither berechneten Werte von 1 — g h bei 0; 0,1; 0,2 mittels der Kepler-Regel. ) 
Die Formel 


k | 
[yds = 34 (— Yo + 13 y, + 13 ye —yg) (k& = Intervallbreite) (6.5) ji} 


my 
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wurde angewandt zur Berechnung der Integrale 


k 
J ¥ dx = 57 (¥o-9 5, F 19 ye + 943) (6.6) 


bestimmt. Im Bereich ¢ = 2,0 sind im Falle A/a ry = 1 die Funktionen 1— gh bereits nahezu 
linear, so da8 im Rahmen der erstrebten Genauigkeit von der Schrittweite k = 0,2 bedenkenlos zu 
k = 1 iibergegangen werden kann. Am Anfang des Bereichs t > 2 wurde die zu (6.6) symmetrische 
Formel benutzt, anschlieBend wieder (6.5). Damit ergaben sich die Werte von Tabelle 7. 


Tabelle 7. eee! » bzw. Eat Ea bei zeitlich linearer Aufenraumaufheizung fiir is ah 
@ ects 98ers Se: % To 
A 
bd 
bs 1 2 4 8 co 
t 
| 
0,0 0 0 0 0 | 0 
0,2 0,00975 0,00743 0,00742 0,00742 | 0,00742 
0,4 07349 05656 05584 05584 05584 
0,6 1879 1491 1456 1456 1456 
0,8 3341 2736 2653 2652 2652 
1,0 5002 4210 4068 4066 4066 
Le 6789 5846 5642 5638 5638 
1,4 8655 7594 7333 7326 7326 
1,6 1,057 9420 9109 9098 9098 
1,8 oY,252 1,130 1,095 1,093 1,093 
2,0 1,448 1,322 1,283 1,281 1,281 
3,0 2,443 2,305 2,257. 2,254 2,254 
4,0 3,442 3,303 3,253 3,249 3,249 
5,0 4,442 4,303 4,252 4,248 4,248 
6,0 5,442 5,003) ae p22 5,248 5,248 
7,0 OMens oe Oncee (heb Osis 


Graphische Darstellungen der hier tabellierten Funktionen finden sich im Anhang. 


Wie bereits friher erwahnt, ist nach einiger Zeit die Differenz zwischen AuBenraumtemperatur 
und Mittelpunktstemperatur und damit auch 
: ae : =F 
sf ioe! oes vay ee | 
9 eer ogc 
A A 


? 


praktisch konstant. Fiir den beidseitig unendlich langen Zylinder kann diese Konstante aus (5.4b) 


Lk: 
berechnet werden. Sie ergibt sich bei A/arg =1, r=0, B= Ty maz + 1) == 0; 150; 


Tabelle 7 liefert hierfiir den Wert 0,752. Auf Grund dieser Ubereinstimmung dirfen wir annehmen, 
da8 der durch die naherungsweise Integration bedingte Fehler héchstens die letzte Stelle unserer 
Funktionswerte geringfiigig verfalschte. 


Mittels Tabelle 7 bilden wir nun fiir t > 0 die in Tabelle 8 angegebene relative Abweichung 
der Temperaturerhéhung im Mittelpunkt des endlich langen Zylinders von der auf der Achse des 
beidseitig unendlich langen Zylinders bei zeitlich linearer AuSenraumaufheizung. 
= ae 


ze 


36 R. Ansorge: Uber gesteuerte Anheizvorgange bei Zylindern Ingenieur-Archiv i 


Te Ts Dio] 
ao 


Tabelle 8. Ties 
Py, 
1 2 4 8 
t | 
0,2 31,4 0,1 We yr 
0,4 31,6 1,3 Pes is 
0,6 29,1 2,40 = a 
0,8 26,0 Sy lel — —— 
1,0 23,0 3,94 — == 
Be 20,4 3,69 ad i 
1,4 18,1 | 3,66 0,1 — 
1,6 16,2, ~| 3,94 0,12 — 
1,8 14,5 3,39 0,18 —_— 
2,0 13-03-20 0,156 — 
3,0 8,39 2,26 0,133 — 
4,0 5,94 | 1,66 0,123 — 
5,0 AST) |ehe29. 253 ee 
6,0 3,70 1,05 ee / an 
7,0 3,10 0,880 — ae 
8,0 2,68 | 0,759 — = 
9,0 pA 5s 0,667 = — 


Tabelle 8 enthalt wiederum nur Abweichungen iiber 0,1%. 
Der Vergleich von Tabelle 8 mit Tabelle 6 lehrt, daB hier wie dort die relativen Abweichungen | 


der Ubertemperatur T; — T'4 von T) — T4 mit wachsendem i (é > 0) zundchst ansteigen, um 9} 


dann nach Erreichen eines Maximalwertes abzufallen. Dabeiist der Anstieg steiler als der Abfall. 
Dieser Vorgang ist jedoch bei unserem linearen Aufheizprogeramm weiter auseinandergezogen als | 
bei konstanter AuSenraumtemperatur. Die Maximalwerte der Abweichungen bei linearer Auf- | 
heizung sind von gleicher Gré®enordnung wie die bei konstanter AuBentemperatur. Diese Er- || 
scheinungen lassen sich ebenso bei A/a rp = 0,1 und A/a rp = 10 nachweisen. 

Somit bleiben im wesentlichen alle friiheren Betrachtungen hinsichtlich der Méglichkeit, den | 
Zylinder endlicher Lange bei konstanter AuSenraumtemperatur als beidseitig unendlich lang be- 
handeln zu kénnen, auch fiir den Fall der zeitlich linearen AuBenraumaufheizung mit D = Ty, be- 
stehen. Auch der EinfluB der Warmeiibergangszahl bleibt qualitativ erhalten. 


f) Aussagekraft der gewonnenen Ergebnisse fiir andere Aufheizprogramme Q(t). 
Die Aussagekraft unserer unter speziellen Annahmen fiir die Falle O(t) = 0 = konst und O(t) = 
Ot + T, erzielten Ergebnisse hinsichtlich anderer nur zeitabhangiger Aufheizprogramme Q(t) ist 
natiirlich gering. Allerdings werden auch bei diesen Programmen die Abweichungen letztlich durch 


die GréBe der Differenzen 1—h us s—t bestimmt, und der EinfluB der Warmeiber- | 
0 0 


gangszahl diirfte der Tendenz nach ebenfalls erhalten bleiben. Bei der Frage, ob ein Zylinder end- || 
licher Lange durch den beidseitig unendlich langen Zylinder ersetzt werden darf, bleibt daher Vor- 


sicht zumindest bei den Parameterkombinationen (x, A/a ry) geboten, bei denen auch die in dieser ||} 


Arbeit untersuchten Falle starke Abweichungen von den Lésungen der entsprechenden ebenen |) 
Probleme zeigten. 


g) Vergleich des Falles x =1 mit dem der Kugel. Gelegentlich (z. B. bei differential- | ; 


thermoanalytischen Untersuchungen?) interessiert der Vergleich des Zylinders Z, mit der diesem |} 


Zylinder einbeschriebenen Kugel. 


Die Lésung von (3.1) fiir die Kugel lautet im Falle konstanter AuBenraumtemperatur % bei im |) 


Innern der Kugel iiberall gleicher Anfangstemperatur T 4 


és oy sin|@, — Lae we 
i ees? > en 1 2 sin @, ( re ey 
— s. : i) 
o— Dy 7 A (@, (oe . Seah 4 ; 
= = — $1n @, COs W — 
a aT, n n A @, To 


1 Lehmann-Hafler, TIZ-Zbl. 82 (1958), 445. 
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fiir den Mittelpunkt (r = 0) also 


Lite I ig oe = 2 sin © — 2 ay 
= ae a - es aS (6.7) 


ae (©, — Sin @,, Cos @,,) 


Dabei sind die Eigenwerte @, die positiven Lésungen der Gleichung 
OCC Os: (6.8) 


Es geniigt, (6.7) mit dem Fall x = 1 der entsprechneden Formel (6.1) fiir nur einen Wert des Para- 
meters A/a rg zu vergleichen. Wir wahlen der Einfachheit halber A/a ry = 1, so da’ die Lésungen 
von (6.8) lauten 


Dee, (eS. tsi) 
(die Lésungen fiir andere Werte von A/a r, finden sich bei Gréber-Erk-Grigull *). Damit folgt 
<0 _ @n—1a- 


T—T, _ 4 (—1r-1 —>— 
= ee Da Te (6.7a) 


(fiir ¢ = 0 erhalt man mit der Leibnizschen Reihe natiirlich T = T,). 


Uber unserem schon friiher benutzten t-Raster ergibt (6.7 a) die Werte von Tabelle 9. Hine gra- 
phische Darstellung dieser Funktion ist im Anhang beigefiigt. 


T—T Pc A 

Tabelle 9. rene T, fir die Kugel, ay == Il 
0,0 0 1,6 0,9754 
0,1 0,05071 1,8 9850 
0,2 2277 2,0 9909 
0.4 5255 

0.6 7103 3,0 9992 
0,8 8231 4,0 1,0000 
1,0 8920 5,0 1,0000 
1,2 9341 - 

1,4 9598 


Mit den entsprechenden Werten der Tabellen 3 bis 5 bilden wir nun die relativen Abweichungen 
der Temperaturerhéhung der Kugel von der des Zylinders Z,, also die Werte von 


(Txuger ii T 4) aie (Tz, ire T 4) 


0 
= 100 [%] 


(aufgetragen wurden nur Werte iiber 0,1%) in Tabelle 10. Die Werte zeigen, da sich der Zylinder 
Z, beziiglich der Erwarmung zunichst viel trager verhalt als die Kugel von gleichem Radius. Die 
Abweichungen fallen jedoch bedeutend schneller als z. B. die Abweichungen der Mittelpunkts- 
temperatur in Z, von der in 2. 


* Siehe FuGnote 1 von Seite 30, 
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(Tugel a T 4) aay (Tz, ay T 4) 


100 [%] fiir 


10. 
Tabelle Ty Seyi, 
i t 
0.1 91 1,6 0,879 
0,2 31,8 1,8 0,603 
0,4 12,6 2,0 0,395 
0,6 7,02 
0,8 4,40 3,0 — 
1,0 2,91 4,0 — 
12 1,94 5 
1,4 1,31 é 


7. Anhang: Tabellen und graphische Darstellungen. 


2 a3 
Tabelle Al. g Es a0 ‘ 
0 


0 1,0000 1,0000 1,0000 3,0 000001 01065 | 5707 
0,1 0,9000 0,9768 0,9974 4,0 002199 | 4695 
0,2 6002 8702 9837 5,0 000454 | 3863 
0,4 2344 6420 9476 6,0 000094 +| 3179 
0,6 090 67 468 6 9114 7,0 000019 2616 
0,8 03506 3419 8765 8,0 000004 2152 
1,0 01356 2494 8430 9,0 1770 
1,2 005 244 1819 8108 10,0 1457 
1,4 002028 1327 TOT 11,0 1199 
1,6 000 784 09682 7499 12,0 098 62 
1,8 000303 070 63 7211 13,0 08113 
2,0 


000117 05152 6936 


Tabelle A 2. h(x: A nae i) 
aT 


0 1,0000 | 1,0000 | 1,0000 | 1,0000} 1,0000  |.1;0000 | 1,0000 | 1,0000 | 1,0000 | 1,0000 1,0000 | 1,0000 
0,1 | 0,9684 | 1,0000 | 1,0000  1,0000} 0,9931 | 1,0000 | 1,0000 | 1,0000 | 0,9992 | 1,0000 1,0000 | 1,0000 
0,2 8293 | 0,9980 | 1,0000 | 1,0000 9506 | 0,9996 | 1,0000 | 11,0000} 9940 | 0,9999 | 1,0000 | 1,0000 
0,4 5574 9610 | 0,9999 | 1,0000 8309 9879 | 1,0000 | 1,0000| 9772 9984  1,0000 | 1,0000 
0,6 3707 8877 | 9994 | 1,0000 T1717 9581 | 0,9999 | 1,0000} 9587 9943 | 0,9999 | 1,0000 
0,8 2464 8043 | 9972 1,0000 6190 9179} 9992/1,0000} 9404 | 9883 9998 | 1,0000 
1,0 1638 7230 | 9919 | 1,0000 5339 8733 | 9972/1,0000] 9224 | 9810; 9996 | 1,0000 
1,2 | 1089 6479 | 9830 | 1,0000 4604 8278 | 9935} 11,0000] 9047 9730| 9991 | 1,0000 
1,4 07239 5798 | 9705 | 1,0000 3971 7831} 9877/1,0000] 8874 | 9646) 9982 | 1,0000 
1,6 04812 9186) 9551 | 1,0000 3424 7399 | 9802/1,0000] 8704 | 9560! 9970 | 1,0000 
1,8 03199 4638 | 9374 | 1,0000 2953 6987} 9709|1,0000} 8537 | 9473) 9955 | 1,0000 
2,0 02126 4147} 9180 | 0,9999 2547 6596 |; 9601|1,0000} 8373 | 9387] 9936 | 1,0000 


3,0 00276 2370} 8100) 9981 1215 4938 | 8928/0,9991] 7601 8959 | 9807 | 0,9999 
4,0 : : : : 05795 3695; 8171) 9957] 6900 | 8549| 9639 9993 
5,0 02764 2766| 7429) 9883] 6264 | 8158| 9453 9979 
6,0 01318 2070| 6736) 9768} 5686 | 7785} 9260 9956 
7,0 00629 1549} 6100; 9618] 5162 | 7428| 9064 9923 
8,0 00300 1159| 5522) 9441] 4686 | 7089! 8870 9881 
9.0 F 2 : 4254 | 6764) 8679 9831 
10,0 3862 | 6455| 8491 9774 
11,0 3505 | 6159} 8307 9712 
12,0 “ 3 3182 5878 | 8127 9645 
13,0 | | : 5 2889 | 5609| 7950 9574 
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Abb. 1. Konstante AuSenraumtemperatur. 
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Abb. 2. Konstante AuGSenraumtemperatur. 
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Abb. 3. Konstante Aufenraumtemperatur. 
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Abb. 4. Zeitlich lineare AuBenraumtemperatur. 
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Abb. 5. Konstante AuSenraumtemperatur. 


8. Zusammenfassung. Ziel der Arbeit war die Untersuchung des Einflusses der endlichen 
Zylinderlange und der Warmeiibergangszahl auf die Temperaturverteilung in Vollzylindern bei ge- 
steuerten Anheizvorgangen. Fragen dieser Art kommen z. B. in der Me8technik (EinfluB der Ge- 
stalt eines Thermometerkérpers auf die Temperaturmessung, Untersuchungen von Erden und dgl. 
mittels Differentialthermoanalyse) und in der Reaktortechnik vor. 


Von der Temperatur des den Zylinder umgebenden Mediums wird angenommen, daB sie nur von 
der Zeit t abgingt. Die Zylinderanfangstemperatur habe in allen Punkten den gleichen Wert. Fir 
die Mittelebene des Zylinders wird gezeigt; daf& die Lésung des dreidimensionalen Problems bei 
wachsender Zylinderlange gleichmaBig in t gegen die Lésung des ebenen Problems konvergiert. 


In den Fallen konstanter Umgebungstemperatur und zeitlich linearer AuSenraumaufheizung 
werden fiir verschiedene Kombinationen der Parameter A/a ry und x die Temperaturkurven fiir den 
Zylindermittelpunkt erstellt und die relativen Abweichungen dieser Temperaturkurven der Zylinder 
endlicher Lange von denen des beidseitig unendlich langen Zylinders errechnet. 
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Aus den Temperaturkurven fiir verschie 
ubergangszahl erkennbar, 
beigegebener zulassiger 
werden darf. 


dene A/a ry bei festem x ist der Einflu8 der Warme- 
aus den Tabellen der Abweichungen die kleinsten Werte von x, fiir die 


Abweichung und gegebenem A/a Ty die Aufgabe als ebenes Problem behandelt 


Fir den Fall ~ = 1 bei konstanter AuBenra: 
linder einbeschriebene Kugel bezii 


Zylinder, daB jedoch alsbald die 


umtemperatur wird gezeigt, daB sich die dem Zy- 
glich der Erwarmung zunichst viel weniger trage verhalt als der 
Abweichungen gering werden. 


(Eingegangen am 1, Februar 1960.) 


Anschrift des Verfassers: Dr. R. Ansorge, Institut fiir Mathematik und Mechanik d 


er Bergakademie Clausthal, 
Clausthal-Zellerfeld. 
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Die Knickung der tordierten Welle mit Einzelkraft 
und kontinuierlichen Langskraft* 


Von H. Leipholz 


1. Einleitung. Das Auftreten einer kontinuierlichen Lingskraft scheint, so weit es mir bekannt 
ist, fiir die tordierte und gedriickte Welle noch nicht beriicksichtigt worden zu sein. Mit diesem 
neuen Problem soll sich daher die vorliegende Arbeit befassen. 

Bei den hier betrachteten Staben soll es sich um ,,lange, dine“ Stabe handeln. Jeder Stab- 
querschnitt habe zwei gleiche Haupttrigheitsmomente, die zudem iiber die Lange des Stabes kon- 
stant seien. 

Belastet sei der Stab durch kontinuierliche, tiber seine Achse nach einem vorgeschriebenen 
Gesetz verteilte, richtungstreue Druckkriafte q, durch eine ebenfalls richtungstreue Einzelkraft P 
und durch ein axiales Torsionsmoment T, welche an den Stabenden angreifen. 

Ferner wird verlangt, daB die Torsion konservativ sein soll. Damit ist dann das ganze Problem 
konservativ; denn bei den richtungstreuen Druckkraften handelt es sich von vornherein um kon- 
servative Krafte. Ist das Problem aber konservativ, so kann man zur Ableitung der Differential- 
gleichungen die statischen Gleichgewichtsbedingungen anwenden. Dies wurde bereits von H. Zieg- 
ler! gezeigt. 

Soll das an den Stabenden angreifende axiale Torsionsmoment konservativ sein, so mu man, 
wie H. Ziegler gezeigt hat, als Stablagerungsfalle die in Abb. 1 gezeigten verlangen, was im folgenden 
beriicksichtigt wird. 


2. Die Grundgleichungen. Wendet man auf ein herausgeschnittenes Stabelement die Gleich- 
gewichtsbedingungen der Statik an, so erhalt man durch die Forderung nach dem Gleichgewicht der 
- Krafte 


d& : Abie 
) Sees f (1) 
und durch die Forderung nach dem Gleichgewicht aller Momente 
dy) 
ot + (tx 8) =0. (2) 


In diesen Grundgleichungen (1) und (2) bedeutet Jt die Resultierende der inneren Momente, & die 
Resultierende der inneren Krafte in einem Querschnitt des Stabes, ferner { die 4uBere Kraft und 
t den Tangenteneinheitsvektor fiir den Punkt der Stabachse, durch den der betrachtete Querschnitt 
geht. 


3. Das Hauptachsensystem und seine Anwendung auf die Grundgleichungen. Als Achsensystem 
soll insbesondere das orthogonale Hauptachsensystem (&, 7»¢) eingefiihrt werden. Diel-Achse zeigt 
in Richtung der Stabachsen-Tangente im Sinne wachsender Bogenlange. Die é- und y-Achsen 
liegen in einer zur Stabachse senkrechten Ebene, die auch Ebene der Stabquerschnitte ist, und sie 
sollen mit den Schwerpunktshauptachsen der jeweiligen Stabquerschnitte zusammenfallen. 

Die Komponenten der in (1) und (2) auftretenden Vektoren seien im Hauptachsensystem 


@ =(0,0,N),  M=(MyM,, 7), 
f = (keg, f,, q)> £—A050, Eis 


Ferner ist noch zu beachten, da8 das Hauptachsensystem beweglich gedacht ist: es gleite mit der Ge- 
schwindigkeit Eins an der Stabachse entlang und drehe sich dabei mit der Winkelgeschwindigkeit tv. 
Im Hauptachsensystem nimmt (1) und (2) die Form 


d’'S 

qt wx R+t=0, (3) 
d'it 
az tw x M+t+H=0 (4) 


: Auszug aus einer von der Technischen Hochschule Stuttgart genehmigten Dissertation (1959). 
H. Ziegler, Z. angew. Math. Physik 2 (1951), S. 265; 3, (1952) S. 96; Ing.-Arch. 20, S. 49 (1952). 
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an, wobei das Symbol d’/ds die Ableitung einer Gré®e im bewegten Hauptachsensystem bedeutet und 
-W = (@g, @,, Te) ist. : 


Schreibt man (3) und (4) in Komponenten, so erhalt man die ,,erweiterten“ Gleichungen von 


Kirchhoff und Clebsch 


dQ: dM. 

tt Ne Oth, = 0; to, — Mc Op 0, ) 

dQ dM, | 

ae + Ost — Nz wy +k, = 0, —, + Mpt— Taz +Q;=0, $ (5) 
dN, dT | 

ae + On Oe — Og O, +9 =0, =, + M,o¢— Mz, =0. | 


4. Zusammenhang zwischen den Gréfen der Stabverformung und den inneren Momenten. Wie 
bekannt, besteht zwischen den Komponenten des Vektors tv (Drehgeschwindigkeit des Haupt- 
achsensystems) und den Komponenten des Vektors der elastischen Verformung ein derartiger 
Zusammenhang, dah Wg, ®, den Stabkriimmungen und t; der Verdrillung entspricht. AuSerdem 
ist bekannt', da zwischen den Komponenten des Verformungsvektors und den Komponenten des 
Stabmomentes die Beziehungen bestehen 

Mz;=aoz, M,=ao,, T=crt,. (6) 
Dabei sind a die konstanten Biegesteifigkeiten und c die konstante Torsionssteifigkeit des Stabes?, 

Die sechs Gleichungen (5) und die drei Gleichungen (6) sind zur Bestimmung der neun Unbe- 

kannten dieser Gleichungen (Komponenten von 9)¢, § und {v) ausreichend. 


5. Der Zusammenhang zwischen dem Hauptachsensystem und einem raumfesten Koordinaten- 
system. Es ist zweckmaBig, eine Beziehung zwischen dem Hauptachsensystem und einem raum- 
festen (absoluten) Koordinatensystem herzustellen. Die elastischen Konstanten und die inneren 
Krafte des Stabes lassen sich im Hauptachsensystem wohl leicht darstellen. Dies wird aber haufig 
fiir die Randbedingungen und die auBeren Lasten nicht gelten. Diese werden meist beziiglich eines 
raumfesten Koordinatensystems gegeben sein, so das man sie auf das Hauptachsensystem umrech- 
nen mu. 

Sieht man von der belanglosen gleichférmigen Translation ab, so kann man das Hauptachsen- 
system als ein um einen festen Punkt kreiselndes Dreibein ansehen. Die Drehgeschwindigkeit dieser 
Kreiselung ist durch den Vektor tv ausgedriickt. In dem festen Punkt mége auch der Ursprung des 
raumfesten Systems liegen. Die Herstellung einer Beziehung zwischen dem kérperfesten (Haupt- 
achsen-) und dem raumfesten System ist mit Hilfe der Eulerschen Winkel moglich. Fir den Ubergang 
vom raumfesten Koordinatensystem zum Hauptachsensystem hat man die Transformationsmatrix 


cos y cos p — sin y sin g cos } sinycosy +cosysingcos} singsin? 
=| —cosysing—sinycospceos? —sinpsing +cosycospcos# cospsind|. (7) 
sin yp sin 3 — cosy sin 0 cos ? 


Die Beziehung zwischen den Komponenten des Drehvektors {y und den Eulerschen Winkeln wird 
gegeben durch . 


W, =cosp i +singsin dp, O, =—sing J +cospsin Fy, T=cosdyp+o. (8) 


Die zum ersten Eigenwert des Problems gehérende Gleichgewichtslage des Stabes wird nur 
wenig von der Ruhelage des unverformten Stabes verschieden sein. Daher wird # ein sehr kleiner 
Winkel sein, so daf man linearisieren und cos # = 1 setzen kann. Somit ergibt sich aus der dritten 
Gleichung (8) der neue Winkel 

(Ce aM 
also 


ety=xaf[rds. 3 : (9) 


1 4, E. H. Love, A treatise on the mathematical theory of elasticity, S. 388. Cambridge 1952. 


2 Zur Verwendung des Hauptachsensystems bei kreisrundem Querschnitt ist zu bemerken: Man markiere 
_ vor der Verformung des Stabes ein (&, 7, ¢)-Kreuz mit einem beliebigen orthogonalen (&, 7)-Achsenpaar, da hier 
ja ein jedes (£, 7)-Paar Hauptachsenpaar ist; dann verfolge man es in Gedanken wahrend der Verformung, so 
da8 man es nach der Verformung wiederfindet. Das auf diese Weise ausgezeichnete Achsensystem macht, 
wie es vom Hauptachsensystem verlangt wird, Kriimmung und Verdrillung des Stabes mit. 
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Beachtet man die Kleinheit des Winkels @ und verwendet man (9), so wird aus (7) 
cosy siny @ 
A, =| — sin CO ye, (10) 


0 0 1 


Dabei hat man noch die Tatsache beriicksichtigt, daB fiir einen lotrechten Einheitsvektor ¢e, der 
im absoluten System die Komponenten (0, 0, 1) und im Hauptachsensystem die Komponenten 
(ez, p> €¢) hat, gilt 
e.g =sing sin? , e, =cospsin id, e =costxl, (11) 
was man erkennt, wenn man auf ¢ die Transformationsmatrix (7) anwendet. 

Kiinftig kann dic Umrechnung der Komponentendarstellung von Vektoren aus dem absoluten 
' Koordinatensystem in das Hauptachsensystem nach 
qHuaAs — YW, abs (12) 


erfolgen. 


6. Aufstellung der Differentialgleichungen. Die Vektorgleichungen (1) sollen zunachst in den 
Komponenten des absoluten, raumfesten (x,y, z)-Systems geschrieben werden, da sich ihre Lésungen 
dann leicht angeben lassen. AnschlieBend sollen diese Lésungen mit Hilfe der Transformations- 
matrix Y{, im Hauptachsensystem dargestellt werden. 


Im (x, y, z)-System hat man die Komponentendarstellung 


K = (Q,, Qys ine t= 40, 0, q) « 
Das gibt fiir (1) 


dQ, AQy __ dN, _ 
Sp POE = it cay meek (13) 
Hierfir lauten die Lésungen 
Q., A = konst. 
Ghabs = Q, ae ‘Be konst. : (14) 
l 
Ne — P—f qds 


und aus R445 — YI, Ms folgt 


1 
A cosy + Bsiny—e(P + f qds) 


/Q 
é 1 
QHAS = | Q, }= | — Asing =P B cos ye. (P -pafogeds). | (15) 
Ne ! 
— P—f qds 
Die Vektorgleichungen (4) lauten in den Komponenten des Hauptachsensystems 
dM: 
gs Mn te + TO, = Q,» | 
dM, 
a +t Mey —To=—Q, (16) 
dT 
qe + Mn eg — Mz, = 0. 


Setzt man hierin (6) und (15) ein und ordnet um, so hat man 


dw } 
a + C—a)to, +¢,(P + fqds)=—Asiny + Boos z, 
_ diy : 
age 4) tee — 0 (P+ fads) = —A cosy — Being, (17) 
c dv =n () es 


“ds” 
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Die dritte Gleichung (17) kann sofort integriert werden und gibt 


c = konst., 


CT; = konst., poe ote = T. (18) 
= nst. 


Jetzt sollen statt der GréBen @zg und w, Ausdriicke eingefiihrt werden, die nur ee, e, und ihre 
ersten Ableitungen enthalten. Nach (8) ist 


Og =cosy)+sinpsin dy. 


Wegen der Voraussetzung cos ? = 1 kann man beide Glieder der rechten Summe mit cos 3 multi- 
plizieren und erhalt 


(oe = cosy cos 99 + sin psin# cos Pp. 


Jetzt werde auf der rechten Seite der Gleichungen noch — sin gp sin 9 und + sin g @ sin # hinzu- 
gefiigt, so daB sich nichts andert. Dann kann man anders zusammenfassen: 


Wg = — sing Pp sin & + cos @ cos OO+ sin y sin } (cos dy + ¢). (19) 
Aus (11) entnimmt man 
ee =singsin?, e, =cosgysind, 
also 
é. = —singgsin? + cosq@ceos 03 
und aus (8) : se 3 


tT =cosdy+o. 
Vergleicht man dies mit dem rechten Teil von (19), so sieht man, daB 

On — e + Cz Te (20 
ist. Ganz entsprechend findet man e 
Wy, = — @g + €, Te - (21) 


Mit Hilfe von (20) und (21) kann man die ersten beiden Gleichungen (17) anders schreiben. Man 
beachte auBerdem, dais wegen der Konstanz von t, gemaf (18) 


x ail Te ds = s (22) 


wird. Man benutze die Abkiirzungen 


und schreibe noch fiir die gegebene GréBe 
fqds =Q(s); (23) 
so erhalt man statt (17) 
? en ateat he 2 P, A) ats ; 
e, +(2—4)e és + (e—1)% ctce | =— A*sint,s + B* costs, 
(24) 


ie—[0— Sets + (G2) + E+ Df m arom ge tBraing, 


Das ganze Problem ist damit auf die Lésung des linearen Differentialgleichungssystems (24) von 
zweiter Ordnung, mit variablen Koeffizienten bei e, bzw. e,, zurtickgefthrt. 
Zu seiner vollstandigen Formulierung miissen nun noch die Randbedingungen angegeben werden. 


7. Randbedingungen. Die Randbedingungen ergeben sich aus der Art der Stablagerung. Wie 
man aus Abb. 1 abliest, sind folgende Falle zu beriicksichtigen: 


a) Die Stabenden sind frei von Querkraften. Das ist der Fall, wenn die Stabenden in der Hori- 
zontalen frei verschieblich sind. Aus (14) ersieht man, daB dies bedeutet 


A=0, B=0. (25) 


Dann werden die Gleichungen (24) homogen. 
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b) Die Stabenden sind fest eingespannt. Im Falle einer Einspannung fallt am Stabende der 
lotrechte Einheitsvektor ¢ mit dem Tangentenvektor t zusammen. Dieser hat aberim Hauptachsen- 
system die Komponenten (0, 0, 1), und da jetzt e = t sein soll, so verschwinden die ersten beiden 
Komponenten des Vektors e: 
e. = 0, &, = 0. (26) 

c) Das Stabende ist unverschieblich, oder es ist in Richtung der urspriinglichen Stabachse ge- 
fiihrt. In diesem Fall mu der Vektor zy, = (x, y) fiir s = I, der in der durch den Stabanfang gehen- 
den Horizontalebene liegt und der auf ein raumfestes (x, y, z)-“System mit dem Ursprung im Stab- 
anfang bezogen ist, gleich Null sein. Nun ist aber 


s 
7, Up = 1% ds 
mit 
dx dy ie 
er) eee | ecteeieet L6) Ve a 
trea (a a) = age 
a 
Der Vektor t#45 — — (eg e,) werde mittels der Matrix 2[, umgeschrieben: 


tb = (x, 7) = UT tHAS, 


Daraus folgt 


Abb. 1. Lagerfall 1 und 2. . = : bs 
x = — (e cos y — e, sin Y) , y = — (esin x + e, cos x). 


Das gibt als Randbedingung ; 
1 
J (@ cos ¥ —e, sin y) ds = 0, 
0 


l 
bp = [tp ds = 05 1 (27) 
: J (@sinz +e, cosy) ds =0. 
f 


8. Die Lésung der homogenen Differentialgleichungen. Es werden die homogenen Gleichungen 
(24) mit der Abkiirzung 
CeO) 
aS G + oe 
geschrieben: 


sg) s[eajedune 


ie—[-—S)ntg+|(E—2) 2 o]a—0, 


und es kann gezeigt werden, da8 sie formal mit den Bewegungsgleichungen eines Punktes auf einer 
Kugelflache in einem zeitlich veranderlichen Kraftfeld mit dem Feldvektor 


. Bs = —m Vit) e, (29) 
iibereinstimmen. Es soll tim folgenden die Zeit bedeuten. 
Die Bewegung werde beschrieben in einem Achsenkreuz (£,77,¢), das seinen Ursprung im Mittel- 
punkt der Kugel hat, dessen¢-Achse durch die Pole der Kugel geht, und dessen (€, 7)-Ebene um die 
¢-Achse gleichférmig mit der Drehgeschwindigkeit 
u = (0, 0, u) 


rotiert. Ist R = 1 der Halbmesser der Kugel, so lauten die Lagrangeschen Bewegungsgleichungen 
fiir das raumfeste Inertialsystem (x, y, z) 


(28) 


MX == 1A.%, my =miy, mz=—mV(t)+mdz, 
wenn : 


, F=>(?@+y4+2—1) =0 
die Zwangsbedingung fiir die Bewegung auf der Kugel und / ein noch unbestimmter Parameter sind. 


Um diese Gleichungen fiir das bewegte (&, y,¢)-System umzuschreiben, mu8 man von den Ge- 
setzen der Relativbewegung Gebrauch machen und erhalt in vektorieller Form 


mb, = 8 —2mu xX v,—mu X (ux rt) + Agrad F. (30) 
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Dabei soll auBer den bereits erklarten GréBen noch b, der Vektor der Relativbeschleunigung, x der 
Ortsvektor in der (§, 7)-Ebene, welcher senkrecht zu 1 ist, und b, der Vektor der Relativgeschwin- 
digkeit sein. Schreibt man (30) in Komponenten, so gibt das 


d? d 
mga mug + mew + mdE, | 
d* dé 
mama sent cman i 
2 


Es ist aber d?¢/di? = 0 und ¢ ~ 1, wenn nur solche Bewegungen betrachtet werden, bei denen 
€ und 7 kleine GréBen bleiben. Diese Annahme ist im Kinklang mit der beim Stab gemachten 
Voraussetzung kleiner Verformungen. 

Dann folgt aus der dritten Gleichung (31) 


A=; 
und damit lassen sich die beiden ersten Gleichungen (31) schreiben 
—&—2uy—(uw+V)E=0, 7 +2ué—(W+V)y =0. (32) 


Die Gleichungen (32) stimmen mit den Gleichungen (28) genau iiberein, wenn 

u =(l —¥4)% , und Fareealas| 
ist. . : 
Von der Relativbewegung des Punktes auf der Kugelflache interessiert die Projektion der Be- 
wegung auf die (€,7)-Ebene. Die Bahnkurvengleichung ist beziiglich des (£, 7)-Systems zu bestim- 
men. Hat man sie gefunden, so hat man mit Riicksicht auf die Ubereinstimmung zwischen (32) 
und (28) auch die Lésung fiir die homogenen Differentialgleichungen des Stabes erhalten. 

Die Gleichungen (32) lassen sich mittels 9 = € + i7 komplex zusammenfassen zu 
do 


ead 


Im raumfesten (x, y)-Kreuz sei der Ortsvektor r = x + iy erklart. Dann ist 
CV Gan (34) 


die Beziehung zwischen den beiden Darstellungen des Ortsvektors im bewegten (€,7)- und im raum- 
festen (x, y)-Kreuz. 
Setzt man (34) in (33) ein, so erhalt man 


For 0; (35) 


was die Bewegung des Punktes im raumfesten (x, y)-System beschreibt. 

Aus (35) folgt: Der Punkt bewegt sich mit seiner Projektion auf dem in der (x, y)-Ebene fest- 
liegenden Ortsvektor r. Die Bewegung gegeniiber der rotierenden (£, 1)-Ebene wird aus der Lé- 
sung von (35) erhalten, wenn man gema8 (34) mit e—*“' multipliziert, was der Uberlagerung einer 
Rotation mit u entspricht. Der im (,7)-System feste Beobachter wird also die Projektion des 
Punktes als Bewegung auf dem Fahrstrahl g sehen, der seinerseits mit der konstanten Drehge- 


_ schwindigkeit 4 um den Ursprung rotiert. 


Die durch (35) festgelegte Bewegung auf dem Fahrstrahl wird nun je nach der Art von 
P,Q) ete\? 
Se tae) 


a a 


- verschiedener Natur sein. Die Gré®en P/a und c7,/2 a sind Konstante. Die Art von V ist also im 
 Wesentlichen durch Q(t) bestimmt. Ist Q(t) im Sonderfall auch konstant, so ist V = konst. und 
_ (35) beschreibt eine harmonische Bewegung. Das gilt natiirlich inshesondere fir Q(t) = 0. Wegen 


der dann auftretenden harmonischen Bewegung der Projektion des Punktes zeigt es sich, dafs der 


~~ Punkt in diesem Sonderfall sich wie der eines Foucaulischen Pendels verhalt. 


i allie 


Im allgemeinen Fall wird Q(t) weder Null noch konstant sein. Dann stellt (35) eine lineare 
Differentialgleichung zweiter Ordnung in der Normalform mit veranderlichem Koeffizienten V dar. 
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Die Differentialgleichung (35) habe die partikuléren Lésungen : 
rs 7, (2) (k=, 2) 
Dann findet man die Lésung von (33), indem man mit e—*“* multipliziert: 
0, = (t) e—*** (h = 1,2). 
Die allgemeine Lésung von (33) ist also 
e=E+in = [Ur (t) + Br] eo (36) 


mit den beiden komplexen Integrationskonstanten 2 und B. 
Mit (36) hat man sowohl die Projektion der Bahn des Punktes im bewegten System als auch die 
Lésung fiir die homogenen Differentialgleichungen des Stabes erhalten. Man braucht fiir den Stab 


nur iiberall statt t die Variable s zu schreiben und 


5 c\te 
_— (2 —_ a > 
zu setzen. 


Wegen ihrer grundsatzlichen Bedeutung fiir die Lésung des Problems soll (35) die Basisgleichung 
heiBen. 


9. Die Eigenwerte fiir Lagerfall 2 und q(s) = 0 bzw. q(s) = q = konst. Im Lagerfall 2 von Abb. 1 
sind die Stabenden frei von Querkraften, so da dic Randbedingungen (26) gelten. Die Differen- 
tialgleichungen des Stabes werden daher homogen und man hat mit (36) die Lésung schon zur Ver- 
fiigung. 

a) q(s) = 0. Das gibt fiir die Basisgleichung 


d’r Je is ee 
Schreibt man zur Abkiirzung 
Pee eceye 
Fy (eef aut i 


so lautet die Lésung von (37) 
ry = 8IN 09S, Tz. = COS VyS. 


Bezugnehmend auf die Analogie zum bewegten Punkt, kann man sagen, daB die Punktprojektion 
im raumfesten (x, y)-System eine lineare, harmonische Schwingung auf dem Fahrstrahl ausfiihrt. 
Beziiglich des mit u rotierenden (&, 7)-Systems liegt also die bereits erwahnte Foucaultsche Pendel- 
bewegung vor. 

Die allgemeine Lésung von 0(s) wird zu 


o(s) = ee +e, = [sin v9 s + B cos vy s] e—**s. (39) 
Auf (39) sind zur Bestimmung der Eigenwerte die Randbedingungen (26) in komplexer Schreib- 


weise anzuwenden: 


@=e tie, =—0 fir.s=0 und s=l, (40( 


was zu der Forderung nach der Singularitat der Matrix 


ea 0 1 
tk . * ° 
cin, te" cos, 1°?! 


also det€ = 0 fiihrt. Dies gibt 
sin vole; oe 0. 


woraus unmittelbar Fe Eigenwertgleichung 
snvyl=0, also nl=a2 (41) 


folgt. hao man noch (38) und ¢ t; = T , so hat man sofort die bekannte Greenhillsche Formel 
vor sic 
1 Ne I \? 
£+(8)-( 
b) q(s) =q =konst. Die Basisgleichung lautet jetzt wegen Q(s) = f qds = q(l—s) 


P q(l—s) 2 
mae lg a +(55) [r= 9. (42) 
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Nun seien die Abkiirzungen eingefiihrt 


Meee, Tt ct\? q 
eae ete s\ Par: (43) 
womit (42) in 
d*r 
+ (a—fs)r <0 (44 
tbergeht. Wendet man auf (44) die Transformation 
g = Peep s (45) 
an, so erhalt man 
dr 
Po +gr=0. (46) 


Aus (46) 1a8t sich mittels einer Lommel-Transformation eine Besselsche Differentialgleichung 
herstellen, so daB man fiir (46) die beiden linear unabhangigen Lésungen 


2, 
i= 37 Cys (s ia : (eee Ueny ie 1/3 ts i") (47) 


erhalt. Die allgemeine Lésung fiir die Differentialgleichungen des Stabes lautet daher in diesem 
Fall 


os) = 6 Pe = ee 3M? Jays i”) +B 3M? Japs gimjenins. (48) 


Auf (48) sind wieder die komplexen Randbedingungen (40) anzuwenden. Fiir s = 0 und s =I 
hat man nach (43) und (45) die Beziehungen 


cao yeas [EA (AYIA, 


a a a 


+ 


Hiermit geben die Randbedingungen wieder die Forderung nach der Singularidt der Matrix 
2 2 
ay" Jun (3 32”) ay" Fan ($0) 


2 —iu 2 —iu 
ay Jiys (= ai”) é 1 au? J_1p3 (5 at”) é l 


(49) 
aes a= 6 Pp. FB — pet = 


Setzt man ihre Determinante gleich Null, so findet man als Kigenwertgleichlung 
2 2 2 2 
Jijs & ai") Jz is (= at”) = Jiys & at”) J= 1s & | : (50) 
Man kann nun noch die Funktion 


__ Jis(2) 
tb = 1/3 51 
ctb 1/3 (z) J-1) (51) 
-einfiihren. Ihre Reziproke hat fiir groBe Werte des Argumentes z einen tg-ahnlichen Verlauf. Ihr 
Funktionsbild ist in Abb. 2 dargestellt. 

Man kann die Kigenwertgleichung (50) mit Hilfe von cthj)3(z) schreiben 


2 de 
cthi/3 & °) — cthi/3 tz a”) (52) 


/ Die Ermittlung des kleinsten Eigenwertes geht daher folgendermaSen vor sich: Gegeben seien 
-T, P, q, a. Auf Grund von (52) erhalt man die zu diesen Daten gehérende Knicklange his, des 


; 04 , 
‘Stabes. Denn die gegebenen GrdBen gestatten es, gemaB (49) den Wert 4 = 3; zu errechnen. 


Die weitere Lésung sei graphisch durchgefiihrt. Man sucht auf der Abszissenachse des Funktions- 
bildes (Abb. 2) den errechneten Wert % auf. Zu diesem 4 liest man die Ordinate etbj/3(z,) ab. Im 
Abstand dieser Ordinate zieht man eine Parallele zur z-Achse. Es interessiert der erste rechts von 
| z, gelegene Schnittpunkt der Parallelen mit dem betreffenden Zweig der Funktion ctbj/3(z); denn 
= 4 


pear \ faa 


a —- os 
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ees 2 35s + ew | 
die zu diesem Schnittpunkt gehérende Abszisse ist % = = 33/2. Mit Hilfe des soeben gefundenen % || 


kann man dann nach (49) ausrechnen: 


wen BSE ATIGT” | 


Zwei Zahlenbeispiele: 
1) T=10!kpm, a=10'kpm’, g= 10?kpjm, P=0, 
2) T=10*kpm:, .@ = 10° kpm?,- ¢ = 60 kp/m, P=10% kp. 


ctby ile) 4 a 


Abb. 2, 


Aus der zweiten Gleichung (49) erhalt man 


3/2 
108 » 108 — 0,008 , 


yi 
1) 4= 3 a7 = (710m -1020) 03 


a 
8 3\ 3/2 5 
ae =F is) ag 15S 


2) a= - ai 4.1020 | jos 


Hierfiir findet man graphisch 
1) 2% 2,92 , 2) % & 4,66, 
und es ist nach (53) 


Pe te 192)2/3 2) 1/3 
i) bw. = |(3 2,92)" — oan (10%) | 209 este Gy 


ayia \( 4,66)" — (1,67 - 108)2/3 im =) (1,67 - 10°) — 29,5 m. 


Ware die Langskraft nicht gleichmaBig verteilt, sondern konzentriert als Einzelkraft an den|| | 
Stabenden angebracht: 


Lv P=2 10k, 9 2 P S60: 22,5 — 1350 kp , 


so ware die hierfiir entsprechende Knicklange nach der Greenhillschen Formel zu ermitteln. Sie | 
wiirde dann nur i) 


1) Unrit, = 26 108 (108 + 4-105. 12- 108)-¥? = 9m, 
2) _ Leis, = 2 7% 105 (108 14.105. 2,35 - 10%)— 2 — 19,5 m 
betragen. ; 
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An Stelle von (53) kann man noch eine sehr brauchbare Naherungsformel angeben. Es sei 
% = € + %, wo € ein veranderlicher Parameter ist. Nun ist aber fiir die als Abszissen von cth1/3(z) 
auftretenden Wertepaare z, und z, der Parameter ¢ stets sehr nahe bei 1. Dies gilt insbesondere fiir 
groBe Argumente z. Es la8t sich namlich zeigen, daB cth;/3(z) fiir groBe Werte von z asymptotisch 
gegen eine Funktion mit der Periode z geht. Dann folgt tatsachlich aus cth1)3 (2) = cthy)3 (z) die 
Beziehung 2 ~ 7 + 3. ; 

Ks ist 


1 3 
dis 1/3(2) Eres Ji/3(2) rare os Nj;3(2) ie ; 
wo Nj/3(z) die Neumannsche Funktion bedeutet, also folgt 


1 = 4 /3(2) ss 1 ae 3 N4;3() 
ctby /3(z) Ji/3(2) 2 2 J173(2) ; 


Es gilt asymptotisch 
nN BE 5x 2 eae 8 
1/3(2) ~~ ae sin = 2 5 J1/3() ~ ee ——— oe 


[etb/3(z)]-! ~ 3 te IB eg (: 5 


und somit 


Das ist aber eine Funktion mit der Periode z, so daB die oben angestellten Uberlegungen richtig 
sind. 
Ks 1aBt sich daher mit 


2 2 
y~aty Hatta ee 


4a P+ T?\3/2 a 
4 a? q 
aus (53) die Naherungsformel ableiten: 


12aqa? 2/3 _)4aP+ T? 
Tirit, Pind {la aP+ T?)3/2 1 7 4aq ° (54) 


Die Giite dieser Naherungsformel kann sofort an den bereits behandelten Zahlenbeispielen 
nachgepriift werden. Man erhalt mit ihr fiir 


: 3 2/3 
Beispiel 1: ee (ec a8 1 oon i 0,25 =12,8m, 

ee 12 z 60 - 101° 2/3 4-108 + 108 
Beispiel PAS Levit. Lod (q-t08--r055 aa 1) ae i 4-102. 60 = 22,9 m. 


10. Die Lésung der inhomogenen Differentialgleichungen des Stabes. Die allgemeine Lésung 
der inhomogenen Gleichungen (24) wird sich aus der allgemeinen Lésung (36) des homogenen Sy- 
stems und aus zwei partikularen Lésungen des inhomogenen Systems zusammensetzen. 

Die Ermittlung der partikularen Lésungen soll mittels der Methode der Variation der Konstan- 
ten erfolgen, wobei auch jetzt wieder die beiden reellen Lésungen zu einer komplexen zusammen- 


gefaBt werden sollen. 
Zuerst wird das inhomogene Differentialgleichungssystem (24) komplex geschrieben, und man 
erhalt mit den schon bekannten Abkiirzungen V und u 


d* . d —its 4 Te a —its 
| “t2iug—(w+V)o=Ate 7 +i B* e—its —Ce-it, (55) 
-  Setzt man (34) in (55) ein, so folgt mit 7 — u = w 
T—Vr=Ce—is, (56) 


: Es seien r, und r, die beiden linear unabhingigen Lésungen und r = Yr, + B rp die allgemeine 
‘Lésung von (35). Dann kann fiir die partikulare Lésung mit den zu variierenden ,,Konstanten” 
Q(s) und B(s) der Ansatz gemacht werden 


r = A(s) r, + B(s) ry. 
Das fiihrt fiir die Funktionen 2(s) und %6(s) auf die Bedingungsgleichungen 


Wr, +Br,=0, Hr, + Br, =—Ce-ivs. - +57) 
. 


; 


oe ae Cry ‘da LO) 
Ber Core TN 
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Nennt man die Determinante 


die die Determinante des Gleichungssystems (57) und auch die Wronskische Determinante des 
Fundamentalsystems der homogenen Basisgleichung ist, kurz W, so liest man fiir 2{ und § aus 


(57) ab 


0 rol ry 0 
& Ce—ivs Ff ity ee r Ce—iws _& —iws 
Ue U (oO = u W Sp "1? U e 


Im vorliegenden Fall ist iibrigens W, als Wronskische Determinante der Gleichung (35), stets 
eine Konstante, da das Gliedr in der Differentialgleichung nicht vorkommt. Deshalb erhalt man 


durch Integration fiir 2{(s) und ¥(s) 


AW(s) = = Fo Gate! a5, 8 (s) = {7 e, tees (58) 
6 6 


Geht man mit (58) in die Gleichung r = 2{{s) r, + B(s) rz ein, so ergibt sich als partikulare Lé- 
sung der inhomogenen Basisgleichung 


rab lin freinds + ry fry ote de (59) 
Mit (59) findet man 
0 = Seniws E rel tg ds ral zr, 6 ** ds|, 
0 0 


und die allgemeine Lésung der inhomogenen Gleichung (55) lautet dann 


: ity s : s : : 
@ =e +ie, ={Ar,(s) + Br,(s) + y us r,(s) ) r,(s) e—*”* ds +- r,(s) i ras) € 8 as| e_ ***_ (60) 
1l. Die Eigenwerte fiir Lagerfall 1 und q(s) = 0 bzw. q(s) =q = const. Fir den Lagerfall 1 
von Abb. 1 gelten dic Randbedingungen (26) und (27), die beide in komplexer Schreibweise anzu- | 


wenden sind: 
1 


o=e+tie,=—0 fir sgleichOunds =1, Slee +ie,) et* ds=—0. 
6 


a) q(s) =0. Hierfiir sind die Lésungen der Basisgleichung wieder 


T, =SIN US, Te = COSU)S, 


und es folgt aus (60) nach Ausrechnung der Integrale und von W 
0 =e +e, = Asin v)9s e—*** + B cos vgs e—*4¥* +C e— its, (61) I 


Die Eigenwerte erhalt man mittels der Randbedingungen. Die Bedingung 9 = 0 fiir s = 0 

fuhrt zu ; 
; 8 + ty aes 0 5 c SS Ds) ry 
womit aus der allgemeinen Lésung (61) 
@ =e +e, = Asin vgs e—'** + B(cos vgs e—*¥* — e— its) 

wird. 

Dann ist aber 

o e'** = W sin v9 s e§”* + G (cos vy s e§”* — 1), 


und die restlichen Randbedingungen ergeben die Knickdeterminantengleichung 


sin Vy 1 e— #4! cos vat e—iv! — ¢g—itl 
H l ai ; 
f sin v9 s e'”* ds f (cos v9 s e&”* —1) ds =0. (62) | 
0 0 
Es gilt noch 
v? — w? — 2 und — = Sage 
7 P A ape Tra alia: 


was im folgenden verwendet werden soll. 
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- Die weitere Ausrechnung bringt zunachst 


i 
. Sint a . . . 

Jesin tose weds = [e”! (i w sin v9 1— v4 cos vy 1) + v9], 

J : ae. : 

1B vos e¥s ds == [e'”! (i w cos vg 1 + v9 sin v9 1) —iw], 

und damit erhalt man aus (62) als Eigenwertgleichung 
; SES © pee i: : ; 
sin vy l e—** p le” (10.608 091 + v9 sin v1) —i we] —1| 
= (cos vy 1 e—*¥! — e— ith) iP [e’”! (i w sin vy 1 — vy cos v4 1) + vy]}. 


Multipliziert man aus und faBt weiter zusammen, so bleibt eine sich als reell erweisende Bestim- 
mungsgleichung fiir die Eigenwerte, wie sie bereits bei M. Beck! vorkommt: 


ag Hsin tl + 2 v9 (608 %1—c08 9 = AUN 


b) q(s) =q =konst. Setzt man zur Abkiirzung 


8 s 


— ra) r,(s) e—*”§ ds + r2(s) J r,(s) e'@* ds = p(s) , 
so schreibt sich (60) 


C 
= jx ry(s) + B ra(s) + wPO)| gi? 
und die Randbedingungen fihren zu den die Bestimmungsgleichungen 
9 r(0) + B r3(0) +7 p(0) =0, 
Ur) + Brel) + 4p) =0, 
I 1 l 
OW | r,(s) e's ds + B | r,(s) e'”' ds + — | p(s) e'”s ds = 0. 
J 4 w 


Hieraus folgt die Knickdeterminantengleichung 
(0) r2(0) p(9) 
r(0) rl) p(t) at 
[rls e Vids [ras e”s ds Sp(s) CGE 


Entwickelt man diese Determinante und beachtet, da8 p(0) = 0 ist, so erhalt man die Eigenwert- 


gleichung 


[r,(0) ro(1) — 72(0) nl p(s) e”* ds —r,(0) p(Z) f ro(s) e'”* ds + 15(0) p() f(s)” * ds = 0. (62) 


0 0 


Zur Durchfiihrung der weiteren Rechnung werden die Abkiirzungen 


s . s 


f r,(s) cos ws ds = SP? , f7,(s) sin ws ds = S??, 
: 0 
} (63) 
if-7As) cose ds = Soo, J r,(s) sin ws ds = S¥o. 
0 i) 


- eingefiihrt. Diese Abkiirzungen der Integrale seien als Funktionen von s aufgefaSt. Will man fiir 
eine dieser Funktionen einen Funktionswert fiir eine bestimmte Ordinate, zum Beispiel s = 1, 
 angeben, so sei dies durch 


l 
J r,(s) cos ws ds = S}? (64) 


angedeutet. — 


1 M. Beck, Ing.-Arch. 23 (1955) S. 246. 
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Es ist unter Verwendung von (63) 
p(s) =—14(s) Seo + ra(s) Ste + # (ra(s) St’ — 72(s) Sze) 
und 


1 1 
: . 0 5,0 . 
f p(s) &”* ds = f [S?) r, cos ws + S*, r. sin ws —S;, r, cos ws —S,, r, sin w s] ds 
i 


+ if [S*° r, sin ws +S; r, cos ws — Si. r, sin ws —S¥, 7, cos ws] ds. (65) 
Ferner ist 
—1,(0) p() r9(s) e"”* ae (Seo + 4 Se )'7,(0) [(Sh2 —t She) ra(l) — (S32 — 7 S32) (I. 
und das gibt | 
— 1,(0) p(2) i rs(s) e”* ds = 1,(0) r4(1) [(S$2)? + (SE2)7] + 14(0) ra) [— (Ste Sz? + Sr? See | 
+ (Sts She — Sit Si2)]. (66) | 
Ebenso erhalt man 


l 


IES a ag a (Sho + 4 St?) r2(0) [ra(l) (St? + ¢ St2) — 14D) (S32 — 7 S32), 


was bei weiterer Ausrechnung auf 

1 

20) p(L) f r(s) e”* ds = r,(2) r9(0) [— (See St? + S32 Ste) + i (S32 Sve — S32 Sty) : 
6 

+ 79(0) ra(Z) (S32)? + (Sr9)?] (67) 
fiihrt. | 
Setzt man die Ausdriicke (65) bis (67) in die Eigenwertgleichung (62) ein, so erhalt man eine Be- |} 
dingungsgleichung, bei der man Real- und Imaginarteil getrennt schreiben kann, und beide Teile || 
miissen fiir sich zu Null werden. 
Auf Grund einer leichten Rechnung 1a8t sich zeigen, daB der Imaginarteil identisch Null ist. | 
Als Kigenwertgleichung fiir den Lagerfall I bleibt also der Realteil von (62) itibrig. Verwendet man 

die weiteren Abkiirzungen 


l 
L(l) = f [Sf 1, cos ws + Sv? r, sin ws — $3? r, cos ws — S;") r, sin ws] ds, (63) | 
0 | 
Tn, (I) = Sho See Sto She, - Ey(D) = (Sh)? 4 (Ses -E.() = See)" + S52)". | 
so 1aBt.sich der Realteil schreiben als ; 
r4(0) ra(Z) [L(2) — Ly2(2)] — r2(0) ry(Z) [L() + Li,2 ()] +740) 4D) Lp(2) + 79(0) ro(l) Ly) =0, (69) | 


woraus man die Verwandtschaft zur Kigenwertgleichung des Lagerfalles 2 erkennt, die einfacher, |} 
namlich | 
r,(0) ra(Z) — r2(0) r,(1) = 0, (70) 
lautet [vergleiche hierzu (50) ]. | 
Fir den Belastungsfall q(s) = q = konst. gibt die Basisgleichung gemaB (47) die Lésungen 
2 2 
1 = 3? Sys & 3 ) te ete a 5 3 ) 


Diese sind in (69) einzusetzen, um fiir diesen Fall die Eigenwerte zu errechnen. 

Zur Erléuterung soll ein Zahlenbeispiel dienen. Gegeben sei T = 10'kpm, a = 105 kp m?, 
q = 10’ kp/m. Die Rechnung wird mit einer geschatzten Stablinge begonnen. Die Eigenwert- 
gleichung (69) wird dann nicht exakt erfiillt sein, sondern man erhalt fiir verschiedene angenom- || 
mene Werte von / eine Fehlerfunktion, deren erste Nullstelle man aufsuchen mu. Man findet diese 
durch Eingabeln. Zu dieser ersten Nullstelle gehért als Abszisse die gesuchte Stablange [,,;:. 

Einen ersten Hinweis auf die zu wahlende Stablange bekommt man dadurch, daB man die Stab- |} 
lange fiir den Lagerfall 2 bestimmt, die sich leicht finden 1aBt. Mit den gegebenen Werten der Auf- || 


gabe ist 
2 gy 2sf/ OL \ /105\2/3)3/2 
{= ore (a {0} fc | = 8,31. 
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Da fiir so groBe Werte von z ohne weiteres die Naherungsformel (54) verwendet werden kann, 
-erhalt man durch diese 


ie ~) 


‘(12 2 103. 1010 ie 120° 
1015 ag lat 


lw [(1,377)22 —1] 25 = 0,24-25 —6m. 


Nun beachte man, daf fiir den Lagerfall 1, wegen der gréBeren Zahl der Bindungen des Stabes 
gegeniiber Lagerfall 2, die Knickfestigkeit gréBer wird. Dann kann man annehmen, daB auch die 
Knicklange des Stabes im Lagerfall J gréfer sein wird als die Knicklange im Lagerfall 2. Es soll 
daher beim Lagerfall I zunachst mit J = 10 m gerechnet werden. Nach (43) ist 


104 105 
10° | \2.105 


Damit erhalt man nach (45) 
4 =a B— 7/3 — BS 5 — 7,54 —0,215s. 


108 


2 
ee ) = 035, = 755 = 0,01. 


Ferner ist 
ct T 105 


as Fa. 2 2 k05 


— ae 


Die weitere Rechnung erfolgt in Tabellenform. Die zum Teil vorgeschriebenen Integrationen 
werden an Hand der Tabellenwerte nach einem der gebrauchlichen Verfahren der praktischen 
Analysis ausgefiihrt. Man sehe Tabellen 1 und 2. 


Tabelle 1 

Cu 0 1 2 ] 4 2 6 a 8 9 10 

= 7,040 | 7,325 | 7,110} 6,895) 6,680} 6,465} 6,250) 6,035) 5,820/ 5,605} 5,390 
z= 2/339? = | 13,87 | 13,27 | 12,62 | 12,14 | 11,47 | 10,94 | 10,41 | 9,94 9,34 8,87 | 8,34 

J1)3(2) = 0,214 0,179 | 0,068 | —0,039|— 0,175] —0,236| —0,234|—0,176 |—0,044| 0,079) 0,204 
J_ 13(2) = 0,110 | 0,198 | 0,220 | 0,177) 0,048 —0,077| —0,186|—0,245 —0,245 |—0,182 | —0,059 

ON 0,0 | 0,5 1,0 | IES 2,0 TAs 3,0 ac0 4,0 4,5 5,0 
cos ws = 1,0 0,878 | 0,540, 0,071 —0,416 —0,801] —0,990 —0,936 |--0,654 |—0,211| 0,284 
sin ws = 0,0 0,479 | 0,841) 0,997} 0,909} 0,598) 0,141'—0,351 |—0;757 |—0,978 | —0,959 
31/2 Jy/(2) cosws | 0,588 0,426) 0,098 —0,007| 0,188) 0,480) 0,579 0,405; 0,069|—0,040, 0,134 
a? Jyjs(2) sin ws | 0,0 0,232 | 0,153 | —0,102| —0,410| —0,358) —0,082} 0,152) 0,080 |—0,183 | —0,454 
3/2 J_ 1/3(z) cos w s| 0,303 | 0,471 | 0,317) 0,033) —0,052! 0,157; 0,460) 0,564} 0,386) 0,091} —0,039 
3¥2 J_13(z)sinws| 0,0 | 0,257| 0,494) 0,464) 0,113; —0,117;— 0,066] 0,212) 0,447) 0,422; 0,131 


Tabelle 2 
s= 0 1 2 3 4 5 6 a 8 9 10 
ae ee en ae eee ee 1,769} 2,261) 2,498, 2,512} 2,559 
gee = 0,0 | 0,116 | 0,309 | 0,334 0,078| —0,306| —0,526| —0,491|—0,375 —0,426 | —0,744 
' Spe = 0,0 | 0,387 | 0,781 | 0,956 0,947) 1,000) 1,308} 1,820; 2,295) 2,534| 2,560 
oe — 0,0 | 0,128 | 0,504 | 0,983 1,271) 1,269) 1,177) 1,250) 1,579) 2,013; 2,290 


SS2° 7r_coswsds =| 0,0 | 0,119| 0,360] 0,495 | 0,485, 0,559 1,004) 2,109) 3,228) 3,824) 3,889 
0 


SS2o rgsinwsds=| 0,0 | 0,015| 0,106 | 0,260 | 0,342) 0,365 0,400) 0,366, 0,230 0,056 —0,082 
0 


c 


SS2°r,coswsds=| 0,0 | 0,082| 0,203| 0,238 | 0,323/ 0,652, 1,270} 2,017) 2,465, 2,493) | 2,614 
J 


fS2°r,sinwsds=| 0,0 | 0,015 | 0,068 | 0,057 |—0,254 —0,741| 1,016] —0,969 —0,811 —0,932 /—1,636 
0 


2s 
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Mit den Werten der Tabellen errechnet man 


) 

r,(1) = 3)? Jays (8,34) = 75,39 - 0,204 = + 0,473 , 
) = 342 J_ 1 (13,87) = //7,54 - 0,110 = + 0,308 , 
) = 34? J_1)s (8,34) = //5,39 (— 0,059) = — 0,137 


und nach (68) 

L(l) = 3,889 — 0,082 — 2,614 + 1,636 = + 2,839, 
= 2,559 - 2,560 — 0,744 - 2,290 = + 4,847 , 
= 2,5592 + 0,744? = 6,548 + 0,554 = + 7,102 , 
L,(l) = 2,560? + 2,290? = 6,554 + 5,244 = + 11,798. 


Die Eigenwertgleichung (69) gibt 
— 0,588 - 0,137 (2,839 — 4,847) — 0,303 - 0,473 (2,839 + 4,847) 
+ 0,588 - 0,473 «11,798 — 0,303 - 0,137 - 7,102 = + 2,046. 


Da sie nicht erfiilJt ist, sondern den Fehler f = + 2,046 aufweist, wird die Rechnung mit einem ' 
neuen Wert von /, namlich mit 1 = 8 m, wiederholt. Es andert sich dadurch der Wert von a: 


_ 8-108 
105 


coats) 8 0,25 = 0,32 
+ (a-a08) = 0,08 + 0,25 = 0,32 , 


und es ist somit 3 = 6,894 — 0,215 s, wogegen ( und w unverdndert bleiben. 
Setzt man wieder ganz entsprechend wie vorher die Rechnung fiir / = 8 m in Tabellenform an, 
so erhalt man 


r,(0) =—0,103, — 1,(1) = 0,609, 7(0) = 0,466, — r,(l) = 0,119, 


Sie == 2,028 , See = — 1,520, oie == 1.620 5 Sho = 1,330, 
1 1 
J Ste 7, cos ws ds = 3,463 , fs S}. 7, sin ws ds = — 0,340, 
0 0 
I 1 
J Sie 7, COS WS ds = 3,105 ? 3 S33 Ty sin ws ds = — 2,709 9 
0 0 


womit nach (68) 


L(Y) = + 2,727,  Ina(l) = +1,278, L() =+6,423, L,(l) = + 4,416 
und nach (69) 
— 0,103 - 0,119 (2,727 — 1,278) — 0,466 - 0,609 (2,727 + 1,278) 
— 0,103 - 0,609 - 4,416 + 0,466 - 0,119 . 6,423 — — 1,076 
folgt. 
Die gesuchte Knicklange des Stabesfindet man aus den beiden Fehlern der Eigenwertgleichung zu 


(10 — 8) 2,046 


Lowe TOO 
writ, 10 2,046 + 1,076 


= 10 — 1,311 = 8,7m. 


(Eingegangen am 17, Februar 1960.) 


Anschrift des Verfassers: Dipl.-Math. Horst Leipholz, Stuttgart-Untertiirkheim, Fiechtnerstr. 51. 
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Die Kreiszylinderschale unter konzentrierten Belastungen* 
Von K. W. Bieger 


: 1. Einleitung. Die Berechnung von Kreiszylinderschalen fiir beliebige Oberflachenlasten wird 
in der Literatur im wesentlichen nach zwei Methoden behandelt. 

Die erstere nach der genaueren Theorie von Fliigge! erfordert die Lésung der drei gekoppelten 
partiellen Differentialgleichungen fiir die Verschiebungsfunktionen u, v und w. Zur Lésung werden 
hierbei meist Doppelreihen? herangezogen, die vor allem bei konzentrierten Belastungen zu erheb- 
lichem Rechenaufwand fiihren. Koepcke® lést das Problem der konzentrierten Belastungen durch 
Uberlagerung zweier Lastfalle am unendlich langen Rohr, wobei die Kontinuitat durch die Be- 
nutzung des Randstérungsproblems erzwungen wird. 

Die zweite Méglichkeit ist die Vernachlassigung mehr oder weniger vieler Glieder in den Zu- 
sammenhangen zwischen Schnittlasten und Verschiebungen‘, so da man die drei gekoppelten 
Differentialgleichungen auf eine einzige achter Ordnung fiir die Verschiebungsfunktion w zuriick- 
fiihren kann. Es hat sich gezeigt®, daB diese Vernachlassigungen vor allem bei konzentrierten Be- 
lastungen zum Teil unzulassige Ungenauigkeiten mit sich bringen. 

In der vorliegenden Arbeit soll dargetan werden, da es méglich ist, durch Einfiihrung einer 
Lésungsfunktion, ahnlich der Spannungsfunktion bei der Berechnung der Scheiben, die drei Fliigge- 
schen Differentialgleichungen zu einer einzigen achter Ordnung fiir die Lésungsfunktion zusammen- 
zufassen. Dadurch ist die Berechnung der Kreiszylinderschalen vor allem fiir konzentrierte Be- 
lastungen wesentlich vereinfacht. Die Zusammenhange zwischen den drei Verschiebungen — und 
damit zwischen den Schnittkraften — und der Lésungsfunktion sind im dritten, vierten und 
fiinften Abschnitt angegeben. 

Fiirden Sonderfall der radialen Belastung und mit der Annahme, daf der Schalenparameter k 
gegentiber 1 vernachlassigt werden kann, was im praktischen Gebrauch fast immer médglich ist, 
sind diese Zusammenhange auch bei Wlassow® angegeben. Da jedoch Wlassow nicht die von Fliigge 
aufgestellten Differentialgleichungen als Ausgangsbasis benutzt, ergeben sich noch einige Unter- 
schiede, so vor allem in der Differentialgleichung der Lésungsfunktion. 

An Hand des Beispiels einer radialgerichteten, im Gleichgewicht stehenden Einzellastgruppe am 
unendlich langen geschlossenen Rohr wird im sechsten Abschnitt gezeigt, daB dieser Lésungsweg 
eine explizite Darstellung der Konstanten der Greenschen Funktion der Kreiszylinderschale er- 
méglicht und damit eine vereinfachte Berechnung von EinfluSflachen gestattet. 


2. Grundlagen. Es gelten die in der Schalentheorie iiblichen und schon von Fliigge angegebenen 
Voraussetzungen. Die Schnittlasten, Verschiebungen, Belastungen und Koordinaten sind nach 
Abb. 1 positiv angenommen. 

Werden fiir die Differentialquotienten zur Abkiirzung die Bezeichnungen 


of ‘ Of mir? 
es und (Frese 


eingefihrt, so ergeben sich mit den von Fliigge aufgestellten Zusammenhangen zwischen den Schnitt- 
lasten und Verschiebungen die drei gekoppelten partiellen Differentialgleichungen der isotropen 
Kreiszylinderschale 

Wey 


pa tix: 2 
wu" + ; WE as ai : Wee aay? 1 5 sd Se (la) 


* Auszug aus der Dissertation des Verfassers, EinfluBflachen der Kreiszylinderschalen, T. U. Berlin 1959; 
Berichter: Prof. Dr.-Ing. W. Koepcke und Prof. Dr.-Ing. J. Szabo. ; 

1 W. Fliigge, Statik und Dynamik der Schalen, 2. Aufl., S. 150. Berlin-Géttingen-Heidelberg 1957. : 

2 7, B. Fr. Dischinger, Beton und Eisen 1935, S. 260, oder C. B. Biezeno u. R. Grammel, Technische Dynamik, 
2. Aufl., Bd. 1, S. 524, Berlin-Géttingen-Heidelberg 1953. ; ; ew) 

3 W. Koepcke, Die Berechnung von Kreiszylinderschalen unter Flachen-, Linien- und Einzellasten, Habili- 
tationsschrift der T. H. Berlin 1949. — W. Koepcke, Uber kreiszylindrische Schalentrager unter Flachen-, 
Linien- und Einzellasten. Unveréffentlichte Arbeit zum 60. Geburtstag von Prof. Dr.-Ing. Dischinger, Berlin 1947. 

4 vgl. z. B. K. Girkmann, Flachentragwerke, 4. Aufl., S. 475 und S. 490, Wien 1956. 

5 vgl. auch die Dissertation des Verfassers. : ; ; 

6 W.S. Wlassow, Allgemeine Schalentheorie und ihre Anwendung in der Technik, S. 210, Berlin 1958. 


f 
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ane Sel l=. 9 eh 3(1—») 3—? r) any lb 
hs is mar aes ee a ea eas D’ (1b) 
st SL yt a ; i : a? 
yut+tv+tw- k(- ee ee = Po tw™ 4207 4+ w 4 Rut tw) = — ZF (1c) 
Hierbei ist y die Querkontraktionszahl und 
{2 
is Da 12a 
der sogenannte Schalenparameter, worin 
pee Ee 
12 (1— 7) 
die Biegesteifigkeit und 
sili ts 
~ 1—? 


die Dehnsteifigkeit der Schale bezeichnen. 
Der Schalenparameter k liegt fiir die im Bauwesen auftretenden Falle zwischen 10-4 und 10-§ 


und kann daher fiir praktische Berechnungen gegeniber | vernachlassigt werden!, 


Abb. 1, Schalenelement. a) Dehnungskrafte, Querkraéfte und Belastungen; b) Momente und Verschiebungen. 


3. Lésungsansatz fiir Lasten in Richtung der Erzeugenden. Die drei gekoppelten Differential- 
gleichungen (1) lassen sich mit Hilfe der Losungsfunktion F(x, y) dann am iibersichtlichsten lésen, 
wenn die Oberflachenlasten in den drei Hauptrichtungen fiir sich untersucht werden. Da sich 
die fiir jede dieser Belastungen aufgestellten Differentialgleichungen fiir F(x, q) allein in den Ab- 
solutgliedern andern, ist hierdurch nur eine geringe Mehrarbeit verbunden. Bei gleichem funk- 
tionalem Aufbau der drei Belastungsarten und der Lésungsfunktion unterscheiden sich die fiir die 
jeweilig betrachtete Lastrichtung ermittelten Funktionen F(«,q) nur durch eine Konstante. Die ent- 
sprechend den folgenden Angaben errechneten Verschiebungen der einzelnen Belastungszustande sind 
dann nur zu summieren. Zum anderen sind in der Praxis meist die Lasten in Richtung der Erzeu- 
genden Null und die Lasten in Ringrichtung klein, so daB sie zum Teil vernachlassigt werden kénnen. 

Setzt man fiir Lasten in Richtung der Erzeugenden die Verschiebungen in Abhangigkeit von 
der noch zu bestimmenden Lésungsfunktion F(x, pm) in folgender Form an 


2 at tn: ae 5 — Sie i ; 
Mae = Coy EE —.2 k FF: ( _ ! . a 


3 — 6v— 7? 


— BQ») +30») k] KEY —|2—s) + 2 


Kk Be Ws 


1— == 

— ar +4k+3i)F’—*> 2 +3k) kere, 
Die 9, of Lay 2 = ein ay eee <n . 2) 
Mie GaaK he tet to ao (ee ee ee L 


1—y. 1+» > , 2+3y—p sate 1 — 
+ a k)F oe a +( pg ear, 


4 
Wan = K Lt RP Fed Ba 


1— . 
Sip Fe 4 (L$ 3k)» F138) k Fer} 


1 Ks darf jedoch nicht allgemein k = 0 gesetzt werden. 
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dann ergibt sich nach Einsetzen der Verschiebungen (2) in die Gleichung (la) die Differential- 
gleichung des Problems 


UU es Pal : 
(1 +2k—3h) F +(4 +k at eee + [6 + (6—3»)k—o2 ke] Fee 


! (4 ! 7 ; Fer = 8) PO + (Lb) + Oy 4 Ook) PU 4 [6-432 —y 08) BM 


1—y7 
k 


+ [8—2y +(7—5r)k +3 (1—») B] FU 42 (1 +k) FF +| + (4—39) +3k/ RF" 


+(4—274 EO ie Sy BlFe + +h) FP = 4 Xq): (3) 


Die Gleichungen (1b) und (1c) werden durch (2) identisch erfiillt, wenn man beriicksichtigt, daB 
wegen der Aufspaltung der Belastung fiir den betrachteten Fall Y(x, 9) = Z(x, py) = 0 zu setzen ist. 


4. Lésungsansatz fiir Lasten in Ringrichtung. Fiir Oberflachenlasten in Ringrichtung sind die 
Verschiebungen folgendermafen anzusetzen: 


2 at 1lt+y eine 1+3y yee 3—4y+7 yrs 
Me Tok {t hr 4 +t?) p +049) PT ake 


4 ere =| O eeetmeena =r aeerea a Reae 


2 2 2 
Sepa! Sey es ZL + ARE ——»)(1 +h)kF*—|(2—») coemer ee Fu x 
= +k)? F?—(2—y+vr)k PF" [yaaa ae 
— [(l— vr?) +k] Fl—2yk PM —( WRF, 
M0) = {Cl + k) tele chee eb (a?) B” ~2kEr J 
In diesem Fall ergibt sich die Differentialgleichung der Kreiszylinderschale aus der Gleichung (1b) 
(l+2k—32) rer gs er + [6 + (6—3»)k—v 2] FO" 


Saar a ae +k) F**** + (2y + 6k) F""" + [643(2—y +7) k] F'" 


1—yr 


(8 —29 +(7—5e)k +30 —vy) BLP" +2(1 +4) F4 - + (4—3 9%) + 3h] FM 


+ (4-294 aay deere) Me (+B) PES Vedeya (5) 


wahrend die Gleichungen (la) und (1b) wiederum identisch erfillt sind, wenn voraussetzungs- 
gemaB X(x, p) = Z(x, p) = 0 angenommen wird. 


5. Lésungsansatz fiir Lasten in radialer Richtung. Fiir Radialbelastungen liefern die Ver- 


schiebungen 


Mag) = REP FTEs + (1+ 3 (4 RE, 


V(x, 9) = (a+ F: —(2 tA Gta RFY +e +9) FV —2k Fr, (6) 


2 


W(x, ¢) =_{-4 + k) em ++>*6@ +32) me (1 a 3K Fr) 


60 K. W. Bieger: Die Kreiszylinderschale unter konzentrierten Belastungen Ingenieur-Archiv 


aus der Gleichung (1c) die Differentialgleichung achter Ordnung fiir die Lésungsfunktion F(x, ¢) 
(1 ae Oy a k?) pugs (4 See = 37), = ae) pun =o [6 Sis (6 — 3) k—y? k?] FU a 


4 (4 ! — Eee +k) Fe! 1 (2y + 6yk)F""" + [6 +3(2—v +) k] F"™ 
7 


+ 2y + (59) k £30 — 9) Pl ee ee) +| "+ (4399) +n) Fr 


i (4 ye Fe +h) FY = + 2Z(x,9). (7) 


Aus den Gleichungen (2), (4) und (6) laBt sich ablesen, daB — wie fiir jeden elastischen Kérper — 
so auch fiir die Kreiszylinderschale der Satz von der Gegenseitigkeit der Verschiebungen und 
Lasten (Maxwell- Betti) gilt, sofern man noch beachtet, da8 es bei dem unendlich langen geschlos- 
senen Rohr hinsichtlich des Verschiebungszustandes des 
Punktes (x,g) nur auf die Differenz der Ordinaten dieses 
Punktes und des Lastangriffspunktes (also nur auf (x — xp) 
und (y—gp)) ankommt?. Denn mit den Bemerkungen zu 
Abschnitt 3 folgt: Fir eine Last X = 1 in Richtung der Er- 
zeugenden wird die Verschiebung v(x,g) in Ringrichtung 
gleich der Verschiebung in x-Richtung u(x, p) infolge einer 
Last Y = 1. Es gelten somit folgende Identitaten: 


v(x, P)x_ 1 = U(x, P)ya1? 
w(x, P)x_ 1 = U(*, P)z_1> 


Abb. 2. Radiale Einzellastgruppe. w(x, P)y_ = v(x, P)z ane 


6. Beispiel fiir eine radiale Einzellastgruppe. Um eine von allen Randbedingungen unabhiangige 
Losung fiir radiale Lasten auf Kreiszylinderschalen zu erhalten, wird eine im Gleichgewicht stehende 
Einzellastgruppe nach Abb. 2 auf das unendlich lange geschlossene Rohr aufgebracht. Aus diesem 
1aBt sich, je nachder Art des zu untersuchenden Tragwerkes, das vorhandene Gebiet herausschneiden 
und mit Hilfe der bekannten Randstérungstheorien die notwendigen Randbedingungen hefriedigen. 

Um die interessierende Kinzellastgruppe als stetige Funktion der Koordinaten x und  darzu- 
stellen, wird diese vorerst in y-Richtung in eine Fouriersche Reihe mit der Periode 2 und an- 
schlieBend in x-Richtung in ein Fouriersches Integral entwickelt: 


/ ~ oe} 
poe Re 
Ze9) = aaa | cosa Fda eae eae ng [cosozas), (8) 
a=0 a=0 


Wird als Beispiel der zweite Summand in der geschweiften Klammer betrachtet — die Behandlung 
des ersten Summanden, der dem in Ringrichtung konstanten Anteil der Last entspricht, unter- 
bleibt hier — und wird fiir die Lésungsfunktion der von den Veranderlichen x und y in gleicher 
Weise abhangige Ansatz _ 


2P © ~ 
Fng)=22, 


R ce 
cosn@ | Cy, cos « —da (9) 
=2,4,6... ae 


a=0 


gewahlt, dann folgt nach Ausfiihrung der Differentiationen und Einsetzen in Gleichung (7), wenn 
man noch zur Vereinfachung k gegeniiber 1 vernachlassigt, 


2P r 4 
Bai a 008 np | Cam CHa, n) cosa da = 4, 5 cos ng | cos a do. (10) 
Hijesuniee a eae 
CH (a, n)=a8 + 2 (2 n? —y) a6 + (6 nt — 6 n? + ae 
+2 [2 n® —(4—») nt + (2 —») nt] o# + (n?§ —2 m8 + nt) (11) 


1 Herrn Priy.-Doz. Dr.-Ing. R. Trostel bin ich fiir einen Hinweis bei diesen Zusammenhangen dankbar. 
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oe mit der charakteristischen Gleichung der isotropen Kreiszylinderschale nach Biezeno- 
rammel}, 


Die Gleichung (10) ist nur dann fiir jedes x/a und y erfillt, wenn 


1 


Can = C H(a,n) 


ist. Setzt man dieses Ergebnis in den Ansatz (9) ein und verfahrt in gleicher Weise mit dem ersten 
Ausdruck der Klammer in (8), dann ist mit der Lésungsfunktion 


A oP(1 par tee r cos a — 
(0) = as 2 CHimie eo, CFT ioe & ) 


=0 n= 2,4,6.. 0 


unter Benutzung der Gleichungen (6) der gesamte Spannungszustand des kreiszylindrischen Rohres 
infolge der betrachteten Einzellastgruppe bestimmt. 
Die Integrale in (12) lassen sich mittels Partialbruchzerlegung geschlossen lésen. Werden die 
_ Wurzeln? der charakteristischen Gleichung (11) mit 


ae = +m tip, und 45,6 = +E Me = ty 

; 7,8 

bezeichnet und zur Vereinfachung 

Fog og (rep my + 2 oe my + 5 2g + eg — 10 oe5 3 — 6 xh xe — 2 it 2 + 6 x8 2 + 2 x2 2) , 


p= + 4x8 Me (og + mi + 2 at + 2g + 5 og — 10 29 wg — 2 ve ne — 6 2 2 + 2 oh? + 2 x2 Ud) | 
ys — 4m wi (Bog + Mi — LO mi ay + eg + oa + 2 20g mg — 6 xy wg — 2 at uh + Arg ut + 6 re Us) » 
OS + A et py (244 + 5 ey — 10 vet at + ed + ng + 2 3 uy — 2 oof meg — 6 3 ug + 6 x3 ut + 22? 2), 


eingefiihrt, dann lat sich nach einiger Zwischenrechnung die Lésungsfunktion in folgender Form 
angeben: 


n= 2,4, 6, 
ps 2y x pe ee x 
+- e fae cos Ho re ety to “)} 5 (13) 
Hierin ist — wieder mit der Naiherung k < 1 — fiir den in Ringrichtung konstanten Anteil der 
Belastung 
4 1 —y? 
=) ae 
zu setzen. : 


7. Auswertung.. Es 1a8t sich zeigen, daB fiir das unendlich lange geschlossene Rohr die Ein- 
fluBflachen der statischen Gréfen fir eine radiale Einzellastgruppe — bis auf das Vorzeichen bei 
einigen Schnittlasten und Verschiebungen — identisch mit den entsprechenden Zustandsflachen 
sind. ; 

Am Lehrstuhl fiir Stahlbetonbau der Technischen Universitat Berlin wurden auf Grund der im 
vorstehenden in groBen Ziigen angegebenen Theorie fiir die beiden extremen Schalenparameter 
k = 10-4 und k =10-* die Einflu®flachen der Verschiebungen und wesentlichen Schnittlasten 
numerisch berechnet? und in Form von Héhenschichtplanen dargestellt* (Abb. 3 und Abb. 4). 
Wegen des grofen Einflusses des Schalenparameters auf die OrdinatengréfBen und die Gestalt der 
EinfluBflachen ist beabsichtigt, in Kiirze auch den zwischenliegenden Schalenparameter k = 10~° 
auszuwerten. 


1 C. B. Biezeno u. R. Grammel, S. 528, FuBnote 2 von S. 57. 

2 Uber die Wurzeln der charakteristischen Gleichung vgl. F. Dischinger, Beton und Eisen 1935, S. 263 und 
J. Moe, Abh. Int. Ver. f. Briickenb. u. Hochbau 1953, S. 283. 

3 Hierbei wurde die Querkontraktionszahl zu » = 1/, angenommen. 

4 Diese sind in der Dissertation des Verfassers abgedruckt. 
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Es bereitet mit den nunmehr bekannten EinfluBflachen keine Schwierigkeiten, fiir beliebige 
radiale Oberflachenlasten die statischen GroBen an interessierenden Stellen zu ermitteln. Die Ver- 
schiebungen und Schnittlasten an den Randern des aus dem unendlich langen geschlossenen Rohr 
herausgeschnittenen interessierenden Gebietes werden nicht mit den Randbedingungen des gegebe- 
nen Schalendaches in Einklang stehen. Sie werden jedoch lings der Rander nach verhaltnismafig 


at kes 


#2 


Abb. 4. EinfluBfliche des Biegemomentes Tre L058 


My (100z-fach) k = 
glatten Kurven verlaufen, sofern nicht eine Einzellast in unmittelbarer Nahe cines dieser Rander 
steht. Es kénnen daher mit geniigender Genauigkeit Naherungsverfahren benutzt werden, die es 
gestatten, an Hand von Tabellen’ verhaltnismaBig schnell die geforderten Randbedingungen zu 
erfiillen. : : 

(Eingegangen am 7. Marz 1960.) 


Anschrift des Verfassers: Dr.-Ing. Klaus-Wolfgang Bieger, Lehrstuhl fiir Stahlbetonbau 
(Prof. Dr.-Ing. W. Koepcke), Technische Universitat Berlin, Berlin-Charlottenburg, Hardenbergstr. 34 


* D. Riidiger u. J. Urban, Kreiszylinderschalen, Leipzig 1955; A. Aas-Jacobsen, Die Berechnung der Kreis- 
zylinderschalen, Berlin 1958. Re eS ENS 
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Beanspruchung und Verformung rotierender Scheiben 
durch axiale Drehmomente 


Von Karl Karas 


1. Einleitung. Die Beanspruchung und Verformung rotierender Scheiben beliebigen Profils 
durch radiale Massenbeschleunigungen ist in weitem Mae durchforscht worden, wahrend analoge 
Untersuchungen hinsichtlich azimutaler Massenbeschleunigungen bisher in der Literatur nur in sehr 
beschrinktem Mafe zu finden sind’. Dabei sind die bzgln. Fragen hinsichtlich radialer und azi- 
mutaler Massenbeschleunigungen vom theoretischen Standpunkte einander véllig gleichberechtigt 
zugeordnet. Wenn trotzdem die ersteren gegeniiber den letzteren bisher eine so tiberragende Be- 
achtung gefunden haben, so mag der Grund dafiir wohl vor allem in dem Bediirfnis nach der Beherr- 
schung des stationdren Dauerbetriebszustandes gesucht werden, der natiirlich fiir alle Turbo- 
maschinen eine viel groéBere Bedeutung besitzt als der nur gelegentlich auftretende instationdre Be- 
triebszustand mit grofen azimutalen Massenbeschleunigungen. Dieser kann insbesondere durch 
Drehmomente bewirkt werden deren Achse mit der Drehachse der Scheibe zusammenfallt. 


2. Herleitung der Grundgleichungen. In Abb. 1 a—c sind die Scheibenabmessungen, ein axial 
witkendes Drehmoment M, mit der (konstant vorausgesetzten) Winkelbeschleunigung 2 = @ und 
die durch M, bewirkten Schnittkrafte am Scheibenelement eingezeichnet”, wobei y das Einheits- 
gewicht, 4 = y/g die raumliche Dichte und G den Schubmodul des Scheibenmaterials bedeuten. 


tryag+a(try)/ap 


Ave 


Wellertachse 


Abb. 1. a) Schnitt durch die Scheibe; b) Ansicht derselben mit durch schwarzes Feld gekennzeichnetem Scheibenelement; 
c) VergréBert gezeichnetes Scheibenelement mit den durch das axiale Drehmoment bewirkten Schnitt- und Tragheitskraften. 


Ist ferner v die azimutale Verschiebung eines Scheibenpunktes, so gilt fiir den Zentriwinkel 
der azimutalen Verschiebung: 


Ga. (1) 


Die Achsensymmetrie des Schubspannungszustandes bedingt die Unabhingigkeit von gm und 
damit die Entkoppelung von den Gleichungen des Dehnungsspannungszustandes. 


1 Mit der vorliegenden Arbeit beriihren sich die folgenden Untersuchungen von Herrn Grammel, denen sie 
auch ihre Anregung verdankt: R. Grammel, Z. angew. Math. Mech. 5 (1925), S. 193; Ing.-Arch. 6 (1935), S. 256. 
Man vgl. auch: C. B. Biezeno und R. Grammel, Technische Dynamik Bd. 2, S.3 und S. 29 ff., 2. Aufl. Berlin- 
Géttingen-Heidelberg 1953. 

2 Die Bezeichnungsweise ist dem in | erwahnten Werk von C. B. Biezeno und R. Grammel angeschlossen. 
Abb. 1 a bc entspricht den Abbn. 31 und 32 auf S. 32 des erwahnten Werkes. 
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Aus der Gleichgewichtshetrachtung eines Scheibenelementes in azimutaler Richtung (Abb. 1c) 
ergibt sich die erste der folgenden Gleichungen: 


1 d(ty) T 2 
a ae Bree +par=0, (2a) 
dv v 7A 
mp ae i 


deren zweite den azimutalen Verschiebungszustand in bekannter Weise beschreibt. Eine der drei 
Funktoinen y(r), t(r) oder v(r) kann gewahlt werden. Gewéhnlich ist die Profilkontur gegeben. 
Nur bei Wahl der Verschiebungsfunktion v(r) ist zur Ermittlung von t(r) und y(r) blo& eine Inte- 


gration erforderlich. 
Eliminiert man aus (2a, b) die Schubspannung 7, so folgt fiir v(r) die Differentialgleichung 


zweiter Ordnung 4 ; 
dv 1 dy 1 v v Ar _ 
ae + (5 eee cae G pr: (3) 


die spaterhin mehrfach zu Kontrollzwecken ausgeniitat wird. 

Eine etwa vom Schaufelkranz und den Schaufeln am AuBenrand der Scheibe herrthrende Mas- 
senbelegung m, je Langeneinheit des Randes wird wegen der geforderten Achsensymmetrie unab- 
hangigvon gy vorausgesetzt. Ist z. B. r, der Abstand des Schwerpunktes der Querschnittsflache f;, 
des Schaufelkranzes von der Wellenmitte und analog r, derjenige der Schaufelschwerpunkte, so er- 
gibt sich mit G, als Gesamtgewicht aller Schaufeln der Lange /, deren Querschnitt konstant angenom- 
men wird }, 


G, r?+ P/12 rp 
ong r8 + th. S- (4) 
Fir das axiale Massentragheitsmoment T der Scheibe mit der linearen Massenlbelegung m, am 
AuSenrand findet man 


T =22m,r3 +20 freydr. (5) 


Die Ubertragung des Drehmomentes M, von der Welle auf die Scheibe erfolgt an der Innenflache 
2 ry Yo derselben mittels der dort wirksamen Schubspannung 1(rp) = T). Mit T in (5) gilt daher: 


Gleichung (6) bietet im Zusammenhang mit (5) eine willkommene Kontrolle fir T und 1). Die 
GréBen t(r) und v(r) werden durch Integration der beiden linearen Differentialgleichungen erster 
Ordnung erhalten. Dabei wird die Integrationskonstante D von (2 a) durch die Forderung festge- 
legt, da die Schubkraft am Scheibenaufenrand mit der Schubspannung 1t(r,) = t, gleich sein 
mu8 der Massenbeschleunigung der Randmasse m,: 


Ta Vo =AM,T, : (7) 


Die Integrationskonstante E von (2b) hingegen wird durch die Befestigung der Scheibe auf der 
Welle bestimmt: 
v(T9) = Uo — 0 e (8) 


3. Die sogenannten klassischen Profile. Wie fiir den stationadren rotativen Radialdehnungs- 
zustand lassen diese Profile auch fiir den instationadren azimutalen Verschiebungszustand geschlos- 
sene bzw. leicht tabellierbare Lésungen zu. 

a) Scheiben gleicher Dicke y = b. 

Gleichung (2a) vereinfacht sich zu 


dt Tt ; 
Potty thAr=0, 


wozu man mit D als Integrationskonstante die Lésung 


_D_ 4A, 
nae (9) 
1 Der allgemeinere Fall verdnderlichen Schaufelquerschnittes sowie die Beriicksichtigung 
Bindedrahte und eines Deckbandes soll in einem anderen Zusammenhange spater erdrtert werden. ier 3 
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Ilo 


findet, wahrend (2b) fiir die Azimutalverschiebung v mit E als weiterer Integrationskonstanten 
D wa 


23 a8 3 
v Fee ses Er (10) 
ergibt. Aus (10) erhalt man leicht 
dy aD 3 UA, dv D 3d 
Poe ae 1 eS de Gr AG 


und durch Einfiihrung in (3) 
D sua lee wa fleas 


Grin AG ''T Gr .4G. ‘Cae 


wodurch die erwahnte Kontrolle der Ergebnisse (9) und (10) bewerkstelligt ist. Die Anpassung von 
(9) und (10) an die Randbedingungen (7) und (8) liefert fiir D und E die Werte 


Daan(e +), (Ila) Bey jom 
0,078 
Uri he ee Le Bar 
. rena ( = ae i (11b) 
tkpret 4 
Dies ergibt nun mit (9) und (10) die endgiiltigen Lésungen ea 
fe ces 500|- G5 
Ho) =A + ija ae (2) 
_ Ar’ (me ova\ Lt 1 wa pe 
BPE) tee. ay 2 
: ; 100- 67 
Bei fehlender Ringmasse m, = 0 vereinfachen sich (12) und (13) zu 
A 0 
Ropers 
<F) = “f (3 EG r) . (12a) 
Au r4 
v(r) = £4 (2 — 7) (3 : r): (13a) He 
Man erkennt itbrigens aus (12a), daBt (r =r,) =T, = Vist, wie 


es bei fehlender Ringmasse m, = 0 gemaf (7) sein muB. Ferner 
folgt nach (5) mit y = b Abb. 2.. Die Schubspannung 7, die Azimu- 


talverschiebung v und der Azimutalwinkel } 
fiir Scheiben gleicher Dicke. Unten ist d 
Fa 3 4 gleicher Dicke en ist der 
T=22 MT cin CE U b te ae 9) ? (14) Verlauf der halben Profildicke y/2 = b/2 
angegeben. 
und man hat zur Kontrolle nach (6) mit (12) fir r =r, und y = b 


3 2 
M, = |2a%m,r2 + eb (re—ri)] A= dang 2 (Ts oe a ra od ‘ (15) 
ro 


Die Gleichungen (12) und (13) zeigen iiberdies, daB sich die Funktionen 1(r) und v(r) erwartungs- 
gem4B proportional mit 1 vergréBern. Dies gilt auch fiir die anderen Profile. 

Zahlenbeispiel 1: Es werde gewahlt: r7 = 10cm, r, = 50cm, w =7,85-10-% kpsek?/em‘. 
n = 4875 U/min. Diese Drehschnelle soll in t = 0,1 sek durch einen axialen DrehstoB vernichtet 
werden. Mit w = n2/30 = 510 sek und A = w/t = 5100 sek erhalt man A = 4-10 kp/cm*. 
SchlieBlich erhalt man mit G = 850000 kp/cm? noch A u/8 G = 0,589 - 10-8 cm. 

Es sei kein Schaufelkranz vorgesehen, somit m, = 0. Damit erhalt man aus den Gleichungen (12 a) 
und (13a) die Zahlenwerte von t(r) und v(r) und aus (1) diejenigen von #9. In Abb. 2 ist der 
Verlauf der Funktionen und unten der Profilverlauf wiedergegeben. Die im Gegensatze zu den 
folgenden Beispielen mit verjiingten Profilen groBen Werte von t in Wellennahe sind durch die 
groBen duBeren Scheibenmassen bedingt. 

b) Scheiben mit hyperbolischen Profilen; alsoy =br-”. Damit ergibt (2a) die Differen- 
tialgleichung 


dt 2— 
eo T 


Eyek war 0. 
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die mit D als Integrationskonstanten zur Lésung 
D wa 
rn A 


i (16) 


i 


fiihrt; sie geht fiir Scheiben gleicher Dicke mit n = 0 sofort in (9) tiber. Ebenso erhalt man nach 
(16) aus (2b) mit E als neuer Integrationskonstanten 


D if bar 
poy rie eo ene ae se 


 — 


woraus mit n = 0 sofort wieder das Ergebnis (10) fiir Scheiben gleicher Dicke folgt. Aus (17) ge- 
winnt man durch Ableitung 

dv Di—n 1 3dr? E d*v DI—n 3dr 

a. Gi—nr—s 2(4—n)6 | ¥ z. = 


Hiermit ist Gleichung (3) identisch erfiullt (wegen : : ab Dee Se ; 


Mit (16) und (17) erhalt man aus (7) und (8) 


1— n\ 
D=iAr8 (e a (18a) 
be Ars Tag a YET ee wre 
Sree 1 = 2 2(4—n)G ash 
und somit 
(r) =A ( fig lee ur 
A ae b EE a ri—n = 4—n]’ (19) 
uy EAS ite Naess: ( r 1 bar 
0) = (2 —n)G (" | pn ay — 2(4—n)G (7? —r6). (20) 
Mit n = 0 gehen (18a, b) in (11a, b) und (19) bzw. (20) in (12) baw. (13) iiber. 
Fur n = 2 folgt wegen i 
r i 
2—n 1—n 
lim —° 5 : =rln— 
n—>2 =F 
sofort 
Ar> /ma 
on) a (5 - a3) 7m = — ar(e AD (21) 
a 0 
Ferner erhalt man aus (5) 
<< 27 ae) 
T One yo Ore oe") (22) 
und aus (6) wegen 
wnafte +) gt pl 
—n ae 4—n 
gemadB (19) ; 
2 — 7 
My=2ahm e+ ge whb(e—"— ri") = 2a rhb ena |(™ 4 Hg 7a 3 AGM 
b 4—n ak 4 


Hier erweist man durch Ausmultiplizieren der rechten Seite leicht die Identitat beider Gleich 


seiten. Fir n = 0 folgt sofort Gleichung (15). Pee 
Fir m, = 0 erhalt man schlieBlich aus (19) und (20) 
TA ote 
u(r) = far = ag r| , (19a) 
o(r) = a Ey r 1 rr 
@—n)G |2—n\2=" — Fin] —7 9 (20a) 


wahrend (21) fiir n = 2 ergibt 


_ Apr 2 r 9 
ae 2 ra ing il (21a) 
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Zahlenbeispiel 2: In den folgenden Unterfallen 1), 2), 3) i i j 
: , 2), 3) ist der Vorfaktor b =——iPiraa i 
so gewahlt worden, da die Scheibendicke stets yy = 10 cee ee te 


1) n =1;6, = 100 cm?. Im ibrigen seien die Annahmen wie im Beispi a 
: 3 ispiell. Manerhalt y, = 2 
- (mit r, = 50cm). Aus (19a und (20a) folgt speziell a" " 


Am (508 a 
t ee ia [oS ae ee Lh r r 
(r) 3 ( : r und U(r) 3G so (75 — i — (r?— 100)|. 
2) n = 2; by = 1000 cm’; y, = 0,4cm. Hier wird v(r) nach (21a) bestimmt. 

3) n = 3; bs = 10 000 cm‘; y, = 0,08 cm. In diesem Falle werden t(r) und v(r) wieder nach 


(19 a) und (20a) bestimmt. Der erste Term in der eckigen Klammer von (20a) bleibt trotz der Be- 
ziehung n > 2 positiv. 


Abb. 3. Die Schubspannung 7, die Azimutalverschiebung v 

und der Azimutalwinkel # fiir hyperbolische Scheiben mit den 

Exponenten a) n = 1, b) n = 2 und c) n= 3, sowie die halbe 
Profildicke. 


9 0 2 3% 40 Soemr0 0 24 3% 40 soomr 0 0 2 3% YO boOemr 


In Abb. 3a) b) c) sind die Ergebnisse dargestellt. Die Werte der Schubspannungen Tt, der Ver- 
schiebung v und des Winkels # nehmen mit dem Betrag des Exponenten n ab 1. 

Die Randschubspannung T, kann betrachtliche Werte annehmen, wenn ein Kranz mit Schaufeln 
vorhanden ist. Sind z. B. in Abstanden von je 5cm Schaufeln konstanten Querschnittes von 
| = 25 cm Lange und 1 kp Gewicht vorgesehen, so folgt — bei Beibehaltung aller tbrigen Konstan- 
ten des Beispieles 1 — mit G, = 62,8kp, r,=55cm, f, = 100 cm? und r, = 67,5cm aus (4) 
m, = 0,001414 kp sek? cm~*. Damit ergibt sich aus (19) fir n=1 und r=r, der Wert: 
T, = 180,5 kp cm. Fiir n = 2 erhalt man (wegen des kleinen Wertes y, = 0,4 cm) sogar den fiinf- 
fachen Betrag, was als Kontrolle dienen mag. 

Die Schubspannung t(r) kann (man vgl. Abb. 3c) u. U. einen Extremwert annehmen. Seine 
Lage r bestimmt man aus der durch Differentiation von (19) folgenden Gleichung: 


met (4 —n) (n—2) tur, |n—2—2 (SY | = 0. (23) 


r 


1 Man kann zeigen, daB auch fiir n = 4 die Ausdriicke (19) und (20) endlich bleiben. Von der Wiedergabe 
dieser Ergebnisse kann aber hier abgesehen werden, da man in der Praxis kaum tiber n = 2 hinausgehen wird. 
Man vgl. hierzu J. Malkin, Festigkeitsberechnung rotierender Scheiben, S. 57, Berlin 1935. 

5 * 
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Fir m, = 0 folgt aus (23) die einfachere Gleichung: 
n—2 


Ta 


die n = 3 und r = r,/2 (Abb. 3c) in der Tat befriedigen. Fiir r = ry = 1,/5 folgt aus (23a) 
5” — 312,5 (n — 2) = 0 
mit nw 2,0929 als Naherungslésung. Bei diesem Profil beginnt somit t(r) bei m, =0 mit horizontaler 


Tangente in r = 1». 


ee Sars 5 
c) Scheiben gleicher Dehnungsfestigkeit o 9, also y(r) = be 23%. = ipa eeee 
Die Einfiihrung von y(r) in (2a) ergibt nach kurzer Rechnung fiir 7(r) 
dt 2 p 
ais (7 —2kr|e + war=0. 
Integration ergibt + 
D a ee 24) 
Fiihrt man nun (24) in (2b) ein, so st6®t man auf das Integral ” 
ek r? >. ek r* Pes kr? 
[a-Si + te Us 
worin li den Integral-Logarithmus bezeichnet. SchlieBlich folgt 
as D kr? Dk . foert BA 1 > aly 
I) eae +e yrlile ) +E ap (brine gz) + Er. (25) 
Wegen 
1 dy Fae 1 
icc ues emi cokes 
erhalt man aus 
De De Db gs eed iS 
Se OT te lie )+Gop|ka +1an +55 +E 
zunachst 
1 dy 1\ /dv Ce 2D 5 a Da a weet Nae wa 
ee | 7) zt Og On hei es eae ieee 
Mit 
Oe 23 3D) ge wa 1 BA 1 
dr? Gr Corr G 2kr Od Si ee 


ergibt die Einfihrung in (3) durch Tilgung aller Glieder die erwiinschte Kontrolle. 
Die Anpassung von (24) und (25) an die Randbedingungen (7), (8) ergibt fiir D und E nach 
langeren Rechnungen 


Dm dree tM — talk +3) mit y,=be "a, (26a) 
Dio tea 4 1 
E 2 Ee — thee | 2 ia (kim ry — 5h). (26b) 


Damit erhalt man aus (24) und (25) 
— 4 |(re\’ —k(2—12) [mara iu 1\} bh 1 ates Bia 
el wet tts) ste aie 
: meee Pay," Tere Ones ae ets bAlTr r LS fr. 1 
v(r) x6 | . ae 0 kr [li (e 0) li (e +e lee tae(a—s]]- (28) 
Hierin ist noch D aus (26a) einzusetzen. 
Ferner findet man nach (5) wegen 3 


fe ek dr = (1 + kr?) Fr. 


+ Man vel. hierzu W. Meyer zur Capellen, Integraltafeln, S. 228, Nr. 1.3, Berlin, Gottingen, Heidelberg 1950. 
* Man vgl. hierzu W. Meyer zur Capellen, Integraltafeln, S. 228, Nr. 2.3 und S. 289, Nr. 3.1. 
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den Wert ; 
T =22m,r3 +"5 [e—*0 (1 + kr?) —e-*G (1 + kr)). (29) 


Geht man nun mit (27) und (29) in (6) ein, so erhalt man 


M, = 20m, + SEPA Te #8 (1 4 kt) — oP + hd] 


= 2 —kr Ta\* 2 re 72 Mata i r2 va 1 1 \ 
navenifafene- alread tl 


0 
Hiermit ist die Ideutitat erwiesen und (27) mit (29) kontrolliert. 
Zahlenbeispiel 3: Wie bei ge 1 sei m = 510 sek; ferner werde oy = 1000 kp/cm? 


510? - 7,85 . 10-6 
gewahlt. Damit folgt k = OE =10-° cm. Die Konstante 6 mége nun so ermittelt 
werden, daB yy) = 10 cm sei; damit findet man b = 10- e%! — v' 
11,052 cm. Dabei ergab sich D gemif (26a) mit m, = 0 zu mat oes 
D =a eG (F +7) 140 : 
4-10? - 2500 3,9 4 i 
SE or eT ae 2,5 35 10-3 = — 5746 kp. 100 ae G005 
Bemerkenswert ist in der zugehérigen Abb. 4 der fast geradlinige ll Ee 4 é es 
Anstieg der Azimutalverschiebung v(r). a G6 7\_ He ea 
: 40} x} g002 
d) Konische Scheiben? also y =a ( —3] (Abb. 5). Mit c 
R 20+ 02 goof 
dieser Profilkontur nimmt (2a) die Form an | 
0 0 2 DH 4 siemr 
dt 2 1 
dr ie . R=)" bAur=0. 
Hierzu findet man mit D als Integrationskonstante die Lésung | 
SJ 
D _ A 22 R—4r 
ee Ras) 00. Re oo) 
Mit (30) gibt (2b) fiir v(r) nach langerer Rechnung und mit Sls 


E als neue Integrationskonstante die Lésung 


v(t) =— a5 ee 42 (R—r) rin(A 


Abb. 4. Die Schubspannung 1, die Azimu- 
a talverschiebung v und der Azimutalwinkel 3 
fiir Scheiben gleicher Dehnungsfestigkeit, so- 

wie die halbe Profildicke. 


ti [2r?—r? R—r R?ln(R—r)|+ Er. (31) 


Aus (31) erhalt man durch doppelte Ableitungen 


dv D {[R?—Fr R—r R 
F=TR| 2 +2—In(~ )+z 
Ap 2 a r R? 
—mo|o"* 27 R— R?In(R ale, 
dv D 


Abb. 5. Axialschnitt durch eine konische 
R? oe Scheibe mit den der Berechnung zugrunde 


ae 0G | + (R a 7 gelegten Abmessungen. 


Die Einfiihrung dieser Ausdriicke in (3) befriedigt diese Gleichung. Weiter hat man durch An- 
passung von (30) und (31) an die Randbedingungen (7) und (8) 


D=1/"s tier! oe (32a) 
Bre aa Ieee <=" In a(S aE JH [218 —1y R—R*In(R—¥)]- (2b) 


1 Man vel. C. B. Biezeno und R. iat Technische Dynamik, Bd. 2, S. 19, Abb. 18. Die dortige Be- 
zeichnung ist hier im wesentlichen iibernommen worden. 


R 
P ~ex | re t r(R—r) ter 
12 = 
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Mit (32a, b) erhalt man aus (30) und (31) 


A ala Ta BS . 
er) = mala aren fs ue 2(5R—4r,)—2(5R tn). (33) 
D R—1H)? (R—r)? (R—rp) an R—1r aa 
Rs en| ee 2r ce ee Bae R= to 
A 2 nals 
+967 —208 =) +RE—7)— ei ae. (34) 
Hierin ist noch fiir D der in (32a) ermittelte Wert einzusetzen. 
Ferner ergibt sich nach Gleichung (5) als axiales Massentragheitsmoment 
T=22 m7 > RIR (r2 — r§) —4(r3 — r)]. (35) 


Beachtet man zur weiteren Kontipite Gleichung (6), so erhalt man die Identitat 


A 
Thplem’ pa vem M, =20m,AA +77" zd R(x — r4) — 4 (r3 — 7?) ] 
90 
a BS “ad = 2218 y¥_% = 201 = (R—M) 
Or # g0o7 
Cor 2 0.006 A mr, R a A 
soe [Mate RE [Te (5 R= 47) 12 (S Re 
50 40 2008 Ral er 20 72 ( ) 0 ( S 


Zahlenbeispiel 4: Es werde wieder ry = 10cm, r, = 50cm, 
0003 = yy = 10cm und y, = 2 cm gewahlt. Dann folgt gemaB Abb. 5 
0002 a=—12cm, R=60cm; weiter sei Ay =4-10-*kp/em* und. 
G = 0,85-10° kp/em?, m, = 0. Fir D ergab sich nach (32a) 
D =1250000 kpem. In Abb. 6 sind die Funktionen 1(r), v(r) 
0 0 2 wo Ww wmp und Pr) dargestellt. 


e) Erstes Exponentialprofil von Malkin, namlich 
y =be—'""*, Die Einfiihrung dieser Profilkurve in Gleichung (2a) 
liefert die Differentialgleichung 


Blco 
dt 2 
ween beget ay 1/3 — 
Ble zOre\c + par = us 
Blea 
die die Lé it 
Abb. 6. Die Schubspannung 1, die Azimutal- e d € sobee ha 
verschiebung v und der Azimutalwinkel 3 3 e 
fiir konische Scheiben, sowie die halbe Pro- 4/3 Amey le} — 2/3 —2 
fildicke. u(r) = > T ag SAP TT a 28 re Are eo) 


Mit (36) findet man aus (2b) nach langerer Rechnung fiir die azimutale Verschiebung v(r) mit E 
als neuer Integrationskonstanten 


D ere ower s 
=| = dr + pate ($ perse —3 rer) 4 Br. (37) 

Das in (37) noch enthaltene unbestimmte Integral kann durch die Substitution 

: y3/2 
— pit 9 
somit 
3 ate d 2 

dr = > BA du , 6 ele ees ae - =a pare (38) 


transformiert werden in 


eb rtl® ge : 3 wes : 3 
[erm [eed peten 3 om [Gee bora, 09 


¢ 
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Man kennt das unbestimmte Integral J,; es soll in u ausgedriickt werden; } 
= 4 uz 1 u2 1 
h= 3 fe du— 6 (+5): (40) 
Durch Ableitung von (37) nach r, wobei man (38) und (40) beachten muB, folgt 


Bead). (14 1 1 
Bese = 3/2 vu? uw dJ, du Sia 5 
= 2 cP {ls [e PaO rks (2 | ) hale, Z| Tens ee tt ee 
worin 
dJ, du 2 9 321 grass 
Stas dp Tage Cree oF 


ist. Nochmalige Ableitung dieser Gleichung nach r ergibt wegen 
d-[.. 1 it eg me 1\ du 
dr e a (2 4 =| re ae 
mit Beachtung von (38) nach langerer Rechnung 


GED cc hb Au (5 1 
mo ee ee 5/3 = t Ro(a Ber 18 er, 


: a ld 1\ /dv 
wahrend man fir 5 ~ . (F — 7) denselben Ausdruck mit umgekehrtem Vorzeichen ver- 


mehrt um —Ayr/G erhalt, so daB man erkennt, daf die Gleichung (3) befriedigt und somit die 
erstrebte Kontrolle erhalten worden ist. Die Anpassung von (36) und (37) an die Randbedingun- 
gen (7) und (8) ergibt weiterhin 


3 3 
(pay, es re (B2 r3/3. 4 2B rt/3 +4 2) a | : (41a) 
Biel BD 3A 3 
Be eg Pi) aca iS peg pty he ro) , (41b) 


worin noch D aus (41a) einzusetzen und J,(u9) = J,(61/? r?/*) ein reiner Zahlenwert ist. Mit (41a, b) 
erhalt man aus (36) und (37) 


u(r) =A {2 Ta 1 aprais 7. E (6278/3 +2 Bri + 2) e—P(ra'®—r*)__ (9272/3 4.9 Br—28 42 rh ; 


| (42) 
u(r) = = = B32 r [Jy(u) — Jy(up)] + ice r Fi B? (r2/3— 72/3) — 3 B (r—2/8 —r— 2/8) — (1-2 — rm) 
(43) 


In (43) ist noch D aus (41a) einzusetzen. 
Ferner folgt aus (5)!: 


Seb 
T = 2m, r3— > [ (6 rh? + 2B re + 2) oP — (Brg + 2B rg? +2) 09H]. (Ad) 


Geht man mit (42) fiir r = ry) und yy = b e 0° und mit (44) in (6) ein, so erhalt man als weitere 
Kontrolle die Identitat 


Mz =Afamm, re —2ocbH [pera +28 1 42) e-Pd — (prgM + 2B HP +2) PN 
—f rr! ats r 3 1 rm i Rell aad 
=2aribe Eee eae 2s + 2g hae 
0 0 
SG ee Tet | 


Zu bemerken ist noch, daf in (44) der Ausdruck in der eckigen Klammer des zweiten Termes 
negativ ist, so daB der zweite Summand in (44) positiv ist. Abnlich den hyperbolischen Profilen 
wird die Konstante f vom Verhiltnis der Scheibendicken yp und y, gema8 der Formel? 


ey a 


1 Man vel. W. Meyer zur Capellen, Integr Itafeln, S. 228, Nr. 4.2.4. Dort ist allerdings irrtiimlich ein neun- 
mal zu groBer Ausdruck angegeben. Fiir das bei der Berechnung von T in (44) auftretende Integral / w? e~” dw 
vel. man W. Meyer zur Capellen a. a. O. S. 224 Nr. 4.1.1.2 mit k = — 1. 

2 Man vel. J. Malkin a. a. O. S. 71 Gl. (86). 
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bestimmt, die sich sofort aus der Gleichung der Profilkontur ergibt, wahrend sich die Konstante b 
nach Wahl von f durch die Forderung einer bestimmten Profildicke z. B. yo ermitteln 1aBt. 


Zahlenbeispiel 5: Wie bei den friitheren Beispielen sei rp = yop = 10 cm, r, = 50 cm, y,® 2em; 
Ap =4,10-2 kp/cm4, G = 0,85 - 10° kp/em?; dann erhalt man aus (45) 


In 5 Ee 1,60944 erp 0,01 nt b = Yo et 0,01 104/38 = 12,40407 cm. 


Pe 5047/3 — 1043 162,056 


Damit findet man die in der zweiten Spalte der Tabelle 1 (linke Halfte) angegebenen Werte der 
Funktion y(r). Ferner erhalt man fiir den Ausdruck in der runden Klammer von (41a) den Zahlen- 
wert 


A, = B2r33 + 2B r43 +2 = 9,076964 


und damit aus (41a) mit m, = 0, D = — 43160,9638 kp. 

Mit m, = 0 verschwindet in (42) der erste Summand. Bezeichnet man dann innerhalb der eckigen 
Klammer die mit der e-Potenz behaftete Funktion in (42) mit t,(r) und den letzten Summanden 
mit 7,(r), so findet man ihre Zahlenwerte in der dritten und vierten Spalte, ihre Summe t(r) aber — 
multipliziert mit dem Faktor — 3 u 4/463 =— 30000 kp — in der fiinften Spalte der linken Halfte 
der Tabelle 1. 

Analog enthalt die rechte Halfte dieser Tabelle in der ersten Spalte die gemaB (38) bestimmten 
Werte von u, in der zweiten und dritten Spalte wieder die beiden Teilsummanden v,(r) und »,(r) 
von (43), wahrend die vierte Spalte deren Summe v(r) und schlieBlich die fiinfte Spalte den Tor- 
sionswinkel ? nach (1) enthalt.1 


Tabelle 1 


T(r) in (42) T,(r) in (42) t(r) kp/em? v,(r) in (43) | v,(r) in (43) v(r) em 


0,01 7846 | — 0,024 773 | 207,825 0 0 0 0 


20 7,207 0,00 6191 | —0,008451 | 67,818 0,736.8 | — 0,005 259 | 0,007271)| 0,002612| 1,006 
30 4,883 0,00 4060 | — 0,005 259 | 35,967 0,9655 | —0,010018 | 0,013 761 | 0,003 744) 1,248 
40 3,158 0,00 3532 | —0,004130| 17,934 1,1696 | — 0,014 871 | 0,020233 | 0,005361)| 1,340 
50 1,966 0,003631 | — 0,003 631| —0 1,3572 | = 0,019 983 | 0,026 810 | 0,006827) 1,365 


Tabelle 2 


T(r) in (52) t,(r) in (52) 1(r) kp/cm? u v,(r) in (53) v,(r) in (53) v(r) cm d+ 104 


3,337002 | — 3,824 607| 119,834 
5,058 1,649030 | —1,870984| 54,547 
2,853 1,299388 | — 1,435 674} 33,494 
1,714 1,216814| —1,289451| 17,851 
1,072 1,244628 | — 1,244 631 0 


1,1604 0 ~10 0 0 

1,8420  —0,009355 0,010 691 0,001336) 0,668 
2,4137 | — 0,019 123 | 0,021757 0,002 634| 0,878 
2,9240 | — 0,029 646 | 0,033512  0,003866| 0,966 
3,3930 | —0,040978  0,045962  0,004984| 0,997 


Es ist also 


A 4/3__ ,4/8 
Tete) = pe FU au y 4 
t(r) =— 75 [zx(r) + 7,(r)] . 
1,(r) = —(f? P23 +28 r—28 4 2 r), 


a(t) => GBM lu) — Sila) 
v,(7) == ae: = Tr - p? oe = 72/3) ae B (fe 25 ee = 2/3) ie (r—2 Sr —2) v(r) = v,(r) oe v,(r) : 


1 Zur Auswertung der in J, enthaltenden Funktion Je du konnten mit Vorteil die Funktionentaf 
Jahnke-Emde, 5.106, 2. Aufl. Leipzig und Berlin 1933, fiir den Integral-Logarithmus in (28) Cation 
beispiel 3), sowie in (50) (Zahlenbeispiel 6) ebendort S. 83 bis 86, aber auch F. Tolke, Praktische Funktionen- 
lehre, Ss. 248 bis 297, Berlin 1943 verwendet werden. Fiir die oftmals gebrauchten Funktionen e* bzw.e-* wurden _| 
die »»Fiinfstelligen Tafeln“ von K. Hayashi, Berlin 1944, und schlieBlich fiir die mehrfach auftretende Funk- 
tion x?/3 die ,,Fiinfstelligen Funktionstafeln“, von K. Hayashi S. 106—110, Berlin 1930 benutzt. 
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Die so gewonnenen Endergebnisse sind in Abb. 7 durch Schaulinien mittels der beigefiigten Mab- 
stabe dargestellt worden. 

f) Zweites Exponentialprofil von Malkin, namlich Von b eb rele 
Die Einfiihrung dieser Profilkurve in Gleichung (2a) ergibt die Differentialgleichung 


ae p Ariz + war=0, 


die mit D als Integrationskonstante nach langerer Rechnung die Lésung 


eae A 
o(r) = ofr 4 SAK ze (Br + 5 A472! + 20 6% + 60 f%r—-2 4 1208-4 412012) (46) 
liefert. Fuhrt man nun ie in (2b) ein, so erhalt man mit E als cipbm? yoo es 
weiterer Integrationskonstanten 200+ — 
D eae 3Au (3 15 780|- 
ed Se AG eo aise AALS) 3 
v(r) 2 (ee dt 36K ra ee aeeg tT + 20 6 rlnr a 
— 90 6? 3 — 90 B 13 — 60 r- ) +Er. (4) er 
reo} 42) — 9008 
Das in (47) noch enthaltene unbestimmte Integral kann nol 40 18 | 005 
durch die Substitution / é 
: 80+ 98 3.004 
/2 / 
San (3) ? 60} 96 ; 1 g093 
somit aa a G.002 
s | 
dr = = wil? du , ee = 2 etaa 3 u== BrP | (48) ar et ill 
di " 010 2 a0 WO tom 


transformiert werden in 


B r2/3 a 
e 3 e 3 
[Sea-ge [saa pen, oe 
worin! me 
eee! eu 4 2 a i or 50 
Jy 2 tute) + Lier. (50) 4 
Beachtet man (49) und (50), so erhalt man durch Ableitung ‘ik eee Pe ca oy 
von (47) nach r verschiebung v und der Azimutalwinkel 9 
* D Ae Se i gemaB beni 1 oF rotleiee ee gen 
dv 3 e A profil sowie die halbe Profildicke. 
er ree|—sa tute) + HO tae 
ms wee (ze 7418. 4.22 Bt 2/3 + 20 68 (1 + nr) — 30 621% + 308 4 + 60 4) iE. 


worin u durch (48) erklart ist. Nochmalige Ableitung dieser Gleichung nach r und entsprechende 
Zusammenziehung der Glieder ergibt 


dv D 2 a 
UL) a ial 4-8 os 2, 8 
7 a ae (zB? ; 


a 3A pb Guz ris += tr M8 + 20 fF 7-1 +20 ptr 58 — 408 78 — 1207-4), 


26 Be 


wogegen man fiir ge a + rales Bese: denselben Ausdruck, jedoch mit verkehrten Vorzeichen 
und auSerdem noch vermehrt um — A fer/G erhalt, so da man nach Einfiihrung dieser Aus- 
driicke in (3) erkennt, daB diese Gleichung identisch befriedigt ist. 

Unterwirft man die Lésung (46) der Randbedingung (7), so gewinnt man 


D=A ee = re (B5 29/3 4 5 Br! 4 20 6% r2 + 60 B® r4l3 + 120 Br + 120)-eF'a"|. (51a) 


1 Ww, Meyer zur Capellen, a. a. O. S, 225, Nr. 2.1.4, mit der Erklarung von E i(x) auf S. 289 Nr. 3.1. 
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Unterwirft man analog die Lésung (47) der Randbedingung (8), so ergibt sich 
‘3 15 Lpyge he — 
i = ie TA Bé ri? + = B4r2/2 + 20 69In ry—90 f?r, *—90 Br, 4? —60r, ) 
(51b) 


Fiihrt man (51a) in (46) ein, so erhalt man fiir die Schubspannungsfunktion 1(r) 


—_ B (72/3 — 2/3 
o(r) =A [Mere 5 oh spe [pr (BETH +5 BUR + 20872 + GO prg2 +120 13° + 120) ¢ p (r3/9— 7913) 


— (5 r4l3 4.5 B4 r2/3 + 20 B8 + 60 2-2/3 + 120 B 43 + 120 ry] . (52) 
Die Einfiihrung von (51b) in (47) ergibt mit Beachtung von (49) 


15 ie 
ofr) = Gy Br Lela) — Jaleo] + ag po? [a BE (8 — 8) +5 BEC — 19) + 20 Bm 


— 90 B2 (r—2/3 — r— 2/3) — 90 B (r— 48 — 413) — 60 (r-2 — cu) Rone os 


Hierin ist noch fiir D der Wert (51a) einzusetzen, was der Ubersichtlichkeit wegen unterbleiben 
mége. u bzw. up bedeutet den Wert der Funktion u(r) fiir r bzw. ro. 
Ferner erhalt man fiir das Massentragheitsmoment gema8 Gleichung (5) 


326 
T= 2am, rt} + 2ayb | re-0" dr = 2amyre + “ye [ whem du 


Ta 


To 
To 


Man findet fiir das unbestimmte Integral 
fu e—” du = — (ub + 5 ut + 20 u? + 60 u? + 120 u + 120) e—*. 
Somit folgt mit Beachtung von (48) 
T=2am, nae G 10/8 4. 5 B4 18/3 + 20 63 r2 + 60 B2 r4/3 + 120 B 12/3 + 120) e749 
— (5 r}9I3 + 5 B4 18/3 + 20 6% 12 + 60 f2 r4/3 + 120 B r2/ + 120) eas : (84) 


Fijhrt man nun in (52) r = ry ein und geht mit dem erhaltenen Ergebnis und mit y, = be? ce 
in Gleichung (6) ein, so erhalt man zur Kontrolle die Identitat 


i Rea P mm, 73 — 22H Bs p10/3 + 5 BA 78/3 + 20 6% r2 + 60 B% r4l3 4 120 6 r2/ + 120) e-P72 


ise 
— (B5 1} +. 5 B4 18/3 + 20 B3 12 + 60 B2 rl? + 120 Br2/3 + 120) e—? sal 
— Brel Mars 32/3 3u fl 
=2nr2be? ea te Oe |e (BP riO + 5 Btn’? + 20 8812 + 60 fF ri 
+ 120 B 72/3 + 120) eB (2/8 — 1313) 
— (5 r4l3 + 5 B4 12/3 + 20 6% + 60 B2r,—2/3 + 120 B ry—4/3 + 120 ro) 


Zu bemerken ist wieder, daf in (54) der Ausdruck in der eckigen Klammer des zweiten Termes 
negativ, der zweite Summand in (54) also positiv ist. Die Konstante f wird wieder wie beim ersten 
Exponentialprofil yom Verhaltnis der Scheibendicken am Innen- und Auenhalbmesser der Scheibe 
gemaf der Formel? 

In 22 


Bf 
=. * _, 55 
B Pls eis ( ) 
abhangen, wahrend der Vorfaktor b dann wieder z. B. durch die Profildicke y, bestimmt ist. 


Zahlenbeispiel 6: Wie beim Zahlenbeispiel 5 sei ry = yy = 10 cm, r, = 50cm, aber y,~ lem, 


ferner Aw =4-10-2kp/cm* und G = 0,85 -10°kp/em?. Dann erhalt man zunidchst aus (55) 


In 5 
i= ie 1073 = 0,258; es wurde gewahlt: 6, = 0,25. Fiir den Vorfaktor b ergibt sich dann 


* Man vgl. J. Malkin a. a. O. S. 71, Gl. (86). 
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= 10- orate 10° = 31,91202 cm. Damit ergeben sich dann die in der zweiten Spalte der Tabelle 2 
(linker Teil) angegebenen Werte der Funktion y(r). Ferner findet man fiir den Ausdruck in der 
runden Klammer von (51a) den Zahlenwert 


A, = Br? + 5 Bt 183 + 20 63 r2 + 60 B? ril3 + 120 6 r2/3 + 120 = 3111,569438 
und damit aus (51a) mit m, = 0 den Wert D = — 25693,8966 kp. 


Wie im Falle des Zahlenbeispiels 9 hat man von (52) und (53) die in der Tabelle 2 ] ie I 
der Tabelle 1 berechneten Teilfunktionen a aes sae 


afr) = PGP), 


e : 3 ud 
T(r) = — (65 r49 + 5 Br? + 20 6% + 60 p2r—28 ptr) =— Ze [(r) + TAr)]- 
+ 120 Br— 4 + 120 r-2) , 
3 D 


(7) => FBP [Ja(u) — Jaluy)] | 
v,(r) = sou s ns B(Ae— rs) = B? (77/3 — 2/8) +. 20 63 In ze pr) = 0,(r) + »(r) . 


— 90 2 (r—2/3 — c 2/3) 90 B (r— 48 — re 4/3) 60 (r-? — ra) , 


Die Endergebnisse der Tabelle 2 sind in Abb. 8 dargestellt. 


Thper’ yx" vem 
Die hier behandelten Malkinschen Exponentialprofile sind 120 


Spezialfalle allgemeiner Exponentialprofile von der Gleichung 100|- 40 0,005 

yre= be Ft", (56) 80\- 08 004 
die fiir n = 2 beibesonderen Randbedingungen wieder auf die in 60) 96 ies 
Abschnitt c) behandelten Scheiben gleicher Dehnungsfestigkeit 40} o# 0,002 


zurickfihrt. Fir allgemeinere Randbedingungen haben diesen 
Fall bereits A. Fischer! und R. Gran Olsson?, letzterer auch die 
Falle fir n =1 und n =4 mittels der konfluenten hypergeo- 0 10 20 3 40 boemr 
metrischen Funktion behandelt und auch tabellarisch weitgehend 
unterbaut. Diese Profile und auch solche von parabolischer 
Querschnittsbegrenzung sollen in einer spateren Arbeit eingehend 
untersucht werden. S 


207 G2 g007 


4. Weitere Profile, insbesondere das Profil gleicher Schub- 


spannung 7). Man kann z. B. nach jenen Profilformen fragen, Blas 
deren Schubspannung dem einfachen Gesetz 
ms 57 Abb. 8. Die Schubspannung t, die Azimu- 
GS To sila ( ) talverschiebung v und der Azimutalwinkel } 


em&B Tabelle 2 fiir das zweite Exponential- 


genugt. Die Betrachtung von Abb. 2 35 6 und 7 zeigt insbe- profil, sowie die halbe Profildicke. 
sondere, daf man Profilkonturen, die sich nach auBen verjiingen, 

praktisch also brauchbar sind, fiir n < 0 zu erwarten haben wird. Fiihrt man den Ansatz (57) in 
(2a, b) ein, so erhalt man die beiden Differentialgleichungen 


dy 2+n Ap 1l—n — 

dr t ( r pee jy=0, (58a) 
dv WEE gin, 
| oe Gee (58b) 
die die folgenden Lésungen ergeben: 

ae uA pen 

Nee ee (2 —n) Tt ‘ (59a) 
1 

ae mOe n+1 : 59b 
Moet ae eval é (59b) 


Hierin bedeuten b und E die beiden Integrationskonstanten. 


1 A, Fischer, Z. 6st. Ing. u. Arch.-Ver. 74 (1922) S. 46. Man vgl. auch T. Suhara, Trans. Soc. mech. Engrs. 


Japan 3 (1937) Nr. 10, S. 1. 
2 R. Gran Olsson, Ing.-Arch. 8 (1937), S. 270 und S. 373. 
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Aus (59b) erkennt man, da auch fiir die Azimutalverschiebung v Ansatze gemaf (57) zu Profil- 
formen nach Art (59a) fiihren werden. Aus (57) gewinnt man eine konstante Schubspannung T) 
fiir n = 0, und damit erhilt man aus (59a) als Gleichung der Kontur y(r) fiir ein solches Profil 


yak 2 T% ae 28% (60) 
Wegen des Vorfaktors b/r? verjiingen sich also Profile konstanter Schubspannung, selbst wenn 
A = w* und Tt) = 6, ist, rascher als solche konstanter Dehnungsspannung oy. (Man vgl. Unter- 
abschnitt 3c.) 
La8t man die Massenkrafte auBer acht und denkt sich die Schubspannung Ty etwa auf stati- 
schem Wege durch ein Drehmoment erzeugt, so folgt aus (60) mit w/A = 0 
b 


ee 


somit also das von R. Grammel' fiir diesen Fall bereits gefundene hyperbolische Profil. Bei all- 
gemeineren Profilen, wie z.B. dem Grammelschen Profil? wird man, wie hier bemerkt werden mége, 
bei der Schub- und Verformungsberechnung i. a. geschlossene Formeln kaum erwarten und daher 
Quadraturen nicht umgehen kénnen. Das Verhalten dieses Profils und weiterer Profile, die etwa 
aus dem Grammelschen Umkehrproblem gewonnen worden sind?, gegeniiber axialen Drehmomenten 
soll Gegenstand einer spateren Untersuchung sein. 


(Eingegangen am 26. Januar 1960.) 


‘ Anschrift des Verfassers: Professor Dr.-Ing. Karl Karas, Darmstadt-Eberstadt, 
bs Carlo MierendorffstraBe 38. 


* Man vel. die in FuBnote 1 Seite 63 angegebene Literatur, insbesondere R. Grammel, Z. Mat 
(1925), S. 194, Abschn. 1: Scheibendrillung; und C. B. Biezeno und R. Grammel, ee a. 6. STG. : 
- ® R.Grammel, Ing.-Arch. 7 (1936), S. 137, insbesondere Gleichung (9). 
8 A, Held, Ing.-Arch. 10 (1939), S. 339. 
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x ZUR INFORMATION 

‘Dieser 2weite Band des Werkes behandelt die Schwingungen der Gebilde mit mebreren Freiheitsgraden, die 
sogenannten Koppelschwingungen. Das Buch will kein Handbuch, sondern ein Lehrbuch sein, und zwar 
eines, das sowohl Anfingern wie Fortgeschrittenen dient; deshalb ist der Stoff nach steigendem Schwierig- 
keitsgrad geordnet. Inalllen Teilen des Buches werden die Anwendungen der Schwingungen, vor allem im 
Maschinenwesen, im Bauwesen und bei Fabrzeugen im Auge behalten. Dariiber hinaus ist der zweite Teil 
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